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代数 方程 组 的 高 效 求解 是 计算 数学 中 一 个 基础 共性 且 十 分 关键 的 问题 。 本 书 根据 作者 近 几 年 在 代 
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在 很 多 科学 与 工程 计算 领域 ,都 会 遇 到 大 规模 代数 方程 组 的 求解 ,有 效 求解 大 规模 代 
数 方程 组 是 科学 和 工程 计算 中 的 重要 研究 内 容 。 

本 书 分 为 两 个 部 分 共 6 章 ,主要 内 容 概述 如 下 。 

第 一 部 分 第 1 一 3 章 首先 介绍 增 量 未 知 元 预 处 理 技术 和 近年 来 的 HSS 类 和 迭代 方法 ， 
然后 对 于 非 线 性 特征 值 问 题 , 基 于 二 阶 增 量 未 知 元 方法 ,提出 修正 的 Marder-Weitzner 方 
法 ;对 于 各 项 异性 反应 扩散 方程 ,基于 和 矩阵 分 块 ,提出 一 种 新 的 类 小 波 增 量 未知 元 
WIU) ,这 里 包含 了 作者 的 一 些 研究 成 果 。 

第 二 部 分 第 4 一 6 章 研 究 几 种 增 量 未 知 元 预 处 理 和 迭代 算法 。 第 4 章 对 非 埃 尔 米 特 正 
定 线性 代数 方程 组 的 数值 解法 进行 研究 ,基于 NSS 迭代 方法 ,提出 预 处理 NSS 方法 和 双 
参数 NSS 方法 ,并 对 方法 的 收敛 性 进行 细致 的 分 析 , 给 出 迭代 格式 中 迭代 参数 的 计算 方 
法 。 第 5 章 首 先 介 绍 复 对 称 线性 代数 方程 组 的 迭代 解法 ,提出 基于 SOR 迭代 的 复 对 称 线 
性 系统 的 MHSS 加 速 方 法 ,分 析 迭 代 格 式 的 收敛 性 ,并 以 数值 结果 证 明 方法 的 优越 性 。 
第 6 章 研究 非 线 性 代数 方程 组 的 数值 解法 ,基于 LHSS 方法 提出 Newton-LHSS 后 退 方 
法 并 分 析 方 法 的 全 局 收敛 性 ,对 于 一 类 特殊 的 弱 非 线性 代数 方程 组 ,提出 Picard-AHSS 
迭代 方法 IER HE AHSS-like 迭代 方法 以 及 Picard-MHSS 方法 ,分 析 方 法 的 收敛 性 ,并 给 
出 数值 结果 。 

本 书 多 数 为 作者 的 研究 成 果 , 在 相关 问题 中 参考 了 伍 渝 江 教授 近 几 年 的 研究 资料 ,在 
此 表示 感谢 。 
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1.1. 增 量 未 知 元 方法 与 微分 方程 数值 解 


许多 物理 运动 都 可 以 用 一 个 偏 微分 方程 的 定 解 问题 描述 ,例如 无 限 长 弦 的 自由 振动 
问题 可 归结 成 二 阶 双 曲 型 偏 微分 方程 的 初 值 问题 ,而 弦 对 平衡 位 置 的 偏 移 就 是 方程 的 解 . 
但 是 绝 大 多 数 偏 微分 方程 定 解 问题 的 解 不 能 用 公式 表达 ,有 时 即使 可 以 用 公式 表示 ,也 往 
往 过 于 复杂 ,只 能 用 各 种 近似 方法 计算 它 的 解 . 

有 限 差分 方法 是 求解 偏 微分 方程 定 解 问题 的 常用 近似 方法 之 一 ， 自从 Courant, 
Friedrichs, Lewy(1928 年 ) 提出 偏 微分 方程 的 有 限 差分 方法 以 来 ,已 有 很 多 人 对 此 做 了 
大 量 研 究 . 尤其 是 当代 高 性 能 电子 计算 机 和 数学 软件 的 高 度 发 达 , 促 使 偏 微分 方程 数值 
解 的 研究 蓬勃 发 展 ， 本 书 研究 的 是 偏 微分 方程 数值 解 领域 中 一 类 现代 数值 方法 一 一 基于 
有 限 差分 的 增 量 未 知 元 方法 ,并 用 于 求解 各 类 线性 和 非 线 性 方程 组 . 

对 于 许多 问题 的 数值 计算 ,例如 满 流 问题 ,经 常会 遇 到 大 规模 线性 方程 组 的 求解 它 
的 三 维 情 形 包 含 几 十 亿 的 未 知 元 ,如 果 将 所 有 的 未 知 元 不 加 区 分 , 即 对 于 占有 绝 大 部 分 能 
量 的 大 波长 部 分 和 能 量 很 小 的 小 波长 部 分 用 同样 的 方法 处 理 ,会 花费 同样 的 计算 时 间 ,这 
是 很 不 科学 ,也 是 很 不 合理 的 . 对 于 谱 方法 ,运用 Fourier 展开 ,可 以 将 这 两 部 分 自然 分 
JF. 但 对 于 有 限 差分 ,并 没有 这 种 好 的 性 质 . 为 了 克服 这 一 困难 ,基于 近似 惯性 流 形 理 
论 ,Temam 首次 提出 了 增 量 未 知 元 (Incremental Unknowns,IU) 方 法 . 该 方法 在 求解 线 
性 椭圆 问题 时 ,可 以 看 成 是 相应 系统 的 预 条 件 子 ,即将 求解 线性 系统 Au =b 转化 为 求解 
线性 系统 Aw =b ,其 中 A 二 S™A,Sb 二 STb. 这 里 ,S 表示 增 量 未 知 元 到 结 点 未 知 元 ( 按 字 
典 顺序 排列 ) 的 转换 矩阵 , 即 u 二 Sw. 已 有 大 量 事实 表明 ,IU 可 以 有 效 降 低 线 性 系统 系数 
矩阵 的 条 件数 ,并 且 对 于 线性 系统 的 求解 ,其 效率 与 多 重 网 格 方法 相当 . 基于 上 述 优点 ， 
以 及 它 的 容易 实施 ,所 以 在 其 提出 之 后 的 这 些 年 有 了 长 足 的 发 展 . 许多 学 者 使 用 TU 方法 
求解 各 种 不 同 的 方程 ,并 且 将 TU 与 其 他 迭代 算法 相 结合 ,以 获得 更 高 效 的 算法 . 下 面 着 
重 介 绍 两 种 增 量 未 知 元 : 二 阶 增 量 未 知 元 和 类 小 波 增 量 未 知 元 . 


1.1.1 一 阶 增 量 未 知 元 


iz o-[0.1] Ad NEN 为 细 网 格 步 长 ,2h 为 粗 网 格 步 长 , 微分 方程 为 Dirichlet 


边界 条 件 ,U 表示 结 点 未 知 元 ,其 元 素 为 usas 1a,B<2N 一 1, 一 阶 增 量 未 知 元 USE X. 


增 量 未 知 元 及 其 预 处 理 迭 代 算 法 


如 下 : 
定义 1.1 增 量 未 知 元 UG Yzi,2j + Z2i+1,2j + Z2i,2j+1 *Z2i+1,2j+1 ( 1<i,j<N—-1 ) 四 类 元 
素 , 其 中 ， 
V2i2; T Uzi,2j 
C2i2jtl — Uzi2jtl T U22j 
Zap. 一 Urea; — Uo 
Zilja 一 Uzij+ T U2s2j 


由 于 zi 为 wi 在 其 某 些 相 邻 点 的 增 量 , 所 以 称 U 为 增 量 未 知 元 ,在 定义 1. 1 下 


为 ui, 在 其 两 个 相 邻 点 处 的 增 量 ,利用 Taylor 公式 ,可 以 算出 在 细 网 格 点 (a B) Wb. ze 


OCh) ,因此 ,我 们 称 U 为 一 阶 增 量 未 知 元 . 
1.1.2 二 阶 增 量 未 知 元 
1. 一 维 情形 


BW O=[0.1] b= NEN 为 细 网 格 步 长 ,21 为 粗 网 格 步 长 , 见 图 1.1. 
o X 9 X O X o 


图 1.1 一 维 区 域 9 的 剖 分 : X 表 示 粗 网 格 点 ,* 表 示 细 网 格 点 


ij 


BU 


假设 微分 方程 为 Dirichlet 边界 条 件 ,U 表示 结 点 未 知 元 ,其 元 素 为 uas 1<a<2N 一 1, 一 


维 二 阶 增 量 未 知 元 U 定 义 如 下 ; 
定义 1.2 一 维 二 阶 增 量 未 知 元 也 包含 w ,zz (1 二 i 过 N 一 1) 两 类 元 素 , 其 中 ， 
ls = Usi 


1 
Za = Uzim — > (Uzin + uni) 


2 


利用 Taylor 公式 ,可 以 算出 在 细 网 格 点 a AÈ, z, =OCh?). A , 称 0U 为 二 阶 增 量 未 知 


元 . 下 面 ,以 一 维 Dirichlet 问题 的 数值 求解 问题 为 例 ,给 出 细 网 格 上 的 结 点 未 知 元 与 
未 知 元 之 间 的 转换 关系 . 
考虑 一 维 Dirichlet 问题 : 


—wW =f, 0«zr«l i 
en 一 u(1)—0 
利用 中 心 差分 格式 将 其 进行 离散 ,得 到 
( 


Uo = uzy 一 0 


ua — un =h’ fj, j=1,2,.…,2N—1 


He tit 


.2) 


Hf, =f Gh). 4 U,b€ RN? H uh fi Go1.2,72N—1) 为 按 字典 顺序 排列 成 的 


列 向 量 , 则 式 (1. 2) 等 价 于 线性 方程 组 


AU =b ( 


.3) 


其 中 ,A 为 (2N 一 1) X (2N 一 1) 三 对 角 和 矩阵 


2 1 
—1 2 =l 


=1 2. =i 


= 
W U RR U 中 元 素 按 偶 指 标 、 奇 指标 排序 后 的 近似 解 , 即 


o= (e) 


这 时 ,存在 排序 矩阵 P 了 满足 U=PU. 
这 里 ,P=(E， E,), 其 中 ， 


0 0 0 0 0 1 
1 0 0 0 0 0 
0 0 0 0 0 0 
B= 0 1 0 0 0 "T 0 
0 0 0 0 0 0 
0 0 0 0 1 0 
0 0 0 O0 © lZumc-oxas-o 0 
因而 ,有 
APU = b. P'APU = P'b 
^ 


从 而 得 到 另 一 线性 方程 组 


AU =b 
Mid 
P" 
U= 
Z; 
由 一 维 增 量 未 知 元 的 定义 可 知 
U — SU 


这 里 ， 


T 


o =O 0 


o cO o o0 o 


o ooo o 


ooo 8 0o 


(GN-DXN 


(1.4) 


(1.5) 
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dle |= 


w 
2 
i 
2 
A NXON—D EE. 于 是 ,线性 方程 组 (1. 5) 又 变 成 
ASU =b, S'ASU = S'b 
最 终 ,我 们 有 AU —b ,其 中 ， 
A =S'AS, b —S'b 
这 里 A EORR TE XE A [A «AT VJ RJ JA B BE HE GEORG. 
2. 二 维 情形 


设 0 一 [0,1 了 必 一 去 为 细 网 格 步 长 ,21 为 粗 网 格 步 长 , 见 图 1.2. 

设 微分 方程 为 Dirichlet 边界 条 件 ,U 表示 结 点 未 知 元 ,其 元 素 为 wp, 1<a,B<2N 一 1， 
二 维 二 阶 增 量 未 知 元 U 定义 如 下 : 

定义 1.3” 增 量 未 知 元 U 包含 yag tuens zagnan (OE EUN — D RIE 


素 , 其 中 ， 
V2i2j; — Uzi,2j 
Xa2pa 一 Uzija 一 z (Uzi,2jt2 Bx Uo; 2; ) 
= 1 [S 
Xia. 一 2a. 一 z Cusicz.zj ED Up; 2; ) 
ER 1 
Zia. 一 U2it1,27+1 7 FE Cuzisz,z T U2i,2j d U2i,2j+2 十 Urit2,2;+2) 


车 记 网 格 点 (2ih,(2j 十 1)h) 为 Fl 型 细 网 格 点 ,((2i 十 1)h,2jh) 为 F2 型 细 网 格 点 ， 
(2+ Dh. Cj HDA Ay F3 型 细 网 格 点 , 则 三 种 细 网 格 点 的 类 型 见 图 1. 3. 


o o o o o 
o x o x o 


o o o o o 


x 
o x o x o H " x x 
x o X o 
o o o o o . x x x 
* o3 s ` : a à Fl 型 F2 型 F3 型 
图 1.2 二 维 区 域 0 的 剖 分 : X 表 示 粗 网 格 点 ， 图 1.3 细 网 格 点 的 类 型 


“表示 细 网 格 点 


这 里 需要 指出 的 是 ,由 于 微分 方程 为 Dirichlet 边界 条 件 , 因 此 , 若 a=0 3X 2N .8—0 
或 2N, 则 wp 一 0. 
下 面 ,以 二 维 Dirichlet 问题 的 数值 求解 问题 为 例 ,给 出 细 网 格 上 的 所 有 结 点 未 知 元 
与 增 量 未 知 元 之 间 的 关系 . 
考虑 二 维 Dirichlet 问题 : 
— Au = f(x,y), (zy) E N = (0,1? 
f - C17) 
这 里 ,30 09 Q 的 边界 , 取 f(x,y) 二 2(x 十 y 一 x? — y*0. 利用 五 点 中 心 差分 格式 将 其 进行 
离散 ,得 到 


a Ui j-i — Up — Up — Ue; hh? fig stef = 11250, 2N—1 


(1,8) 


Uo, = uzw,j 一 0， uio = tin = Ovisj = 1,2,°°,2N—1 


HP, Sfi f Gh.jlo. BU BE RON 表示 wj sh? fi Gf 1.2, 2N— 1) 按 自然 顺 
序 ( 从 左 往 右 ,从 下 往 上 ) 排 列 成 列 向 量 , 则 式 (1. 8) 等 价 于 线性 方程 组 


AU =b (1.9) 
其 中 4 为 (2N 一 1)2X(2N 一 1): 和 矩阵 ， 
C —I 
—I C 一 T 
a= —I c —I | 
—-I C -I 
-I C 
这 里 ,IT 为 (2N 一 1) X (2N 一 1) 单位 矩阵 ,C 为 (2N 一 1)X(2N 一 1) 三 对 角 矩 阵 ， 
4 一 1 
—1 4 —1 
b= =j 4 =Ï | 
—1 4 一 ! 
—1 4 


id URR U 中 元 素 按 粗 网 格 点 .Fl 型 .F2 型 .F3 型 细 网 格 点 顺序 排序 后 的 近似 解 , 即 
a 
U; 


U = PU 
这 里 ,P= (EE E,@E, E, ®E, E;G)E;).E, . E; ARA. 4) ,因而 有 
APU = b. P'APU = P'b 
4 A—P'AP.b —P'b. Kil ARE BA 
AU—b (2.10) 


这 时 ,存在 排序 矩阵 己 满足 
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xieó-( ) wise 1.3 可 知 
Z; 
U = SŪ 
这 里 ， 
ra 0 
E Gy ® Ia Iyon 0 
o Iya G Gyn 0 Ius-oxsw-» 0 
G 69 Gap 0 0 Lea 


于 是 ,线性 方程 组 (1. 10) 又 变 成 
ASU =b, S'ASU = Sb 
最 终 , 有 
AU-b Q.1D 

Hp, A —S"AS.b =S"b. 
1.1.3 类 小 波 增 量 未 知 元 

类 小 波 增 量 未 知 元 (Wavelet-like Incremental Unknowns, WIU) fH Temam 首次 提 
出 . 与 上 述 一 阶 ` 二 阶 增 量 未 知 元 所 不 同 的 是 , WIU 可 以 将 近似 解 空间 分 成 两 个 L? 正 交 
子 空间 ,使 得 在 求解 微分 方程 的 过 程 中 能 够 自动 消 掉 一 些 项 ,大 大 简化 了 计算 . 理论 分 析 
表明 ,基于 WIU 的 数值 格式 的 稳定 性 条 件 , 要 明显 好 于 标准 算法 的 稳定 性 . 下 面 ,我 们 给 
出 在 一 维 ,二 维 情形 下 WIU 的 定义 . 


1. 一 维 情 形 


设 Q=[0,1j,h= NEN 为 细 网 格 步 长 ,24h 为 粗 网 格 步 长 . 一 维 类 小 波 增 量 
未 知 元 定义 如 下 . 


定义 1.4 类 小 波 增 量 未 知 元 UBA wi ,zzim(1 过 iN 一 1) 两 类 元 素 , 其 中 ， 


1 
2N+1 


yzi = i Cura + usi) 
(1.12) 


1 
Xe = z Cuz = uzi) 
2. 二 维 情形 


设 Q=[0,1J,h 一 FT*NEN 为 细 网 格 步 长 ,21 为 粗 网 格 步 长 ,二 维 类 小 波 增 量 
未 知 元 定义 如 下 . 


定义 1.5 二 维 类 小 波 增 量 未 知 元 UU 包含 Yziszj * T2127 s T2121» T2121 1 i C 
N) 四 类 元 素 , 其 中 


1 
V2i2; 一 E (Uzis2; 十 ziii 十 Uzija 十 Waeaszi? 


ED 
C2i1,2) 一 Ur — Y2i,2j 4 (uiia; 一 Uzi,2j 一 Uzi, 2j Uzi,2j1 ) 
1 (1. 13) 
2i,2j—1 Uzi,2j-1 — Y2i,2; 4 (3uzi,2j-1 Wi = Uiiig — Wa) 
1 
Z2i1,2j-1 U2i-1,2j-1 V2i2j 4 Curia oig — Uaia.ej — Usi) 
由 定义 1.5, 对 于 1 入 ;入 六 ,可 以 得 出 结 点 未 知 元 和 增 量 未 知 元 的 转换 关系 为 

Uzi.2j Y2i2j X2i—1,2j 22i,2j-1 Z2i-1,2j-1 

U2i-1,25 一 y2i21 | Zzi-1,2j 

M2i,2 产 1 一 y2i,21 | TZ2i,2j—1 

U2i-1,2j-1 = Yzizj s Z2i-1,2j-1 


以 上 三 种 增 量 未 知 元 方法 均 为 对 区 域 O 进行 二 层 离散 . 若 对 求解 区 域 进 行 多 层 离 
散 , 可 以 类 似 地 给 出 三 维 情形 下 二 阶 增 量 未 知 元 的 定义 和 转换 矩阵 结构 . 


1.2 数值 代数 方法 


科学 与 工程 计算 的 许多 重要 领域 ,如 流体 力学 计算 材料 模拟 与 设计 、 电 磁场 计算 、 电 
力 系 统 的 优化 设计 、 数 值 天 气 预报 及 核 爆 数值 模拟 等 都 离 不 开 微分 方程 的 求解 ,而 求解 微 
分 方程 可 以 通过 有 限 体积 有限 元 或 有 限 差分 法 等 对 微分 方程 进行 离散 ,化 为 求解 非 线 性 
方程 组 或 稀疏 线性 方程 组 . 非 线性 方程 组 的 求解 又 可 通过 Newton 迭代 或 简单 迭代 等 线 
性 化 ,化 为 稀 踢 线性 代数 方程 组 的 求解 . 因此 , 稀 朴 线性 代数 方程 组 的 高 效 求解 是 计算 数 
学 中 十 分 重要 的 课题 之 一 . 

目前 ,对 于 线性 代数 方程 组 的 求解 方法 一 般 分 为 两 类 : 直接 法 和 和 迭代 法 . 直接 法 的 
工作 主要 集中 在 20 世纪 六 七 十 年 代 , 直 接 法 需要 对 系数 矩阵 进行 分 解 ,在 不 计 舍 入 误差 
的 情况 下 经 过 有 限 次 运算 后 可 求 得 方程 组 的 精确 解 . 因此 ,直接 法 又 称 为 精确 法 . 对 于 
一 些 中 小 型 线性 代数 系统 , 常 选用 直接 法 求解 . 在 实际 应 用 中 ,直接 法 具有 高 精度 性 、 稳 
定性 和 可 操作 性 ,从 而 使 得 人 们 在 很 长 时 间 内 都 倾向 于 使 用 直接 法 求解 线性 代数 系统 . 
但 是 人 们 逐渐 发 现 , 当 线 性 代数 系统 中 系数 矩阵 的 条 件数 很 大 时 , 舍 人 误差 就 会 影响 求解 
的 准确 性 , 求 得 的 解 往往 与 精确 解 相 差 甚 远 . 另外 由 于 存储 的 问题 ,直接 法 需要 耗费 大 量 
的 计算 时 间 和 存储 空间 ,计算 效率 非常 不 理想 . 

EAC EIA 20 世纪 后 半 叶 开始 逐步 发 展 起 来 的 . ff Be UE [CI YE Gauss. SEE 
上 ,Gauss 意识 到 在 实际 计算 中 当 变 量 数 很 大 时 ,直接 法 是 很 难 完成 的 . 1823 年 ,他 在 给 
Gerling 的 一 封 信 中 很 兴奋 地 向 其 推荐 一 种 自修 正 算法 迭代 算法 ,他 描述 说 ,对 于 这 种 非 
直接 法 迭代 法 ,人 们 甚至 可 以 在 半 睡 半 醒 或 考虑 其 他 问题 时 完成 . 此 算法 后 来 又 被 Seidel 
于 1879 年 独立 发 现 ,这 就 是 后 来 闻名 于 世 的 Gauss-Seidel 迭代 法 . 

迄 代 法 是 采用 逐次 逼近 的 方法 ,从 一 个 初始 向 量 出 发 ,按照 一 定 的 计算 格式 ,迭代 计 
算出 逼近 未 知 向 量 的 无 穷 序 列 ,序列 极限 就 是 线性 代数 系统 的 精确 解 .迭代 法 通常 在 经 
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过 较 少 的 迭代 步 数 之 后 就 可 以 获得 非常 精确 的 近似 解 ,从 而 在 整个 问题 求解 过 程 中 大 大 
减少 了 和 迭代 求解 时 间 并 降低 了 和 迭代 求解 代价 ,尤其 是 迭代 法 具有 充分 利用 系数 矩阵 稀 玻 
性 的 潜力 ,只 需 知 道 系数 矩阵 与 向 量 乘 积 的 计算 法 则 就 能 求解 相应 的 线性 方程 组 ,这 是 直 
接 法 无 法 比拟 的 . 因此 , 随 着 计算 机 的 出 现 和 计算 机 技术 的 快速 发 展 ,迭代 法 与 直接 法 相 
比 具有 更 广泛 的 应 用 空间 和 优越 性 . 

目前 常见 的 迭代 法 有 定常 迭代 法 和 非 定 常 迭 代 法 . 对 于 定常 迭代 法 ,最 常见 的 是 单 
分 裂 迭 代 法 , 它 是 由 系数 矩阵 的 分 裂 

A=M—N 
导出 的 迭代 法 : 
x^? = MO Nx +M"b, k— 02s 

其 中 MON 称 为 迭代 和 矩阵 ; xO 称 为 迭代 初始 向 量 , 对 M 取 不 同 的 值 就 会 得 到 不 同 的 选 
代 算 法 ,例如 ,Jacobi,Gauss-Seidel,SOR,SSOR 等 经 典 迭 代 法 . 定常 迭代 法 形式 简单 . 易 
FREE. Varga 和 Young 在 相关 文献 中 对 此 进行 了 详细 的 论述 . 从 这 些 和 迭代 法 经 典 著作 
中 可 以 清楚 地 看 到 ,迭代 矩阵 的 选取 对 于 和 迭代 法 的 收敛 性 具有 决定 性 作用 ,只 有 当 所 选 迭 
代 和 矩阵 的 谱 半 径 小 于 1 时 ,才能 保证 选 代 法 收敛 ,而 且 和 迭代 矩阵 的 谱 半 径 越 小 ,和 迭代 收敛 
速度 就 越 快 . 因此 ,对 于 同一 个 线性 代数 系统 ,可 以 通过 比较 和 迭代 矩阵 谱 半径 的 大 小 来 判 
断 不 同 迭 代 方 法 收敛 的 速度 ,进而 为 迭代 方法 的 选择 提供 依据 ， 非 定常 迭代 法 的 典型 代 
表 是 子 空间 法 , 主要 包括 Conjugate Gradient (CG), Minimal Residual method 
(MINRES) „Generalized Minimal Residual method(RMRES) , BiCGStab 等 迭代 法 . JE 
常 迭 代 法 的 收敛 速度 依赖 于 系数 矩阵 的 谱 分 布 ,分 布 越 集中 ,收敛 速度 越 快 . 但 当 谱 分 布 
很 分 散 时 ,一 般 收敛 速度 很 慢 , 其 至 不 收敛 ,在 此 种 情况 下 ,可 以 对 线性 代数 系统 采用 预 处 
理 技 术 使 系数 矩阵 谱 聚 集 . 预 处 理 技术 是 解决 该 收敛 性 问题 的 有 效 途 径 , 成 为 迭代 法 研 
究 中 热点 问题 之 一 . 

目前 ,对 于 基于 和 矩 阵 分 柳 的 定常 迭代 法 虽然 直接 的 应 用 比较 少 ,但 是 结合 其 他 方法 的 
应 用 还 是 很 多 的 ,而 且 在 理论 上 有 重要 的 意义 . 这 里 ,主要 介绍 近 几 年 , 白 中 治 和 Golub 
等 提出 的 HSS 类 和 迭代 方法 . 


1. HSS 迭代 方法 


考虑 下 列 大 型 稀 玖 线性 方程 组 
Ax =b (1.14) 
Hp ACC” AK hibit FF ox DeC". 
求解 式 (1. 14) 的 和 迭代 方法 有 很 多 ,例如 经 典 的 Jacobi 迭代 和 Gauss-Seidel 迭代 方法 . 
FEF HERE A 的 对 称 / 反 对 称 分 裂 , 即 
A=H+S 


其 中 H— (RAT S= AA" ), 白 中 治 和 Golub 等 提出 了 对 称 /反对 称 分 列 


(Hermitian/skew-Hermitian Splitting, HSS) 方 法 . 
HSS RAR: 给 定 一 个 初 值 xO EC", 对 于 k= 二 0,1,2,…, 直 到 {x 咏 ) 收 敛 , 计 算 


A 
GI + Hox?) = (af 一 S)x +b 


(1.15) 
(al -- S)x €? = (al — IDx 3) +b 
其 中 a0. 
关于 HSS 迭代 方法 的 收敛 性 ,有 下 述 定理 . 
定理 1.1 4ACCU 是 一 个 正定 矩阵 ,H La FA' ,S Lu At) 分 别 是 它 


的 Hermite 部 分 和 反 Hermite 部 分 ,a 是 一 个 给 定 的 正 数 , 则 HSS 3E (CU i f 3x OB I 
M(a) 可 以 写成 


Ma) = (al +S$) (aI — H) GI + H) aI — S) 
它 的 谱 半 径 满足 下 面 的 不 等 式 
p(M(a)) <o(a) <1 


JEH ola) 一 max ek | canere H iit. 即 对 27-0. HSS 送 代 方法 收 仇 到 线性 


方程 组 (1. 14) 的 唯一 解 x* € C". 

定理 1. 1 说 明了 HSS 收敛 速度 以 ol(a) 为 上 界 ,而 ola) R5 RUER A 的 Hermite 
部 分 H 的 特征 值 有 关 ,而 与 4 的 反 Hermite 部 分 S 的 特征 值 以 及 和 矩阵 4 «H «S. 的 特征 向 
量 无 关 . 对 于 迭代 方法 中 参数 a 的 计算 方法 ,有 下 面 的 定理 : 


定理 1.2 令 AEC”" 是 一 个 正定 矩阵 ,HH L(A+A").S La At ) 分 别 是 它 的 


Hermite 部 分 和 反 Hermite 部 分 ,a 是 一 个 给 定 的 正 数 ,yx 和 xun 分 别 是 矩阵 H 的 最 大 
和 最 小 特征 值 , 则 


a” 三 arg min { max a—À | j = VY min mex 
@ Un» ET 
AL 
. Ymax — V min «(H) —1 
ola" ) 
Yeux 十 V uin (H) +1 
其 中 ,x(H) 是 和 矩阵 H 的 条 件数 . 
2. NSS 和 迭代 方法 


利用 矩阵 4 的 正规 /反对 称 分 裂 (Normal/skew-Hermitian Splitting. NSS) ,将 HSS 
方法 进一步 推广 , 即 
A=N+S 
其 中 ,NEC"" 为 正规 矩阵 ,SEC"" 为 反对 称 和 矩阵 ,得 到 了 NSS 方法 . 
NSS ERAGE: 给 定 一 个 初 值 x” EC" ,对 于 有 & 王 0,1,2,… EA GC } 收 敛 ,计算 
(al + NOx 3) = (al — Sx? +b 
(al + S)x**? = (al — N)x**? +b 
HP a >0,NEC™" H EREE SEC" IR Xt PE. 
KF NSS 迭代 方法 的 收敛 性 ,有 下 述 定理 . 


(1. 16) 
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定理 1.3 令 4EC”…" 是 一 个 正定 矩阵 ,NEC"” "是 一 个 正规 矩阵 ,SEC"” "是 一 个 反 
XY As ME IF ELL ASNES. 6; 为 N 特征 值 的 集合 . a 是 一 个 正常 数 . 此 时 NSS 方法 
和 迭代 和 矩阵 的 谱 半 径 满足 下 面 的 不 等 式 : 


oMa) Sola) 1. Vac 0 (1.17) 
其 中 ， 
|a—Al (a — Y;)* +45 
O PE Tata | sme GET Eq ices 


这 里 ,i /—1. 即 由 NSS 方法 得 到 的 迭代 序列 (x) 收敛 到 线性 方程 组 (1. 14) RO 
x. 除 此 之 外 , 令 yon 和 yoo Pin HI mo 分 别 表示 N 的 特征 值 实 部 绝对 值 的 最 小 值 和 最 大 
值 , 虚 部 绝对 值 的 最 小 值 和 最 大 值 ,2= [Yin aes ] X Erin] ,有 
二 

A Fein F pius 5 Wes ch ee e 


了 
2 


p= 24 2 
a' =argmin) max (S I +i |= 
a Op pea GT y; + Nj 


(1.19) 
相应 地 ， 
1 
min max min /max ^. m = 
uin + Yaax + 2 A uin Yous. — Mg 
aa") = a. 20) 


VV hin T- Tux — e] . > V I TN 
( VY mint Nmax 十 Yin dn NS 
定理 1.3 说 明了 ola) 只 与 正规 矩阵 N 的 特征 值 有 关 ,而 与 4 的 反 Hermite 部 分 S 
的 特征 值 以 及 矩阵 4,N,S 的 特征 向 量 无 关 . 


3. AHSS 和 迭代 方法 


2004 年 , 白 中 治 等 提出 一 种 预 处 理 HSS (PHSS) 人 迭代 法 用 于 数值 求解 Hermite 半 正 
定 线性 鞍点 系统 . 2007 年 ,为 了 提高 PHSS 迭代 法 的 收敛 速度 ,他 们 在 HSS 方法 中 引入 
双 参 数 , 用 来 求解 鞍点 问题 , 称 为 加 速 HSS(AHSS) 迭 代 法 . 

AHSS 迭代 方法 : 给 定 一 个 初 值 XO EC", 对 于 ==0,1,2,…, 直 到 {x | IE 


"s 
(al -- IDx *). = (ql —S)x® +5 


(BI +$)x* = (gr ix C3 45 
其 中 a,8 为 给 定 正 数 ,I 为 与 矩阵 A 同 阶 的 单位 矩阵 . 

SB ACC" AE Hermite 正定 矩阵 时 ,可 以 证 明 , 若 a,B 满 足 一 定 的 条 件 ,AHSS 
方法 的 迭代 和 矩阵 的 谱 半径 小 于 1. 杨 爱 利和 伍 渝 江 等 对 该 方法 的 收敛 性 进行 了 更 加 细致 
的 讨论 ,并 详细 给 出 了 两 个 参数 的 计算 方法 . 数值 实验 证 明 ,AHSS 方法 较 HSS 方法 更 
加 有 效 和 稳健 . 


4. LHSS 和 迭代 方法 
Li 等 在 2007 年 基于 HSS 格式 提出 一 种 偏向 一 侧 的 HSS (LHSS) 和 迭代 法 ,其 定义 


如 下 . 
LHSS 和 迭代 方法 : 5E — 1 WI fl x" EC”, 对 于 二 0,1,2,…, 直 到 {x 中 }) 收 敛 , 计 算 
Hx C? =— Sx? +b 
E (1.21) 
(al -- S)x*? = (aI — IDxÜ ?) +b 


其 中 a 为 给 定 的 正常 数 . Li 给 出 了 LHSS 方法 的 收敛 定理 和 参数 o 的 计算 方法 ,并 证 明 
在 一 定 条 件 下 ,LHSS 方法 的 谱 半 径 的 上 界 小 于 HSS 方法 谱 半径 的 上 界 . 

HSS 类 方法 也 可 以 用 来 求解 非 线性 代数 方程 组 . 2010 年 , 白 中 治 等 利用 HSS 迭代 
方法 求解 Newton 方程 组 ,提出 了 Newton-HSS 方法 .大量 数值 实验 证 明了 Newton-HSS 
方法 在 迭代 步 数 和 CPU 时 间 上 都 要 优 于 Newton-USOR, Newton-GMRES 以 及 
Newton-GCG 方法 . 2011 年 , 郭 学 萍 等 分 析 了 不 精确 Newton-HSS 方法 的 半 局 部 收敛 
性 ,并 通过 引入 向 后 线 搜索 技术 ,给 出 了 全 局 收敛 的 不 精确 Newton-HSS 方法 . 除 此 之 
外 ,对 于 弱 非 线性 代数 方程 组 ,为 了 进一步 提高 Newton-HSS 方法 的 求解 速度 , 白 中 治 等 
利用 Picard 迭代 作为 外 迭代 ,HSS 方法 作为 内 迭代 ,提出 了 求解 弱 非 线性 代数 方程 组 的 
Picard-HSS 和 非 线性 HSS-like 迭代 方法 ,这 两 类 方法 在 计算 的 过 程 中 均 不 需要 在 当前 
迭代 步 x Mbit AU WK Jacobi XB RE, AL ZEA GEAR HSS 方法 中 出 现 的 两 个 子 线性 代数 方 
程 组 的 系数 矩阵 不 随和 迭代 步 数 的 变化 而 变化 . 因此 ,与 Newton-HSS 方法 相 比 ,这 两 类 方 
法 大 大 减少 了 算法 的 工作 量 ,提高 了 问题 的 求解 速度 . 这 期 间 ,很 多 博士 研究 生 的 毕业 论 
文 也 对 非 线性 方程 组 的 解法 进行 了 大 量 的 研究 ,对 数值 代数 理论 的 发 展 起 到 了 一 定 的 促 
进 作用 . 


1.3. 本 书 的 主要 工作 


本 书 主要 包含 以 下 两 方面 的 工作 . 

第 一 部 分 是 对 二 阶 增 量 未 知 元 和 类 小 波 增 量 未 知 元 预 处 理 方法 的 研究 . 

对 于 非 线性 特征 值 问题 ,结合 二 阶 增 量 未 知 元 方法 和 Marder-Weitzner 算法 ,提出 了 
4 种 修正 的 Marder-Weitzner 算法 ,给 出 了 算法 的 收敛 性 分 析 , 数 值 实验 证 明了 新 算法 的 


正确 性 和 有 效 性 . 对 于 各 向 异性 的 反应 扩散 方程 ,因为 导数 了 项 含有 非常 小 的 参数 , 通 


常 使 用 类 小 波 增 量 未 知 元 (WIU) 求 解 这 类 方程 不 是 特别 有 效 . 因此 ,本 书 提出 了 一 种 基 
于 矩阵 分 块 的 新 的 增 量 未 知 元 : 分 块 类 小 波 增 量 未 知 元 (WBIU)，WBIU 一 方面 保持 
WIU 的 优点 ,通过 将 近似 解 空 间 分 解 成 两 个 L^ 正 交 子 空间 来 自动 消 掉 有 限 差分 离散 等 
价 方程 组 的 一 些 项 ,以 简化 计算 . 另 一 方面 ,新 的 分 块 类 小 波 增 量 未 知 元 仅 在 一 个 方向 
(系数 较 大 的 方向 ) 引 入 WIU ,因此 与 一 维 的 WIU 类 似 ,这 样 可 以 有 效 地 减少 系数 矩阵 的 
条 件数 ,节省 求解 的 时 间 . 除 此 之 外 ,本 书 还 建立 了 新 的 范 数 估计 ,给 出 两 种 WBIU 型 数 
值 格式 : 显 格式 和 半 隐 格式 ,并 证 明了 在 新 的 范 数 估计 下 ,格式 是 稳定 的 . 最 后 ,用 数值 
算 例 对 理论 分 析 进 行 了 验证 . 
第 二 部 分 研究 了 几 种 增 量 未 知 元 预 处 理 迭 代 算法 . 


labs 


gl 
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Ts. dE TK LE Hermite 正定 线性 方程 组 的 预 处 理 NSS 迭代 方法 ,然后 在 
该 方法 的 基础 上 ,引入 双 参 数 , 提 出 了 两 参数 预 处 理 NSS 迭代 方法 . 理论 分 析 表 明 ,在 一 
定 条 件 下 ,新 的 迭代 格式 是 收敛 的 . 除 此 之 外 ,还 精确 地 给 出 了 迁 代 格式 出 现 的 两 个 参数 
和 迁 代 矩阵 谱 半 径 的 最 小 上 界 的 计算 方法 . 对 于 复 对 称 线性 代数 方程 组 ,为 了 提高 其 求 
解 效率 ,基于 NSS 超 松弛 加 速 欠 代 格式 ,提出 一 类 求解 复线 性 代数 方程 组 的 超 松弛 迁 代 
MHSS 加 速 迭代 格式 ,然后 推导 了 该 加 速 迭 代 方 法 的 收敛 条 件 以 及 加 速 迭 代 格 式 中 参数 
co 的 选取 办 法 ,最 后 用 数值 实验 证 明 方 法 的 正确 性 和 有 效 性 . 对 于 大 型 稀 玻 非 线 性 代数 
方程 组 ,提出 求解 非 线 性 代数 方程 组 的 Newton-LHSS 方法 、Picard-AHSS 方法 、Picard- 
MHSS 方法 和 非 线性 AHSS-like 方法 ,并 给 出 了 4 种 方法 的 收敛 性 定理 ,数值 实验 证 明 
了 理论 的 正确 性 和 方法 的 有 效 性 . 


^2 
基于 增 量 未 知 元 方法 的 非 线性 
特征 值 问题 的 研究 


2.1 非 线性 特征 值 问题 简介 


非 线性 特征 值 问题 在 科学 计算 的 很 多 领域 中 都 有 广泛 的 应 用 ,如 流 固 耦 合 振动 数值 
模拟 问题 ,加 速 器 设计 中 的 电磁 场 数值 模拟 等 . 很 多 学 者 对 非 线 性 特征 值 问题 进行 了 广 
泛 而 深入 的 研究 ,并 且 取 得 了 丰富 的 成 果 . 目前 ,求解 非 线 性 特征 值 的 方法 有 自 适应 多 尺 
度 方法 ,多 重 网 格 算法 迭代 投影 法 等 . 

本 章 主要 考虑 如 下 形式 的 非 线 性 特征 值 问 题 : 


— Au = YF (u), QE (0,1)".n — 1.2 
| (2.1) 


u| — 0 


其 中 y>0.F; R>R 是 C? 映 射 , 满 足下 列 条 件 : 
Qro-o Qr(»-1 © lim (EGO: Soo 


tuto dull 

Picard 和 迭代 方法 是 求解 不 动 点 方程 的 有 效 方法 之 一 . 对 于 方程 (2. 0 ,虽然 单纯 利用 
Picard 方法 并 不 能 正确 得 到 其 解 ,但 车 在 连续 两 次 进行 Picard 迭代 后 ,引入 松弛 因子 a, 
通过 适当 选取 a 的 值 , 便 可 以 计算 出 方程 (2. 1) 的 解 ,这 就 是 所 谓 的 Marder-Weitzner 77 
ik. 众所周知 ,IU 方法 是 进行 有 限 差 分 时 通 近 惯性 流 形 的 一 种 近似 方法 , 当 对 求解 区 域 
进行 多 水 平 离散 时 ,用 粗 网 格 点 的 增 量 代替 相应 细 网 格 点 的 未 知 函 数 的 值 ,从 而 , 增 量 未 
知 元 函数 在 不 同类 型 的 点 产生 不 同 的 矩阵 结构 。 这 种 特殊 的 矩阵 结构 对 原 代 数 系统 起 到 
了 预 处 理 作 用 ,从 而 缩短 了 计算 时 间 . 我 们 的 目标 是 利用 增 量 未 知 元 的 多 尺度 性 质 ,对 
Marder-Weitzner 方法 进行 修正 . 

本 章 在 已 有 成 果 的 基础 上 ,对 非 线性 特征 值 问题 进 行 研究 . 首先 回顾 了 与 本 章 相关 
的 线性 Richardson 方法 和 Marder-Weitzner 方法 ,然后 通过 引入 一 种 新 的 A 范 数 , 给 出 松 
弛 因子 a 的 4 种 计算 方法 , 称 为 修正 的 Marder-Weitzner 方法 ,并 给 出 收敛 性 定理 . 最 后 
在 数值 算 例 中 ,利用 增 量 未 知 元 方法 ,对 网 格 进行 多 水 平 离散 ,数值 结果 证 明了 新 方法 的 
正确 性 和 有 效 性 . 
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2.2 线性 Richardson 方法 和 Marder-Weitzner 方法 


利用 有 限 差分 方法 将 方程 (2. 1) 进 行 离散 ,得 到 下 列 非 线 性 代数 方程 组 : 
AU = yF U) (2.2) 
这 里 ,4 为 一 A 算 子 离散 后 得 到 的 矩阵 ,下 : R"—R" JE C? 映射 ,y 0. 除 此 之 外 ,为 了 计 
算 方便 ,固定 参数 7, 则 可 以 得 到 与 方程 组 (2.2) 等 价 的 不 动 点 方程 : 
U = TQ) (2. 3) 
式 (2.3) 中 ,T=74 TF.» ). 


2.2.1 线性 Richardson 方法 


考虑 下 列 线性 方程 组 
Mx =b, xe R (2.4) 
HpMER""ADIEERSIExEAXBP. 
线性 Richardson 方法 : 给 定 一 个 初 值 xO ER" ,对 于 80.1.2, ERI (x9 SH, 
计算 x ** P 一 x +a(b—Mx™ ) ,其 中 非 零 实 参数 a 称 为 迭代 参数 . 
Richardson 迭代 格式 收敛 的 充分 条 件 是 


2 
0 IM. 
+o D 


其 中 ,oCM) 为 矩阵 M 的 谱 半径 . 
2.2.2 Marder-Weitzner 方法 


考虑 非 线 性 特征 值 问 题 方程 组 (2. 2) ,我 们 知道 ,仅仅 利用 Picard 迭代 
ue = T(U®) 
不 能 得 到 方程 组 (2. 2) 的 不 稳定 解 . 事实 上 , 若 设 U" AIEA. 20 A a e =U” 
—U* 为 第 m 次 Picard 迭代 的 误差 , 则 利用 Taylor 公式 ,可 以 得 出 
gn 一 We” + o(e/? ) 

Hp WARE THEU 处 的 Jacobi 矩阵 . 为 了 保证 Picard 3E (Cr iMi Sx. MEA oC) — 1. 
但 在 实际 问题 中 , 当 U * 为 不 稳定 解 时 ,这 个 条 件 很 难 满足 . 因此 ,Marder 和 Weitzner 提 
出 了 下 面 的 三 步 迭 代 格 式 . 

Marder- Weitzner 方法 : 给 定 一 个 初 值 U ER ,对 于 k= 二 0,1,2,…, 直 到 (UU) 收敛 ,计算 


y Co) = TU”) 
yt) =- (vt?) — 
U^ = ye emu? —poe -vyen) 


此 时 ,可 算出 Marder-Weitzner 方法 的 误差 方程 
e" — (E —aQU— 9!) &'-F-o(g') 
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为 了 保证 Marder-Weitzner 方法 收敛 ,应 有 
pO —aO— v)» )-1 (2. 6) 
因此 ,只 要 参数 a 足够 小 ,p( 亚 ) 可 以 大 于 1. 
对 于 上 述 Marder-Weitzner 方法 ,有 下 面 的 收敛 性 定理 . 
定理 2.1 BRAT TEMU PS. iW AT FEU" bY Jacobi 4, IERIE 
亚 的 特征 值 为 不 等 于 1 的 实数 . 那么 , 当 0 二 a<a. 时 ,存在 U” 的 某 个 邻 域 V., 当 U® € V, 
时 ,由 式 (2. 5) 定 义 的 序列 (UP SSB UB 


limu =u", i20,1,2 


k- oo 


除 此 之 外 ,如 果 记 ， 
a sup |1—y|, b inf |1— y] 
yE € spo 
则 有 
Ce 一 一 Q, 2 
a* " +h 


这 里 ,ao 代表 参数 a 的 最 优选 取 . 
2.3 Marder-Weitzner 方法 的 改进 


令 AER "为 对称 正定 矩阵 ,向 量 a 履 关于 A 的 向 量 积 定义 为 
(a.b), — a Ab. VYa,b € R* 
相应 地 , || © lA 为 A 内 积 的 从 属 范 数 . 


2.3.1 修正 的 Richardson 方法 


修正 的 Richardson 方法 : 给 定 一 个 初 值 xO € C" ET. 0.1.2.7 HERI GC } 收 
S SE 
xD = yO Eg (b — Mx?) 
= (Mr™ re, 
(Mr ,Mr® )4 
这 里 ,选取 松弛 因子 w 使 得 第 &A 次 迭代 的 残 量 函数 ea ll r® —a,Mr™ || ‘Buh. 
a, 的 值 随 着 迭代 步 数 & 的 不 同 而 不 同 . 


2.3.2 修正 的 Marder-Weitzner 方法 

类 似 地 ,在 Marder-Weitzner 方法 中 ,用 ww 代替 a ,我 们 提出 修正 的 Marder-Weitzner 
Jrik. 

修正 的 Marder-Weitzner 方法 : 给 定 一 个 初 值 0 ER" ,对 于 &A 一 0,1,2,…, 直 到 
(U? jugo ,计算 


Qk 
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y = qe) 


(HE) (#4) 
ee (v } (2.7) 


U^» 一 yp amt? 一 U -uD ) 


下 面 ,给 出 a 的 计算 方法 . 
WE k PERRE H rO =U" — TQU? ) ,利用 式 (2.7) ,可 以 得 到 re RI pU? fx dE 
KA: 


r*P —(I— a, 01 — V )*)r'? d o(r'? ) 
=(I—a,N)r™ + 0(r™) 


Az (I — a, NOr'? 


f 
其 中 N Sa sw y, didi || Oa NO | 取得 最 小 ,可 得 出 a 的 4 种 计算 方法 . 
Q 当 A 一 N 时 ,有 
(ON? P D 
ay = inte INEO] (2.8) 


© «A =N Rt, A 


a (Nr P or? Ya (ri? NU ) TP 
t0 (ONrP NEP) ONE? ur ) U 
© 当 A 一 [时 ,有 
(Nr? or”) 
a = (NEP NF?) £210) 
QD 当 A SNAN 时, 有 
(ON? p? pO) 十 Go wr) 
= (NE EP Nr) 4- (Nr r9) C2 m) 


ELC +. +) 表示 Euclid 范 数 , 也 就 是 A 范 数 中 A =I( 单 位 和 矩阵) 的 情形 . 
为 了 计算 w 的 值 , 需 要 计算 更 了 ,i 二 1,2,3,4 的 值 . 事实 上 ,补充 两 步 Picard 迭代 


u=) = T(v c) 
uH) = r(u c) 
可 以 得 到 
UHS) — U* PeH o (e?) 
US) =U" 4+ He + fe) 
J- r(vC*) =U" + e+ of) 
由 上 式 可 以 得 到 Wir ,i 二 1,2,3,4 的 近似 计算 方法 ,也 就 是 
Vir = UH —U* S tort), i=1,2,3 
Eo =U" —J+olr) 
利用 上 述 结果 可 以 得 到 Mr” 和 M2r® 的 近似 计算 方法 . Wi Mr? .W; =M r” , 则 有 
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Wi —ü—2 4 w*)r? 
= —2 Fr! re 


yu? —u (+) (ue ve Jue u) Fofr”) 


~Ut +3(Utt — uH) —ut 
W, r? 4 vr’ + 6 Pr” 4 rr”? + Pir” 
Uw —u(#s) (vue vc) | sue vc) 


(ve ue) EUCH) —J +o) 


AU +5(20 $) =u) ) —s(gu e ve) 一 了 


2.3.3 Al1 方法 

下 面 ,我 们 给 出 4 个 修正 的 Marder-Weitzner 方法 ,统称 为 Al 方法 . 

假设 U" 是 方程 (2. 1) 的 离散 形式 AU=F(U) 或 者 U=T(U) 的 解 , 初 值 UER" 是 解 
U' 的 一 个 近似 值 ,对 于 8—0.1.2.77 

步骤 1: UG) 二 T(U 695) ) ,i=1,2,3,4 B J=T(UC+#) ) 

步骤 2. r® =y® —y t+ 


WwW, =U” eve -ven) —u (Hs 
W: =U” ex(me -ut ) —s{2u 4) -veb)-; 


步 又 3: 计算 松弛 参数 aP 1 1.2.3.4. 
w (W, .r? ) 


a = (WW) 
a? " (ro wr? ) 
: (W, rw ) 
PE (2.12) 
o? — 315 
(Wi Wi) 
o — OV rm) GO pr) 
Qk o» 
(V; Wi) + (Wisr® ) 


步骤 4: Ut = Faf? (20 63) —u® —U (*$) ) 51152,8354, 


k=k+1 转向 步骤 1. 
下 面 给 出 方法 的 收敛 性 定理 . 
定理 2.2 BEAT TEM U 处 下 可 导 . WP AT EU 处 的 Jacobi 和 矩阵 ,并 假设 
思 的 特征 值 为 不 等 于 1 的 实数 . 那么 , 当 0 二 a 二 a. 时 ,存在 U” HED BM VL. 4 UO E 
V, 时 ,4 种 修正 方法 产生 的 序列 {DU } 收 敛 到 UU* , 即 
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limu (=+) zt. 4051.2 


除 此 之 外 ,如 果 记 
a sup |l—y|. b inf |1—yl 
yes vE DD 
则 有 
ae 一 2 Qop = 2 
oe a 二 


这 里 ,oow 代表 参数 a 的 最 优选 取 . 
2.4 数值 实验 
本 节 结 合 修正 的 Marder-Weitzner 方法, 分别 利用 CSIU( 六 阶 紧 格 式 增 量 未 知 元 )、 
LRA,MLR,MW.,AIIU 和 MWIUa 等 算法 计算 非 线性 特征 值 问题 . 
2.4.1 Dirichlet 问题 
考虑 下 面 经 典 的 Dirichlet 问题 


—^u-f, Q= [0,1F 
| E f ls [o (2.13) 
u=0, a 
为 了 方便 , 令 h=h, =h, WABCO. 13) 离 散 , 得 到 
Ung — 2Uij — auno us; — Zu — Ui j- " 
f (2.14) 
hi hi d 


APU — sin(162y(1 —a) (1 — 3) ,利用 自 适应 Richardson 方法 计算 . 3€ 2.1 给 出 了 
式 (2. 14) 系 数 和 矩阵 最 大 特征 值 和 最 小 特征 值 的 比值 ,而 当 我 们 用 TUS. 方法 的 时 候 , 和 矩阵 
条 件数 明显 降低 . 


表 2.1 系数 矩阵 条 件数 的 比较 


Initial grids Without CSIU With CSIU Initial grids Without CSIU With CSIU 
N=7 9. 4835 2.9951 N=25 35. 8304 2.3188 
N=17 24.5760 2. 3810 N=30 42. 0890 2. 2811 


2.4.2 非 线 性 特征 值 问题 


考虑 下 列 问题 ,找到 u E H* CO) T) HS G0. 满足 下 列 边界 值 问题 : 
ue SENS n= 0.1 (2.15) 
u=0, IN 

HeoBoy.c0.0—ex2. 我 们 的 目的 是 计算 方程 (2. 15) 的 不 稳定 解 . 

当 参 数 y 穿越 一 人 的 特征 值 时 ,方程 的 解 会 产生 分 歧 , 给 求解 带 来 困难 . 下 面 列 出 一 
些 解 的 相关 性 质 . 

S Ass v (p? Iq. pA 是 一 A 的 一 个 特征 值 , Doa = sinl prr)sinlgrr) 是 相应 
的 特征 函数 . 
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当 7y<Aii 时 ,方程 (2.15) 只 有 一 个 稳定 的 平凡 解 二 0. 
当 Ain 达 7 二 Ai,s, 平 凡 解 u 二 0 变 得 不 稳定 ,方程 (2. 15) 存 在 两 个 稳定 解 , 记 为 


K(1,1) ,它们 是 特征 函数 $1,1 的 变形 . 


2 Ap CY p rq 024& f. KOSD 型 解 ,所 有 解 (包括 平凡 解 ) 都 是 不 稳定 的 . 


S a TIL b HB a? +b? <p? aq? ,这 些 不 稳定 解 是 下 列 两 类 函数 的 变形 . 


ep 


一 A 算 子 的 特征 值 函数 D, ,相应 的 不 稳定 解 是 K Ca 0078 08. 为 了 计算 它们 , 取 


Uo 一 k@,s 作 为 初 值 . 


(2) 


Au 型 解 . 为 了 计算 这 种 类 型 的 解 , 取 初 值 
U,—60,,. O, = sin(anz)sin(bry)Z(x,y) 


其 中 ZCz,y) 在 两 条 平行 线 y= 二 x 和 yy 二 一 + 上 值 为 零 . 
在 Cs 网 格 O 上 给 出 数值 结果 , 令 参 数 y=120,v==120,e==1. 取 初 始 值 Ô’ = 


sin) 


sin(x|x— y|). $ 200 次 迭代 . 图 2. 1 和 图 2.2 给 出 200 次 迭代 的 数值 结果 ， 


图 2. 3 和 图 2.4 给 出 50 次 迭代 的 数值 结果 . 


x10" 10-9 
5.1711 3.3533 — : ; 
54711 33533 
5.1711 3.3533 
5.711 3.3533 
5.1711 3.3533 " 
5.711 3.3533 
5.1711 3.3533 
5.711 3.3533 
5.1711 3.3533 
5.1711 im aa a aien 3.3533 + +4111. a 
0 20 40 60 80 100 120 140 160 180 200 0 20 40 60 80 100 120 140 160 180 200 
EC stat 
X107? x10 
5.1847 T T T T 5.711 
5.1847 5.711 1 
5.1847 5.711 
5.1847 5.1711 上 
5.1847 5.1711 f 
5.1847 5.1711 上 
5.1847 5.1711 f 
5.1847 5.1711 上 
5.1847 5.1711 上 


5.1847 
0 


" " " " 5.1711 4 4 " " " " 1 " " 

20 40 60 80 100 120 140 160 180 200 0 20 40 60 80 100 120 140 160 180 200 
srap Bua 
图 2.1 四 种 情况 下 ANU 方法 迭代 步 数 VS 和 迭代 误差 的 比较 


C ”所谓 Cs.: 网 格 , 指 的 是 区 域 的 每 个 方向 上 初始 粗 网 格 点 数 为 1, 加 细 5 次 ,最 后 ,每 个 方向 细 网 格 点 数 为 25X 


GTD-1 


个 , 即 网 格 点 数 为 61. 
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在 数值 计算 的 过 程 中 ,对 网 格 进行 两 层 训 分 . 数值 结果 见 表 2.2 — 3x 2. 17. 这 里 注 
意 , 表 2. 9 中 数值 残 量 由 单 次 迭代 获得 . 注意 到 4 种 情况 (casel ,case2,case3,case4) 相 应 
于 A1 方法 中 4 种 松弛 参数 af? (i 二 1,2,3,4) 的 选择 . 这 里 ,N 表示 初始 网 格 点 ,TT 代表 
CPU 时 间 ,单位 是 秘 . 

表 2.2 LRA 方 法 的 迭代 残 量 


N T Res N T Res 
7 0. 0000 0. 5220 8 0. 0000 0. 4462 
9 0. 000 00 0. 4615 10 0. 0000 0. 4311 
11 0. 0000 0. 3842 12 0. 0000 0. 3746 
13 0. 0000 0. 3575 14 0. 0000 0. 3556 
15 0. 0000 0. 3378 16 0.0160 0. 3180 
17 0. 0150 0. 3208 18 0. 0000 0. 3177 
19 0. 0160 0. 3222 20 0.0160 0. 3016 
21 0. 0150 0. 2919 22 0. 0150 0. 2840 
23 0. 0160 0. 2870 24 0.0310 0. 2820 
25 0. 0320 0.2677 26 0.0470 0. 2675 
27 0.0780 0. 2532 28 0.0470 0. 2530 
29 0.0780 0.2479 30 0. 0940 0. 2498 
31 0. 2030 0. 2430 
表 2.3 MLR 方法 的 迭代 残 量 

N T Res N T Res 
rd 0. 0000 0. 2350 8 0.0160 0. 2057 
9 0.0320 0. 1900 10 0.0160 0. 1811 
11 0. 0000 0.1740 12 0.0150 0. 1690 
13 0. 001 50 0.1655 14 0.0160 0. 1690 
15 0.0310 0. 1604 16 0. 0310 0. 1604 
17 0.0310 0.1614 18 0.0470 0. 1667 
19 0. 0630 0.1631 20 0.0780 0. 1694 
21 0.0940 0.1662 22 0.1090 0. 1724 
23 0.1720 0.1702 24 0. 5620 0. 1764 
25 0. 7650 0.1764 26 0. 7500 0. 1807 
27 0. 9540 0.1793 28 0. 9690 0. 1854 
29 1.8900 0. 1841 30 0. 9530 0. 1902 
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N T Res N F Res 
7 0. 0000 0. 1157 8 0.0160 1. 1129 
9 0. 000 00 1.1111 10 0. 0160 1. 1097 
1i 0.0160 0. 1090 12 0.0160 1. 1083 
13 0. 0150 1.1078 14 0. 0150 1. 1074 
15 0. 0160 1. 1071 16 0. 0150 1. 1069 
17 0. 0310 1. 1067 18 0. 0160 1. 1066 
19 0. 0160 1.1064 20 0. 0470 1. 1063 
21 0. 0460 1.1062 22 0. 0470 1. 1061 
23 0. 0780 1.1061 24 0.0780 1. 1060 
25 0. 1100 1. 1060 26 0. 0930 1.1059 
27 0. 3910 1. 1059 28 0. 2500 1. 1059 
29 0. 4060 1.1058 30 0, 4530 1.1058 
31 0. 4220 1.1058 
表 2.5 SAMa hi Al 方法 的 迭代 残 量 
N(e=1) T Res N(e=1, 2) T Res 
7 0.3120 | 1.1475 X 10"! 7 0.0780 | 3.349910 
9 0.3750 | 1.105610" 9 0. 3280 1. 2464 X 10!* 
17 14.7500 | 2.9365:X10'5 17 14.7340 | 9.932610" 
25 151.8120 | 9.5161X10" 25 151.8750 | 6.8619 10% 
表 2.6 FHM? BF Al 方法 的 迭代 残 量 
N(e=1) T Res N(e—1.2) T Res 
f 0. 0940 1. 1627X 10" 7 0. 0930 303880 X 10" 
9 0.3120 | 1.1199X 10" 9 0. 3120 1. 6105 X 10" 
17 14.7340 | 3.2768X10'5 17 14. 6720 1. 0057 X 10° 
25 151.7660 | 7.6001:X10" 25 151.2810 | 6.9470X10^* 
表 2.7 SAMa h Al 方法 的 迭代 残 量 
N(e=1) T Res N(e—1.2) T Res 
1 0.0780 | 1.1483X10" z 0.0780 | 3.3513X10" 
9 0.3130 | 1.106410" 9 0. 3280 1. 6093 X 10" 
17 14. 6720 3. 2455 X 105 17 14. 6880 9. 9388 X 10?! 
25 151. 3280 9.5163 X 10” 25 152. 1410 6. 8661 X 105 
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R28 SAMa” 时 Al 方法 的 迭代 残 量 


N(e=1) T Res N(e—1.2) T Res 
7 0. 0940 1.1475 X 10"! 7 0. 0780 3. 3499 X 1055 
9 0.3120 | 1.105610" 9 0. 3130 1. 6092 X 10" 
17 14. 6560 3. 2437 X 10" 17 14. 6870 9. 9326 X 10?! 
25 151. 2500 9.516110!" 25 151. 9530 6. 8619 X 10% 
42.9 在 N=7 时 MWIU 方 法 的 迭代 残 量 (1 一 3 步 迭代 ) 

Scheme MWIU literation 2iterations 3iterations 
e=1 0. 9734 5. 3498 X 10"? 2.5112X 10" 
0. 9722 4. 5519 X 10" 4. 8392 X 10” 


52.10 SAMa h AUU 方法 的 迭代 残 量 


N(e=1) T Res N(e=1, 2) T Res 
7 0. 1090 1. 1086 X 10" 7 0. 1090 1. 6571X 108 
17 15.3120 | 3.6844X 10" 17 15.2650 | 4.7046X107 
25 100.1400 | 1.2251X10!* 25 157.4220 | 3.0070X102 
表 2.11 SMa? 时 ATIU 方法 的 迭代 残 量 
N(e=1) T Res N(e=1, 2) T Res 
7 0. 0940 2. 0743 X 10" 7 0. 0940 1. 1233 X 10!! 
17 14.9380 | 3.816710" 17 14.9530 | 4.99216 x10" 
25 89.7350 | 1.2356X 10" 25 156.4060 | 3.068510? 
R212 SAMa” h ATIU 方法 的 迭代 残 量 
N(e=1) T Res N(e=1, 2) T Res 
7 0. 0460 1. 7321X 105 7 0. 0940 1. 1094 X 10" 
17 10.0470 | 3.2437X10' 17 11.5000 | 4.7699xX10' 
25 91.1870 | 1.3451X10" 25 84.7810 | 3.0117X10^ 
R213 SMa” 时 ATIU 方法 的 迭代 残 量 
N(e=1) T Res N(e=1. 2) T Res 
7 0.0470 | 1.108610" 7 0. 1250 1. 6571 X10" 
17 15. 0930 3. 6844 X 10 17 14. 9220 4. 7046 X 10 
25 159.0160 | 1.2251X10!5 25 166.1410 | 3.007010” 
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R214 参数 取 al? 时 MWIUa 方 法 的 迭代 残 量 


N(e=1) T Res N(e—1.2) T Res 
7 0. 1410 3. 6548 X 10? 7 0. 0940 6.0978 X 10!* 
17 16. 6720 1. 2513X 1055 17 15. 7190 1.5269 X 107! 
25 178. 7810 7. 2213X 10" 25 177. 4380 1. 1438 X 105 


R215 参数 取 a 已 时 MWIUa 方法 的 迭代 残 量 


N(e=1) T Res N(e—1.2) T. Res 
7 0. 1560 3. 7204 10" 7 0. 1100 6. 2541 X 10" 
17 15. 7030 1. 2715X 10^ 17 15. 7030 1, 5482 X 107! 
25 197. 0470 7. 2258 X 10” 25 195. 8280 1.1587 X 10” 


R216 SMa H MWIUa 方法 的 迭代 残 量 


N(e=1) T Res N(e—1.2) T Res 
7 0. 1090 3. 6548 X 10" 7 0. 0780 6.0978 X 10'5 
1% 15. 6090 1. 2539 X 10" 17 15. 7340 1. 5282 X 10” 
25 188. 5940 7. 2213 X 10" 25 186. 9840 1. 1448 X 10” 


#217 Ma” H MWIUa 方法 的 迭代 残 量 


N(e—-1) T Res N(e—1.2) T Res 
7 0. 1090 3. 6548 X 10'? 7 0. 0780 6.0981 10!5 
17 15. 2190 1. 2513 X 105 17 15. 1870 1.5269 X 10^! 
25 186. 1880 7. 2213 X 10” 25 189. 3750 1.1438 X 10 


由 数值 结果 可 以 看 出 ,在 这 里 ,AlIU 和 MWIUa 方法 比 其 他 方法 更 有 效 . 在 相同 的 
计算 条 件 下 ,我 们 进一步 比较 了 这 两 种 方法 , 见 图 2. 2. 这 几 幅 图 给 出 了 数值 残 量 随 着 迭 


代 次 数 和 CPU 时间 的 变化 情况 . 这 里 实 线 表 示 MWIUa, 虚 线 表示 AlIU 方法 . 
5 5 - x 
0 0 
-$ -5 
-10 -10 1 
-15F. -15 
-20 -20 
-25 -25 
-30 -30 
-35 -35 
-40 EM -40 
-45 -45 
0.1 02 03 04 05 06 07 0 20 40 60 80 100 120 140 160 
SUB CPUN VSEE ( 取 对 数 ) SAMa Ib i OPE SERIA ( 取 对 数 ) 


图 2.2 四 种 情况 下 ALU 和 MWIUa 方法 的 比较 


23; 


mim 
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0.30 T 


全 

Ss 

a 
Sao 


~0,05 H+ +» +» M —À 
0 0.1 02 03 04 0.5 0.6 0.7 08 0.9 1 


参数 取 o 名 时 CPU 时 间 VS 和 迭代 误差 ( 取 对 数 ) 


0 20 40 60 80 100 120 140 160 
SRM WERE VSR (HOBO 


5 5 

of; 0 1 
-s[ -5 1 
-10r -10 1 
EE -15 1 
-20 -20 1 
-25 -25 1 
-30 -30 1 
-35 -35 
-40 -40 
-45 -45 


0 01 02 03 04 05 06 07 
AS tpa P I CPUIS SI VER (SURE ( 取 对 数 ) 


20 40 60 80 100 120 140 160 
5908 UTE OP V SE (IURE 〈 取 对 数 ) 


5 + + + T + 5 T T T T T 
or 0 
E -5 
-10| -10 
-15f -15 
-20 -20 
-25 -25 
-30 -30 
-35 -35 
-40 -40 
-45 -45 
0 02 04 06 08 10 12 14 0 20 40 60 80 100 120 140 160 
SAR IE CPURSE ISI VS (CURE ( 取 对 数 ) SRM GL! WRG MVS 〈( 取 对 数 ) 
图 2.2 (8) 


2.4.3 结论 


由 经 典 Richardson 方法 和 Marder-Weitzner 方法 ,建立 了 4 种 有 效 的 修正 Marder- 
Weitzner 方法 ,应 用 六 阶 紧 格式 增 量 未 知 元 更 加 有 效 . 

根据 数值 结果 ,MWIUa、AlIU 和 Al 方法 比 其 他 的 修正 的 Marder-Weitzner 方法 有 
Ak. 如 果 和 迭代 之 前 给 出 松弛 参数 的 值 , 每 种 迭代 方法 的 CPU 时 间 与 计算 机 计算 能 力 有 很 
大 关系 . 当 松弛 参数 随 着 每 次 迭代 自动 计算 ,收敛 速度 加 快 ,CPU 时 间 降 低 . 
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一 类 反应 扩散 方程 的 多 层 分 块 类 
小 波 增 量 未 知 元 


3.1 多 层 分 块 类 小 波 增 量 未 知 元 


考虑 下 面 的 反应 扩散 方程 ; 


a ay , a " 
u fe ut JUR Os (ay EQ, t€ (.T] 


EP FP REN 
uCr,y.0) = uo: @y EQ, t0 G.D 
u(xyyst) = 0, (x,y) EIN, t€ (0,T] 
其 中 ,vy 二 0,0<e<1,0Q 是 R" 中 的 一 个 开 集 ,并 且 边 界 足 够 光滑 ， 
g(s) = Mw. baa > 0 (3.2) 
求解 方程 (3. D E 一 种 方法 是 在 一 致 网 格 上 将 方程 离散 ,然后 构造 类 小 波 增 量 未 知 元 


(WIU) 型 数值 格式 . 但 是 数值 实验 中 ,我 们 发现 ,使 用 这 种 方法 得 到 的 稀 玖 矩阵 在 取 值 
很 小 的 情况 下 ,条 件数 增长 很 快 ,因此 ,使 用 WIU 方法 求解 上 述 方程 很 浪费 时 间 . 

为 了 更 加 有 效 地 求解 方程 (3. 1) ,我 们 提出 了 一 种 基于 和 矩阵 分 块 的 新 的 增 量 未 知 元 : 
分 块 类 小 波 增 量 未 知 元 (WBIU). 在 引入 新 的 向 量 范 数 后 ,我 们 建立 了 新 的 范 数 估计 和 
两 种 WBIU 型 数值 格式 : 显 格式 和 隐 格 式 , 并 证 明了 在 新 的 范 数 估计 下 ,格式 是 稳定 的 . 
最 后 ,数值 实验 证 明了 WBIU 的 有 效 性 . 

下 面 , 就 来 仔细 地 分 析 研究 这 类 反应 扩散 方程 的 多 层 分 块 类 小 波 增 量 未 知 元 . 

为 了 简单 , 令 Q= [0,1]:.， 现在 ,考虑 方程 (3. 1) 的 有 限 差分 离散 ,选取 步 长 h= 
1/(24N 十 1) ,NEN. 有 


or FXA,U;, 0 g(0,,) — 0 (3.3) 


其 中 U,, ERO? 是 在 节点 处 的 近似 值 组 成 的 向 量 . Au, 是 阶 数 为 24N)? 的 和 矩阵, 下面 
id: Aa=An, ,Us 一 Un 
B, =i ale =€ 
ii-i E E =ý 
Ae ig ho Road m 


=f B, —e 2+2€ 


UP6- yg kART GEN EE 


沿 工 方向 将 未 知 元 排序 ， 


U,— . |» uj = Cady udp oes? sud )T 


[S 
根据 矩阵 分 块 ,可 将 式 (3. 3) 重 新 写 为 


dug 
2t 


对 于 7 一 1,2,…,24N, 定 义 分 块 类 小 波 增 量 未 知 元 (WBIU ) 为 


FvC— uda + Buus — uda) + gu) — 0. B= 1,2,0,2 N (3.4) 


y = (ug uL/2 Ga 
z% = (uh. — u% )/2 i 
反 过 来 ,有 
uj, = yi — za "Y 
[m = yj Tua i 
引入 下 面 记号 : 
Usven 
U, = 
K ps | 
其 中 
ui ui 
ui us 
Use = ‘ 和 Usa = 7 
U2dN [M 
fil 
= y* 
v (| 
其 中 
y zi 
yi yi . ee zi 
Yen Bly 
容易 看 出 


Us = PjSU,. U, = SU, 
这 里 ,Ps 是 阶 数 为 (24N): 的 排序 矩阵 ,有 下 列 形式 ， 
P, = (Ea Vi OL 
Ey 和 Va 都 是 阶 数 为 2N X 277 N 的 矩阵 


0 0 0 0 0 1 0 0 0 
1 0 0 0 0 0 0 0 0 
0 0 0 0 0 0 1 0 0 
E, =| 0 y 0 0 0], V,2|0 0 0 0 
0 0 0 = 0 0 0 0 Q0 
0 0 0 * 0 I 0 0 0 + 0 
Ss 是 阶 数 为 (O^ NO? BS eO E. 
La oa 
Sa = QI 
i bay ^ 
I 是 24N 阶 的 单位 阵 . 因此 我 们 有 
Ua = PSU 
将 上 式 代 入 方程 (3.3) ,并 且 在 方程 两 端 同 乘 以 CP4S 0 , 则 成 立 : 
a T U = E 
TPS SU VOS "AP aS Ug + CS, g P,S,U,) = 0 


注意 到 PIP = LG, SIS, —21,C)I, P, Mg 可 交换 ,所 以 有 


à 


2 


Ts LYCOS, AP,SU, + Sig (SDs) = 0 


我 们 得 到 2 WBIU 型 格式 . 
类 似 地 ,在 YY 上 引入 下 一 层 WBIU. 4 


sic ya . c yo 
VAa za 


其 中 
yo wl 
4 4 
ye ys z^ 26 
YEN ZHA-. 


与 方程 (3. 6028 M, MF j—1,.2,7:,2^?N, 定义 

yG = yG ip 

bs. = y tus 

因此 

PPS 
其 中 Prim Err Vri Ol 是 22 N? fr AS HERE AR E 

li =f 
Sai = ( Je I, 


lj; li; 


这 里 ,Po-: Sa PH Pa ,Su 结构 相似 ,但 是 阶 数 为 (2 ND? 


* oo 0 0o 
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(3. 7) 


(3.8) 


(3.9) 


(3.10) 


(3.11) 


U28- yý kiA TAAR RNR JE k 


NE i 0 ) Saz Sai 0 
= o Lig "= 0 ha GIL 


其 中 Py ASE (24N)2 阶 矩阵 ,我 们 看 到 


BAPL-LGL. $4812 p " ) 
0 Zian, 
和 
U: = PS aUa (3.12) 
将 上 式 代 入 方程 (3. 90 ,两 端 同 乘 以 PaSa)" ,可 以 得 到 
2 es -FvSTA,SU, , +STg(SU, 1) 一 0 (3.13) 


xt .S—P,S,P oS a1. 
一 般 地 ,对 于 /二 d 一 1,d 一 2,…,1, 定 义 


= Y' 
y" = (5) (3.14) 
很 明显 
y = ps, Y? (3.15) 


这 里 ,已 ,S, 是 27 N? ORE ep. 
P, 一 (E, Vj) C I 


I1 =k 
S,— OL 
I, 


lug ig Y^ —U,.Y ^? =U, , 则 方程 (3. 15) 对 于 =d 也 成 立 . > 


Y! 
z 
l gm 
U, = Ja à 其 中 Y! € 243 N? Al € 243 y? 
1 
Jg z’ 
然后 有 
Un = PSU, (3. 16) 
其 中 


E E 0 z Sin 0 
"U (0 L@L) '" lo 4LGIL 
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其 中 上 = (24-2) N, 因 此 
~~ (LGI 0 E LGI 0 
Pêr- |" ‘ ) sr = 2(" E ) 


0 Oi, 0 LOL 
依次 取 /二 d 一 1,…,1 将 式 (3.16) 连 续 代 入 (3. 90 ,我 们 得 到 4d+1 层 WBIU 

2 Oth | STA SU + Sg (SUL) 一 0 (3.17) 
、 ~~~ ~ ~ ~ 1 1 
这 里 ,$=P,S aP 1S4 1,…,P1S1,Y=Y! xn z=(2 Z O zz 

"ry o5 P4 
3.2 逼近 格式 及 其 等 价 形式 
车 S=PsS4, 方 程 (3.17) 等 于 
aU = = 
S'S “Us + ySTA,SU, +Sig(S.U.) = 0 (3.18) 


经 过 计算 ,有 
sg) 二 p Z1) 
EO OU Z|) 


BH OUZ |?) A OC|Z4 |) ,方程 (3.18) 等 于 


20 wy 
2 dt OSIPA us + 2|" 0 J-» (3.19) 
若 S—P,S P 4S4 ,方程 (3. 17) 等 于 
sts Ue + ST AUS. +S" ST) = 0 din 


利用 上 面 的 结果 ,得 到 


2U qe) 
2? To + STA SU, + 2? ( = 


0 


-5 
x 
3j 


ST gCSU,. 1) =} PE SIP Tg (PuS uP a Ses iUas) 
—SI,PI,SIPIg(P,S Uy) 


=S P} Sig (Say) (3.21) 
ARABS DLE EE TEA Sy. ,P。; 的 定义 ,有 


E Si; (Ph 0)(g(¥*) 
Sg SU) -| i ( H ( ) 
o o IA 0 


- 0 (v | 
0 I 0 I 0 


m pu me » adi ) 
0 0 
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这 样 ,(d 十 1) 层 WBIU iE SCA : 


3 [Y = (Y) 
2 B (7) +vSTA,SU, + 2* g } =0 (3. 22) 
其 中 YE RNX2N ,ZE RO DNX2N. 

对 时 间 进 行 离散 ,我 们 提出 两 种 格式 : 
。 格式 工 : 显 格式 

2: (YH —Y* =. gY") 

= = DL) esaet: 42: ( ó ] =0 
* RR: 隐 格 式 

2 (Y? —y" "T gY") 

s m E --vSTA,SUT" + 24 ë =0 


4 Vy Gk VO SE th JE ROC w m GO + M= Chas jh) 张 成 的 空间 ,aa w(x) 在 单元 K = 
Lihas Go DA) X jh GEDA PEF 1E Hf EF 0. 令 ur GO Ji N V, 中 的 
步 长 函数 ,对 于 VxEK5 i 1.2. 24N  j 1,2, 2N ,有 un G0 uy, (M). BORE 


24N 2N 


un (x) 一 >) D>) ui, MD wn, OO 


i-1 j=l 


引入 下 面 两 个 差分 算 子 
Vin, Bd = GG lue — 9G) 
d 
Vow B(x) = db + her) — B(x)) 
l 
其 中 ,ei 二 (1,0) ,es 二 (0,1), 定 义 离散 的 向 量 内 积 : 
(Cun, Un) iy = CVn Un, Ving Vh IF (Von us, + Von, vn ) 
2 
Jtr CH, 0 EL CO PHL. S| e luu = (C 0: ,容易 看 出 上 ， | AM] + | 是 
Vi 中 的 Helbert 范 数 . 
将 差分 格式 (3. 3) 写 成 变 分 形式 
a 
( EO iu d E VRE (3.23) 
at d d d 
引入 记号 Mia (aha Bha) ,选择 均一 wu GO ,可 以 得 到 格式 (3.3). 定义 Y* 是 由 基 酉 
BEV on ,ms 《x) 张 成 的 空间 ,和 zm, 在 单元 [ihasGiTDhs)X[G2j 一 Dha,(2j 十 Dhs) 上 
等 于 1, 在 其 他 地 方 为 0. 因此 ,对 于 Vy € Y^ ,都 有 


24N 241N 


y(x) = >) 2 Y Miz) Va > LEN 
i=] j= 
定义 Z FE h HE PRL 


Xha My oia GO = w 
TK AR AY H, Alte p FT V 2€ Z^ ,都 有 


2îN 24 IN 


ZG) = 2 >) Moa Xu Os 2 EO 


i=1 j=l 


— w, 
a Miji UM M, oj 
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对 于 VsEZ V y! CY" Hi Vsus, OO 和 xi,,m.s_,(*) 的 定义 ,我 们 有 三 个 结论 : 
PE doe 
(1) E dr = 0. 


(2) | ztytdz — 0. 
a 


(3) Y4.Z4 实际 上 就 是 式 (3. 5) 所 定义 的 WBIU 空间 . 
由 结论 (3) ,空间 Va 可 以 被 分 解 为 
V,-Y' Oz 
IF VY u EV A 
ug = yt + zts y E Yt, E 

利用 上 述 分 解 , 我 们 有 下 面 的 定理 : 

定理 3.1 令 M=(ihs,2jha),Ai=(ihas (2) —1)ha) HF i=1,2,-°,2°N, j=1, 
2, ,27 1 N24 y — Vy, (0) Z= Xna GO. 时 , 变 分 形式 


ayi ~ = " » 

E 3)+ (FHT HEO = 0, V y EY (3.24) 

az 

(& 4) v ((I Hz Z) -0VzeZz (3. 25) 
也 就 是 式 (3. 19). 


证 明 ; (1) 对 于 ;一 1,2,…，,24N ,一 1,2,…，,24 1N 
4 F= V, GO ,可 以 证 明 式 (3. 19) 的 第 一 部 分 的 逐 点 形式 为 


dye, 
2 ith SEC Dyas — eyta + 2(2 + 2e) ly — Dey — hays 
E d 
zio t zaps Viaje) t 2g (yt2;) 0 (3.26) 


下 面 只 需 证 明 当 y — Vo, CY, FE PTE SK C. 24) EIE (3. 26). 事实 上 ,有 


d. 4-1 
ay! ~ SS A ay! (Mia) ayia 
(x .了 = P EE Va, ai OO Va, GO )— 275A (3.27) 


下 面 , 计 算 vCCy! nz" 00, 的 值 . 引入 下 面 的 记号 : 
Ay = (ua Gc + haer) s Von a Gc + haei)) ,As = Cua Gc + haer) «Vos ac OX) 


A; = Cua Gc + haez) Vo, Gc + hae2)) Ag = Qua Gc + hae) Voy Or) ) 
A; = (tg GO V, a Oc + hae )) Ay = (ug GO 5 Boi, OX) 


A; = Cua OO s Porm (X + hae2)) As = (tg GO Vj C) 
利用 步 长 函数 定义 ,得 
A, = Ay = As = A; = yh Az = yshi,As = yizhi. 
Ag = (yiz + yira — Ziga + che has 
A; Ciz + Yije + zhe — heja AG» 
AEI e a BERL O H EF 
Sj Cua Qc + haer) — ua (x) Woy, m Qc haer) — Voas, OO) + 


hi 
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gz Gta æ + hee) — ua Cx) Vy, Gc hues) = V», M OO? 
d 


以 2y£sj 2ey hao b 22 + 2e) yds 2eyt ao; XXL 好 ai t 22a ease) 


(3. 28) 
最 后 计算 
| 
gG' = (BOS) 3 y! Mia) Wo, an, GO) PaM) = 28 hah 
k-1 [-1 
(3.29) 


HRG. 27) . 式 (3.28) 和 式 (3. 29) HAIN ,可知 变 分 形式 就 是 式 (3. 26). 
(2) 可 以 算出 式 (3. 19) 的 第 二 部 分 的 逐 点 形式 为 


22 ua ii yioo + feja — eztia +2(2+ 26) 13; — erfa + 22234 + 
bua F Yiz) 一 0 (3. 30) 
下 面 只 需 证 明 当 Z— ja, GO 时 , 式 (3. 30) 就 是 式 (3.25). 事实 上 ,第 一 项 等 于 
(¥.z)= BAK. (3.31) 
第 二 项 等 于 
jz a Æ + haer ) — wae) qp, a, (E F haer) — nua GO) + 


PCIE 十 huez) — ua GO s Zapa, (X + haer) — Yaya, G0) 
d 


下 面 ,计算 等 式 右 端 部 分 . 引入 记号 ， 
B, = (ata Cx + hae) Xia (x + h461)) B; = Cua Gx + haer) + Xn ,.A, (x)) 
B; 
Bs 


Cua Cx + hae2) ,Xia Oc F haez)) ,Be = Cua Gx + haez) «yu, a, GO? 


Cua Cx) s Xanga, Gc + hae) ,Bs = Gua GO s Xaya, GO 


Br = Ga GO -Xa a, (XH haes)) ,Bs = Cua C0) s Xa a, QOO 


利用 步 长 函数 定义 ,得 
B, = 224 9;-1h3 Be = zhygahi.B. = 2241.21 h4 
B; = 2z2 2 hà B, 一 zZai ha s Bs 一 2225 hi 
Bs Cyto; 2 zea ziz )hi 


B; (— xa; ye — Ziaja — 2132 hà 

经 过 计算 ,第 二 项 等 于 

v(— Ryža- + zioa — ezta + 2(2+ 26) 214a — ezhi H Yim + zim) 

(3. 32) 

将 式 (3. 31) 和 式 (3. 32) 相 加 , 易 知 式 (3. 30) 就 是 变 分 形式 (3. 25). 定理 证 毕 . 

类 似 地 ,将 Y! ,1 二 d,…,1 作 如 下 分 解 

Y" -Y OZ 

就 等 到 了 多 层 分 块 类 小 波 增 量 未 知 元 (WBIU). 这 里 ,要 注意 的 是 Ye+ =V. 因此 ,对 于 


bog 


第 3 草 


Vu, € V4 RNA 

u= y+z, yEy EY, z€z—-Z QOZOQ OZ 
和 定理 3. 1 类 似 , 我 们 可 以 证 明 式 (3. 22) 等 价 于 下 面 的 变 分 形式 : 
Fz),5)20,V 


el 
< 
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dy ~ m 
(2.5) v (Qr 53)», + (OF 
(2.2) v(Q 62). =OVEEZ (3.33) 
3.3 关于 范 数 的 三 个 引 理 
证 明 稳 定性 分 析 之 前 ,首先 介绍 几 个 引 理 . 
引 理 3.1 对 于 任意 函数 g ,都 能 找到 这 样 的 两 个 正常 数 c 和 cs ,使 得 : 
g(s) > Lomi” a 
gG «2062s? ce. Ws 
引 理 3.2 il HRR, ITF 0d ,任意 函数 un € Vi. 
A420 +e) lu |< Mu, la, S Se CO + du, | (3. 34) 
其 中 S, D) Sha (e+ 4-44!) -2, 
引 理 3.3 任意 函数 VEY, 
S2(hi sha) | y IS «x y |l*, SeChisha) = huh. 
Sh sha) Mylau < yl Sho = — t | 
这 里 |。|- 表 示 最 大 范 数 . 
c > 
3.4 显 格式 和 半 隐 格式 的 稳定 性 估计 
定理 3.2 (格式 工 的 稳定 性 估计 ) 假 设 对 于 某 个 固定 的 m ,有 cmn. 
T 3 n 24 T Que) 
hå 161 +2624)” hä 4Ab Ms?” “ae 
其 中 
Mo =| ui, X 十 mcz) | Q| 
那么 ,对 于 任意 的 n 宇 0, 我 们 有 下 面 的 估计 式 成 立 : 
| uk, |? = |y" 1? + [2"1? <M (3. 36) 
WEBB. 容易 看 出 ,格式 工 等 价 于 下 面 的 变 分 形式 : 
(=) OHE HEO 一 0， WHEY 
(3. 37) 


T 


0, Vr €Z 


ME SN us a 
( z =) WOE Za, 


3ROBULUL UTRAM E SEN RUE 


在 方程 (3. 37) 中 令 y —2cy" E —2ez ,把 这 两 个 式 子 相 加 ,利用 关系 式 2(a—b.b) = |a]? 


— |b|*— |a—b|* ea 


D—|y"—y|*—|z" —z|*-2o | y +27 Hi t 2:GGD2 9") = 0 


SE Dn-|yt|*—]|y [2+ [2 | — |e" 2. 利用 引 理 3. 1, 得 到 


2 
Dn-|y"—y]|*—]|z?—2z|*-c2o || yz I DEEP Í, (y) dr s 2u, {a| 


ERG. 37) 的 第 一 个 方程 中 令 y mu"! 一 y"), 有 


| y — yt l* Ew CO! z^, — y», + rlg(y") y — y") 一 0 
利用 Cauchy-Schwarz 不 等 式 和 引 理 3. 3, 有 
[y ? —y'|* &w l y zu, Ey —» lla, te let ly? —» 
i 


<a =——_ | yte lla ly —y" 十 
Sia) i 


e [gone yty]? 


(3.38) 


1 
«1 jyt -yl te [gy |?+— 24 — | ye ls, 
2 |y | | S, Qi sha)? d 
因此 ， 
|y? — y |? < 27 |g( Dpt l y +l? (3. 39) 
FUCHS SEIT ERE C CEN i 
现在 在 (3. 37)» WPA für =r 2) 
| 27 — z^ | 二 rw((y Ha —2")),, —0 
使 用 引 理 3.2 4 1—d. 
gg NS n Lon Vm an 
| 2 I s cv ly" rz lla lle "la, Sg lo tel le | 
(3. 40) 
因此 
Jea Sor ly ten ll, (3.4) 


将 式 (3. 38) 5X3. 39) MRG. 40) 结 合 在 一 起 ,使 用 S, Ca A BRA SC ,我 们 得 到 


2 
D +P, Ce) lb yz | + tba | em “da < 2x, |A| +22 |gCy |? 


其 中 
4(e 十 1)z — 88 (2^ 十 1) 
Ti(e) = 20 m othe 
利用 (3. 35) 的 第 一 个 条 件 , 有 (Ce) > 0. 利用 引 理 3. 1 和 引 理 3. 3 ,我 们 有 


27 bi A n 
Ë EOD I < pas [nr f Oade e [a 


(3. 42) 


" 
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MAW Il yr z^ i, 21y z^. 
To 4-21 Ce) ly He |? + eb HR f. (y')*dzr <T; 


其 中 
2tbaa 


E QHD 27D 


由 O<eX1, nT pL fg s 


i= |y (7, Ty = rey +272) | 2 | 


DHD |y + I + rb A HD| oras c ns 


下 面 使 用 数学 归纳 法 证 明 (3. 3360 5X. 
CIO & 一 0 成 立 因为 |y% |? + | 2° |* Mo 
CI) 假设 (3.36) 对 于 =n HE REE |y |+ |e |! M, 
CIID XT kR=n+1. AACS. 35) 的 第 二 个 条 件 和 (3. 33) ,我 们 有 
ly |? — [yt |? + fe |* — [et [Pov dy te"? <P 
注意 到 [otter |?=ly" +l], 
|y? |* + |z? | << 1— 0C] y |? + | | Hr 


也 就 是 


|y? |? + fe |? <a — 007? CL y! Iz ID 4- TS Dw) 


i-0 


«a-o»? yr ||?) - ir, 


-ü-n"Qy |? + [2 [D - 3c 十 acz) lal 


<M, 
稳定 性 估计 (3. 36 TEE. 
为 了 证 明 格 式 开 的 稳定 性 ,首先 介绍 下 面 的 引 理 . 
引 理 3.4 (离散 的 Gornwall 引 理 ) 假 设 c0 和 1+cA>0.a" 6b" 分 别 是 两 个 非 负 序 
列 ,满足 


ntl 


an «br. Bh <b.Vn>S0 


那么 有 
^ l „a, l+a 1 
a < Fa t= (1 axum)^ Vn>0 
定理 3.3 《格式 下 的 稳定 性 估计 ) 假 设 对 于 某 个 固定 的 mrsm ,着 
T 240 
AA = 25, MT (3. 43) 
其 中 


(2c 十 rocz) | Q | 


1+4 
M, = ji, | +S 
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那么 下 面 的 估计 式 成 立 : 
la, I — |y P+ le <M. 对 所 有 的 n 之 0 (3. 44) 
证 明 : 要 证 明 格 式 开 是 稳定 的 ,只 需 证 明 下 面 的 变 分 形式 ， 
Gs eornm +275), GOD) =0, Vy €Y. 
+ 
gH. ds (3. 45) 
(S22 eon Ta,2),—0, VECZ 


T 
是 稳定 的 . 在 方程 (3. 45) 中 令 y 一 2ry” E= 2 ,然后 将 二 者 相 加 ,利用 208 b,a) = 
lal*—|b|*-- |a—b|* ,得 到 下 式 
D + ly"? —y" |? + lot — 2" |? +20 y yn ann i, 


=— 2r(g(y") y"! — y!) —2r(g CQ y») 
= Sr(g(9), a" — yn) — vba, Í, Cy? dz + 2x, |A| 


ixH DS|y H| ly Hetl 3838 E 5x18] 
D + |y — y|? + |z — z" |? +20 Il yt +z” i, + 


Thri Í, (办 )2dz 十 2rcl |A| &—2c(gCy y" — y") +2rcı | Q| (3.46) 
由 于 (3. 42) 和 
—2r(g(y) y — y') «2c | gly") Il y?! y" | 
x|»"—»[l-c-zàlbgn |? 
有 下 式 成 立 : 
Dy |e" —2"|? +20 || y"? 4- 27? || hc 
Tba 1 -r> Í, Cy dx < 2rci la| +rřo lnl 


利用 引 理 3. 2 和 条 件 (3. 43) ,我 们 得 到 
Do 4v |y? 27 | < 2r |A| +c |a| 


因为 [yn eem [1 Lyr [24 [en |? E MAS ps 
Pe p asc y |? + [e | « 25 [0| +e? 1| 
使 用 离散 的 Gronwall 2| E nT Ail 
ele er BSE com m m) Ca tae Lal 
因此 ， 
lx ine pep < EXER ISI EIUS (2o trce) [0] 


稳定 性 估计 (3. 4140 vr. 
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3.5 数值 结果 


考虑 下 面 的 方程 
9 a 9? 
= G JZ - =} w= f. Gy) € Q.t€ (0,1) 


uCr,y,0) = (1—3)sin(xy)(e&7^7 —1), (r,€Q0.4€0 
uCr.y.0) 一 0 (x,y) € 9N,t E (0,1] 
其 中 
f =Q —2)sin(ry)e! Ce — 1) — (esin( xy) e^ + (1 — z)sin(xy) e^7^ ) + 
(1 — x)sin(xy)z! e (eo — 1) + C(1 — zx)sin(xy)e (e 一 1))3 
Q= (0.1)? ,方程 的 精确 解 是 (1—3)sin(zy)e' (eE — 1), 

分 别 使 用 WIU 和 WBIU 求解 上 面 的 方程 . 我 们 发 现 , 当 。 很 小 的 情况 下 ,使 用 
WBIU 方法 得 到 的 系数 和 矩阵 的 条 件数 要 远 远 小 于 使 用 WIU 方法 得 到 的 系数 矩阵 的 条 件 
数 . 表 3.1 和 图 3. 1 分 别 给 出 了 当 s 和 层 数 4 分 别 取 不 同 值 时 ,系数 矩阵 条 件数 的 比较 . 
从 表 3.2 MR 3. 3 可 以 看 出 ,使 用 WBIU 方法 的 CPU 时 间 要 远 远 小 于 WIU 方法 使 用 的 
CPU 时 间 . 图 3.2 和 图 3. 3 分 别 给 出 了 在 s=0.1,0.01,N=1 和 qd=4 时 ,使 用 半 隐 格式 
和 显 格式 的 误差 估计 比较 . 从 图 形 可 以 看 出 ,虽然 ,使 用 者 两 种 方法 达到 的 精度 基本 相 
同 ,但 是 ,基于 WBIU 的 数值 格式 节省 了 大 量 的 CPU 时 间 . 


表 3.1 系数 矩阵 条 件数 的 比较 


E 0.1 0.01 0. 001 
WBIU 23. 622 21.859 21.682 
WIU 706. 67 1247. 5 1357. 3 


53.2. 半 隐 格式 CPU 时 间 的 比较 


E 0.1 0.01 0. 001 
WBIU 115. 33 97. 391 83. 359 
WIU 371. 86 450. 98 453. 98 


53.3. 显 格式 CPU 时 间 的 比较 


€ 0.1 0.01 0. 001 


WBIU 1193 1218 1197 


WIU 2363 2393 2338 
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600 ij 

—+— WBIUs 
500 f —e— WIUs 
400 F 


图 3.1 e=0.001.r = l 


200 


0.035 
0.03 
0.025 
0.02 
0.015 
0.01 
0.005 
0 


1 站 > 7 
0.8 - 02 
0.6 - S 
04 02 24] x 
y 


€=0.1, t=0.5 


£—0.01, =0.5 


图 3.2 c—1/200.d—4.N—1 时 ,WIU 和 WBIU 在 半 隐 格式 下 的 误差 曲面 图 


0.035 
0.03 
0.025 
0.02 
0.015 
0.01 
0.005 
0 


09 08 07 06 
y 
£-0.1, 1 


99 08 07 06 05 04 03 02 
3 0 
y 


£-0.01, {=1 


时 WIU 和 WBIU 系数 矩阵 条 件数 的 比较 


05 04 03 02 "mi 


0.07 
0.06 
0.05 
0.04 
0.03 
0.02 
0.01 
0 
1 


5 
09 08 07 06 05 04 03 93 fo x 


第 3 章 一 类 反应 扩散 方程 的 多 层 分 块 类 小 波 增 量 未 知 元 039 


0.03 
0.025 
0.02 
0.015 
0.01 
0.005 


0 
1 


08 0.6 


01 0 02 st] x 
y y ý 


£—0.01, (=0.5 €=0.01, t=1 


图 3.3 


t =1/2000.d=4.N=1 R}. WIU 和 WBIU 在 显 格式 下 的 误差 曲面 图 


YT 
非 埃 尔 米 特 正定 线性 系统 的 迭代 方法 


4.1 简介 


在 科学 与 工程 的 很 多 重要 领域 中 ,经 常 要 求解 下 列 大 型 稀 玻 线性 方程 组 
Ax =b, x.bcc" (4.1) 

Hep A€ CH" JAKAR it P 38 A HEIL Wr JUR BEE 

对 于 这 种 大 型 稀疏 线性 方程 组 ,我 们 通常 采用 迭代 法 求解 . 24 R AUE IE A 为 对 称 正 
定 矩 阵 时 ,最 有 效 的 方法 之 一 就 是 共 力 梯度 法 . 而 对 于 更 一 般 的 对 称 系统 , Paige 和 
Saunder 提出 了 MINRES 方法 , 即 利用 对 称 Lanczos 基 来 计算 一 个 方程 的 近似 解 , 并 使 得 
残 向 量 在 Krylov 子 空间 上 最 小 ，1986 年 ,Y. Saad 和 M. H. Sehultz 进一步 推广 了 此 方 
法 ,提出 了 GMRES 算法 ,用 来 求解 非 对称 线 性 系统 .而 在 解决 实际 问题 时 ,常常 会 遇 到 
这 样 的 情况 : 尽管 迭代 矩阵 的 谱 半径 小 于 1 但 却 非常 接近 于 1. 此 时 定常 迭代 法 的 收敛 
速度 会 非常 缓慢 ,效果 并 不 好 ,甚至 不 收敛 ; 当 系 数 和 矩阵 的 谱 分 布 很 分 散 时 , 非 定常 运 代 法 
的 收敛 速度 也 很 慢 ,甚至 不 收敛 . 因此 ,办 代 法 的 求解 速度 经 常 被 线性 代数 系统 系数 和 矩阵 
的 性 质 所 制约 而 不 能 实现 快速 求解 . 为 了 解决 这 一 问题 ,人 们 提出 了 预 处 理 技 术 , 预 处 理 
技术 是 指 对 线性 代数 系统 的 系数 矩阵 进行 处 理 , 进 而 使 处 理 后 的 线性 代数 系统 易于 求解 
的 技术 . 经 常 采用 的 方法 就 是 在 方程 两 边 同 时 乘 以 一 个 非 奇 异 和 矩阵 的 逆 , 然 后 再 用 迭代 
法 来 求解 . 利用 预 处 理 技术 可 达到 使 定常 迭代 法 迭代 和 矩阵 的 谱 半径 变 小 、 线 性 代数 系统 
系数 矩阵 的 谱 分 布 更 加 集中 等 目的 . 近年 来 ,人 们 不 断 提出 各 种 新 的 预 处 理 技术 ,而今 ， 
预 处 理 技术 和 迄 代 法 已 密 不 可 分 , 预 条 件 技术 作为 加 快 迭 代 法 收敛 速度 的 有 效 方法 ,在 求 
解 线 性 代数 系统 的 过 程 中 扮演 着 越 来 越 重 要 的 角色 . 

本 章 在 已 有 的 成 果 基 础 上 ,基于 大 型 稀 朴 非 埃 尔 米 特 正 定 线性 系统 的 正规 /反对 称 分 
裂 (NSS) 方 法 ,提出 了 预 条 件 正规 /反对 称 分 裂 (Preconditioned Normal/skew-Hernitian 
Splitting, PNSS) 夫 代 方法 及 其 变形 和 两 参数 预 处 理 NSS 迭代 策略 理论 分 析 表 明 ,在 一 
定 条 件 下 ,这 些 和 迭代 格式 是 收敛 的 , 除 此 之 外 ,给 出 了 和 迭代 格式 中 参数 和 和 迭代 矩阵 谱 半 径 
的 最 小 上 界 的 计算 方法 . 在 数值 实验 中 ,选取 增 量 未 知 元 (IU) 和 对 称 逐 次 超 松弛 
(Symetric Successive Over Relaxation,SSOR ) 两 种 预 处 理 和 矩阵 . 数值 结果 证 明了 收敛 定 
理 的 正确 性 和 方法 的 有 效 性 . 


第 4 章 非 埃 尔 米 特 正 定 线性 系统 的 选 代 方法 


4.2， 非 埃 尔 米 特 正 定 线性 系统 的 预 条 件 NSS 方法 


4.2.1 预 条 件 正规 /反对 称 分 裂 (PNSS) 和 迭代 方法 的 建立 
为 了 提高 NSS 方法 的 收敛 速率 ,我 们 将 NSS 方法 进行 改进 ,将 NSS 迭代 格式 用 于 求 
解 下 面 的 预 条 件 线 性 方程 组 
Ax =b (4, 2) 
Jep A-R C AR 8 =Rx.b=R*b.REC™" EK P48 HE HY Mh EAE ME, RT = CR)" 
是 R HIRIE RE EL. 386 IE P= RR. AS SUB FE B/S ERAT 8 PNSSO REX. 
PNSS 迭代 方法 : 给 定 一 个 初 值 x”EC", 对 于 k= 二 0,1,2,…, 直 到 {x ) 收 敛 , 计 算 


E 
(aP 4- NOx) = (ap —$)x® +b 
(4.3) 


(aP -- S)x*? = (aP — Nx CP +b 
其 中 ao>0,PEC"" 是 一 个 埃 尔 米 特 正定 矩阵 . 
注 4.1 从 迭代 格式 可 以 看 出 , 当 PI, I 为 单位 矩阵 ,PNSS 和 迭代 格式 退化 为 NSS 
迭代 格式 . 
下 面 我 们 重点 讨论 PNSS 方法 的 收敛 性 . 


4.2.2 PNSS 迭代 方法 收敛 性 分 析 


PNSS 迭代 格式 可 以 等 价 地 写成 如 下 形式 ， 
x" = M(a)x® +N(a)b, RE 三 0)1, 2 (4. 4) 
其 中 


M(a) (aP + S)! (aP — N) (aP + N) (aP — S) 

N(a) = 2a(aP + S) ! P(aP +N)” (4. 5) 
这 里 ,M(a) 为 PNSS 迭代 的 迭代 矩阵. 由 经 典 的 迭代 理论 可 知 , 要 证 明 PNSS 3E XE 
收敛 的 ,只 须 证 明和 迭代 矩阵 的 谱 半径 oM CODO 1. 

由 于 S HIRERE, P SSR R") S=R? (CR) SR OR, AA PS 所 有 的 
特征 值 均 为 虚数 . 设 O, 为 矩阵 PS 特征 值 的 集合 ,P71S 的 第 j 个 特征 值 为 ie Cj— 1. 
2,…,n) i= /—1. 我 们 可 以 得 到 PNSS 迭代 方法 的 收敛 性 定理 . 

定理 4.1 设 4EC"” "是 正定 矩阵 ,NEC"” "是 正规 矩阵 , SEC” "是 反对 称 和 矩阵 ,并 且 
满足 4 一 N 二 SPEC”" 是 埃 尔 米 特 正 定 矩 阵 . OX PN 特征 值 的 集合 . a 二 0 为 正常 
数 . PNSS 方法 迭代 矩阵 的 谱 半径 则 满足 下 面 的 不 等 式 

e(M(a)) <ola)< 1, Vac0. (4. 6) 


La—A,; | (=y) +45 
= E | i x 
dde etal ae A EU mon 


这 里 ,i 一 V 一 1. 即 由 PNSS 方法 得 到 的 迭代 序列 {xz” } 收 敛 到 线性 方程 组 (4. 1) 的 解 x*. 


其 中 
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增 量 未 知 元 及 其 预 处 理 迭 代 算 法 


BR EZ Pb s Vmin RI. Ys Menin 和 mx 分 别 表 示 PN 的 特征 值 实 部 绝对 值 的 最 小 值 和 最 大 值 ， 
虚 部 绝对 值 的 最 小 值 和 最 大 值 ,.Q 二 [Yon t Vimar ] X tni Tha ] ,有 


pas 2 | 42 
a” = arg min{ max eee 


Aui Voss — Tuas? han d Viia Yea 
o.»eaN (a+ 3)? + 7° | 


A Yain F Brux . Tmax = M Vmin Ymax 


4. 8) 
相应 地 ， 
ye E 
min 十 Ymax 一 2 Yminmax 一 Nmax yay 
[2 i yal Ta | , Ws < min Ymax 
Yosin + ee de sy nx — max 
ola") = (4.9) 


Yao ase = Yun |. 
[ee rm = i 
WEBB. 由 和 矩阵 范 数 的 相似 不 变性 ,有 
oCM(a)) =pl aP — N) (aP + N)™ GP — S) aP -- S) 
< || C&I — P? N) (aI + P? N)?)(a1I — P?^S)(al 4-P^?S)?) |l; 
< || C&I — P?! N) (aI 4- P?! N)?) |l; || (Ql — P?! S) (aI 9- P!S)?) |l; 


Toa: = VY min max 


MAH || (aI — P S) (aI HPS) || ;— max || —1. 因此 
ie €@, |a +iej 
e(M(a)) < || (al — P^ N) (aI + P> N)? |l; 
=å a— (Y; dig) 
< Wy _ aT AY T j? 
eo, a Àj ort Co, at Oi; d ig) | 
— (a — Y;)* +95 


E Cat y +95 
由 于 PON 为 正定 矩阵 ,可 知心 的 实 部 y, > 0. 又 因为 >0, 易 看 出 o(M(a))<ola). 
下 面 计算 使 得 PNSS 迭代 矩阵 谱 半径 上 界 (a) 取得 最 小 值 时 参数 a 的 值 . 


ixl Ric DT. 为 关于 了 的 单调 递增 函数 ,我 们 有 
aty)’ Fy 


mar |S = ae D Ly 三 max Ca — aa Toes 
peaN Cat)? HN Tainta N Ca +)? + sa 


(a — Yun)! + Quas [Ca — Yeux). F Pina Tpit. 
m . ?Max 一 NV Yuin Y max 
23 NC + Yu)? + D M lat Ya 十 oe | 


Ca — Yasin)” Ye ， Lm 
Ca F Yain)? + ax Nmax — V Vmin Ymax 


如 果 a* 是 ola) 的 最 小 值 点 ,那么 我 们 有 下 列 结论 成 立 : 
A) Hi Numar < V Vmin Ymax 时 ,有 


(a = Ha) FE sx es + as 
(4.10 
(Gate la F Your)? + SAN 


解 式 (4. 100 ,得 到 a” = VY min Ymax — T] max 


BAF 非 埃 尔 米 特 正定 线性 系统 的 迭代 方法 


(2) 当 Mma V Vmin Ymax 时 ， 


d ( (a — Yuin)? + qus 
erm : ne 


解 式 (4. 1D ,得 到 a* = Sin pu. 
综 上 所 述 ,我 们 有 


. ld Tax C. A Yuin max qud 
w= ; 
VN Wan. Ws /Nai 
die ca ER m 
LE RA , [Fe Teed Tie ET 
z7 ME 
min | /max min/max ^  ?]max 
Tuin + Yuax + 2. A/ Juin max — inex 
oa") 一 (4.13) 


$ 1. we RES » T 
| ET me | 和 
ma T Ja in 


定理 证 毕 . 

注 4.2 对 于 预 处 理 和 矩阵 忆 的 选取 ,首先 应 保证 它 是 对 称 正定 矩阵 ,其 次 ,尽量 使 得 
PNSS 迭代 和 矩阵 谱 半 径 越 小 越 好 . 在 数值 试验 中 ,我们 选取 了 增 量 未 知 元 和 SSOR 两 种 
TUAE FOR PE ,给 出 了 和 迭代 矩阵 谱 半径 在 预 处 理 前 后 的 比较 . 


4.2.3 不 精确 的 预 条 件 正 规 / 反 对 称 分 裂 (IPNSS) 和 迭代 及 其 
收敛 性 分 析 

在 PNSS 迭代 过 程 中 ,每 一 步 和 迭代 需 要 求解 两 个 子 线性 方程 组 . 当 系 数 和 矩阵 为 大 型 
稀 玖 矩阵 时 ,如 果 要 精确 算出 它们 的 解 ,费时 费力 . 为 了 提高 PNSS 迭代 格式 的 效率 ,我 
们 对 这 两 个 线性 方程 组 进行 迭代 近似 求解 ,并 把 这 种 迭代 方法 称 为 不 精确 的 PNSS 
CIPNSS) 38 Jr ik. 

不 精确 的 PNSS(IPNSS) 和 迭代 方法 : 给 定 一 个 初 值 x*“” EC", 对 于 ==0,1,2,…, 直 到 
{中 ) 收 化 ,以 5 外 为 初 值 ,利用 某 种 办 代 方法 (比如 CG 方法 ) 近 似 计算 (aP 十 NE OE) 之 
(aP —S) x +b BU fg. IRE x: C) LER JD x C10 为 初 值 ,利用 某 种 迭代 方法 (比如 
Krylov 子 空间 方法 ) 近 似 计算 (aP 十 S)x **" ~ (aPN) x C3) +b YAR EB x 6». 

为 了 简化 数值 计算 和 收敛 性 分 析 , 将 IPNSS 迭代 格式 写成 如 下 等 价 形式 : 

给 定 一 个 初 值 XEC', 对 于 k= 二 0,1,2,…, 直 到 {x 中 ) 收 敛 . 

CD HEADER RAR HET FBR CaP HNZ er” =r =b—Ax® HARE pe =r? — 
(aP--NDz ^ il || p^ le lr || ,然后 计算 run +) 一 元 十 去. 

(2) E BUR E PETI REAL PHOT? ar OD qp OU a AUD ) 直 到 残 
ga UD =O) atoz wie | GED pec FOP | ,然后 计算 m = 
x(t B E (#4), 

利用 和 矩阵 分 裂 
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增 量 未 知 元 及 其 预 处 理 迭 代 算 法 


A=M, —N, = (aP + N) — (aP — S) 
M; — N: = (aP +S) — (aP — N) 


引入 下 面 范 数 
Ill x Hla = Il Mx Il. 
其 诱导 范 数 为 
lll X lll = || MXM“ |. 对 于 所 有 XE Co 

这 时 , 若 令 M=aP +S. WES) FAY IPNSS 收敛 定理 . 

定理 4.2 令 4EC"*”" 是 一 个 正定 矩阵 ,NEC"" 是 一 个 正规 矩阵 ,SEC"x" 是 一 个 反 
对 称 和 矩阵 JH HW E ASNES. PE C" "是 一 个 埃 尔 米 特 正定 矩阵 . x' CC 为 线性 方程 
组 (4.1) 的 解 . 车 {x 中} 为 IPNSS 方法 得 到 的 迭代 序列 , 则 有 下 式 成 立 : 

||| P — x* ||| Cola) + om) (lt) Il x? —x" Hl. £0.52. 

FEB Il xl] RÆ Ill xe lll apts yo= || CaP +S) CaP3- ND * I2.0— || A&P +S)" || 2. 

特别 地 ,如 果 (c(a) + Opus? Cl F Oemar) <1 JE A EAR AE Ce? ) 收敛 到 线性 方程 
组 (4.1) 的 解 x* XE sem = max {e } Hes = MAX { Me m 


4.2.4 数值 算 例 


例 4.1 求解 下 面 的 反应 扩散 方程 
— (us d uy) bu. — f» Gy EN 
bes. (x+y) € AN 
其 中 020 是 一 个 常数 ,Q 二 (0,1)*. 使 用 网 格 步 长 为 的 五 点 离散 格式 : 


1 b . 
ni Gui ui UiH,j Ui j- Ui jl ) + 2i Cuna uia? fu , 


这 里 uy FL fy DIRIR u 和 了 的 近似 值 ,得 到 下 面 的 线性 方程 组 
AU — f (4.14) 
和 和 矩阵 的 HS 分 裂 不 同 , 对 于 给 定 的 矩阵 4, 和 矩阵 的 NS 分 裂 不 是 唯一 的 ,这 里 选取 


N LC FAT S s LA Ai-sa S, Sia; 


其 中 i= VaI, 是 一 个 实数 . 

为 了 有 效 求解 线性 方程 组 ,我 们 采用 增 量 未 知 元 (CIU) 预 处 理 和 对 称 逐 次 超 松弛 
(SSOR) 预 处 理 矩 阵 ,分 别 从 和 迭代 矩阵 的 谱 半 径 .和 迭代 步 数 和 迭代 时 间 三 方面 对 NSS 方法 
和 PNSS 方法 进行 比较 . 


1. NSS 方法 和 PNSS 方法 迭代 矩阵 谱 半 径 的 比较 


选取 c—3.0—1 和 10, 分 别 用 NSS 方法 .NSS(IU) 方 法 和 NSS(SSOR) 方 法 进行 计 
算 . 由 定理 4. 1 可 知 ,无 论 NSS 方法 还 是 PNSS 方法 ,都 有 nox > I Yus Ymar ,所 以 有 a= 


VYmin 十 Wmax， 表 4.1 和 表 4. 2 分 别 给 出 了 三 种 方法 在 N 取 不 同 值 时 参数 a 和 谱 半 径 
PCM(a)) 的 比较 . 


RAK 非 埃 尔 米 特 正定 线性 系统 的 选 代 方法 


表 4.1 b=1 时 ,参数 a 和 谐 半 径 op (Ma) MLR 


NSS NSS(IU) NSS(SSOR) 
a p(M(a)) a p(M(a)) a p(M(a)) 
3 3.0475 0.9141 3. 2208 0. 8244 3. 0606 0. 8996 
4 3.0154 0. 9500 3.1370 0. 8599 3. 0397 0. 9138 
5 3. 0064 0. 9674 3.0740 0. 9038 3.0255 0. 9310 
6 3. 0031 0.9772 3.0403 0. 9200 3.0166 0. 9444 

54.2. b—10 时 ,参数 a 和 谱 半径 p (M Ca )) 的 比较 

NSS NSS(IU) NSS(SSOR) 

N 
a p(M(a)) a p(M(a)) a p(M(a)) 

3 3.0475 0. 8810 3. 2208 0. 8066 3.0137 0. 9214 
4 3.0154 0. 8972 3.1370 0. 9463 3.0108 0. 9333 
5 3. 0064 0. 9292 3.0740 0.9141 3.0071 0. 9418 
6 3.0031 0. 9498 3.0403 0. 8726 3.0049 0. 9444 


2. NSS 方法 和 PNSS 方法 收敛 速度 的 比较 


选取 c—3.0—1 和 4 二 10, 分 别 用 NSS 方法 .NSS(IU) 方 法 和 NSS(SSOR) 方 法 进行 
计算 ,定义 
NSS 方法 所 用 时 间 
时 间 比 = PNSS 方 法 所 用 时 疝 
从 表 4. 3 和 表 4.4 中 可 以 看 出 , 随 着 网 格 数量 的 增加 ,时 间 比 逐渐 变 大 ,PNSS 方法 优势 
逐渐 增强 ,无论 是 迭代 步 数 ,还 是 CPU 时 间 ,PNSS 方法 都 要 优 于 NSS 方法 . 


表 4.3 0 一 1 时 ,迭代 步 数 和 和 迭代 时 间 的 比较 


NSS NSS(IU) NSS(SSOR) 
IT CPU IT CPU 时 间 比 IT CPU 时 间 比 
3 105 0.5938 47 0.1875 0. 1875 85 0. 1803 3.30 
4 179 1. 9844 62 0. 5313 0. 5313 106 0. 3705 5.34 
5 274 6. 6094 84 1.5625 1. 5625 133 0. 9313 7.10 
6 389 14. 188 113 4. 0625 4. 0625 166 2. 7540 5.15 
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表 4.4 b—10 Bd GARD BUND FCR i8] B Ee g 


NSS NSS(IU) NSS(SSOR) 
di LT CPU IT CPU 时 间 比 IT CPU 时 间 比 
3 62 0.5313 42 0.25 2,12 88 0. 2904 1.83 
4 94 1. 1875 52 0. 7812 3.52 104 0. 4106 2.89 
5 139 3. 7344 60 1. 5313 2.44 122 0.9514 3,93 
6 195 5. 2576 7l 3.0781 1.71 142 2.5637 2.05 
4.2.5 总 结 


本 节 研 究 了 求解 非 埃 尔 米 特 正定 线性 方程 组 的 预 条 件 正 规 / 反 对 称 分 裂 (PNSS ) 方 
法 ,理论 分 析 表 明 , 对 于 任何 给 定 的 初 值 ,PNSS 方法 都 收敛 到 原 线性 方程 组 的 解 . 除 此 
之 外 ,还 推导 出 了 和 迭代 矩阵 谱 半 径 的 上 界 以 及 当 上 界 取得 最 小 时 参数 a 的 计算 方法 . 在 
求解 每 个 子 线性 方程 组 时 ,为 了 节省 时 间 ,文章 给 出 了 不 精确 的 PNSS(IPNSS) 方 法 及 其 
收敛 定理 . 在 数值 试验 中 ,我 们 选取 预 处 理 矩 阵 忆 为 IU 和 SSOR ,取得 了 较为 满意 的 效 
AR. 今后 的 工作 可 以 考虑 将 PNSS 方法 用 于 求解 系数 矩阵 为 非 埃 尔 米 特 半 正定 矩阵 的 线 
性 方程 组 ,例如 求解 鞍点 问题 ,并 进行 相应 的 理论 分 析 . 


4.3” 非 埃 尔 米 特 正定 线性 系统 的 两 参数 预 处 理 NSS 
迭代 方法 
4.3.1 两 参数 预 处 理 NSS 迭代 方法 的 建立 


在 NSS 迭代 方法 中 ,每 迭代 一 步 需要 求解 两 个 子 线性 方程 组 ,并 且 这 两 个 子 线性 方 
程 组 的 系数 矩阵 aI +N fla --S 的 性 质 各 不 相同 ,因此 ,从 优化 的 角度 来 讲 , 我 们 可 以 采 
用 不 同 的 参数 , 即 在 NSS 近 代 格式 中 引入 两 个 参数 ,分 别 记 为 a 和 BB, 这 样 得 到 的 迭代 方 
法 称 为 两 参数 NSSCGNSS) 和 迭代 方 法 . 

GNSS 迭代 方法 : 给 定 一 个 初 值 x" EC", 对 于 k 二 0,1,2,…, 直 到 {x } 收 敛 ,计算 


(af + NOx 3) = (aI — S)x'? +b 
(4. 15) 


(BI -- S)x*? = (gr — Nox C? +b 
其 中 620,80, N € C" J TERE TE S € C" BAT PRAE P 
显然 ,GNSS 迭代 方法 是 NSS 迭代 方法 的 推广 , 当 a— t, GNSS 方法 退化 为 NSS 
迭代 方法 . 我 们 将 GNSS 迭代 格式 用 于 求解 预 条 件 线性 方程 组 (4. 20 ,就 得 到 两 参数 预 条 
件 NSSCGPNSS) 和 迭代 方法 . 
GPNSS 迭代 方法 : 给 定 一 个 初 值 xo EC", 对 于 k= 二 0,1,2,…, 直 到 {x 中 } 收 敛 ,计算 


(aP +N)x("?) = (ap —$)x® +b 


(4. 16) 


(BP + S)x**» = qp — Nox C? 45 
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其 中 az0.870.N€ C" 为 正规 矩阵 ,SEC” "为 反对 称 和 矩阵 . 

注 4.3 当 预 条 件 矩 阵 P= 工 时 ,GPNSS 退化 为 GNSS Wik. 

注 4.4 当 预 条 件 矩 阵 P— 1.2 0 a— it .GPNSS 方法 退化 为 NSS ERE. 

下 面 , 我 们 证 明 当 两 个 参数 满足 一 定 的 条 件 时 ,GPNSS 和 迭代 方法 是 收敛 的 ,并 且 
GPNSS 的 收敛 分 析 与 NSS 和 PNSS 方法 的 收敛 分 析 是 相 容 的 ; 除 此 之 外 ,我 们 给 出 了 使 
得 迭代 和 矩 阵 的 谱 半 径 取 得 最 小 时 参数 a 和 8B 的 计算 方法 ;最 后 ,我 们 选取 IU 和 ILU 分 解 
做 预 处 理 ,进行 数值 试验 ,数值 结果 证 明了 GPNSS 方法 在 理论 分 析 上 的 重要 性 和 实际 计 
算 中 的 优越 性 . 


4.3.2 两 参数 预 处 理 NSS 和 迭代 方法 收敛 性 分 析 


为 了 方便 ,我 们 将 GPNSS 和 迭代 方法 写成 下 面 的 等 价格 式 : 
xH = Mla, x? +N(a,B)b, k —0,1,2,-— (4. 17) 
其 中 


M(a.B) = (BP 4- S) (BP — N) (aP +N) (aP — S) 

N(a.B) = (a4- [D CgP + S) PaP +N) (4.18) 
这 里 ,M(a,B) 为 GPNSS ARAIRE RE. 要 证 明 GPNSS 迭代 格式 是 收敛 的 ,只 需 证 明 
和 迭代 和 矩阵 的 谱 半径 oCM GB) —1. 

由 于 5 为 反对 称 和 矩阵 ,P S=R (OU) S =R? GR? ) SR OR, a PS 所 有 的 
特征 值 均 为 虚数 , 设 P71S 的 第 j 个 特征 值 为 iej Cj — 1.2. 0i V —1. 下 面 给 出 
GPNSS i1 (C RS lic Stc E E. 

定理 4.3 令 AEC"** 是 一 个 正定 矩阵 ,NEC"" 是 一 个 正规 矩阵 ,SE CU" à — 
XP KEE Jf Fili i£ ASNES. PEC*" 是 一 个 埃 尔 米 特 正定 和 矩阵. O FO, SPE PN 和 
P'S 特征 值 的 集合 ,对 于 V4 二 Xi 十 iy EO1,k 二 1,2,…,n,7 表 示 特 征 值 的 实 部 ,和 表示 
特征 值 的 虚 部 ,i= V 一 1 ,引入 下 面 的 记号 : 


Yu = max {Ye}s Ya» = min (y, 


Js = max (|m|)« qms = min {| me |} 
对 于 0. 中 的 特征 值 ,引入 下 面 的 记号 : 
eux = maxi |e |}, min = min { |e; |} 
DEL DEL 


那么 ,我 们 有 pCM Ca) —oCa [D ,其 中 ， 
re 


no ee, N Ca + Y)? 十 Nk eee, 


如 果 参 数 a 和 有 满足 下 面 两 个 条 件 之 一 : 


n 
(1) qu rum GE Ua, ,其 中 

NQ = (Ca. | at B BK), Qs = (Ca. D | B< min(a.Ka)) .. Pilaf) <0}, 

Os = (Gf | BGO <B<a, WylasB)<0}, Q, = (G4 | BS max(a GO) ) . 
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函数 W Ca B) V, Ca D) 和 Ba) 定义 如 下 : 
V, Casp) = Ol, — els) Ca — B) — 208% max — Zein Ymax + Ca — fous (4. 19) 
V, (a. = (Pin — el) (a — B) — 208% min — 2els Yon + (a — Bus (4.20) 


2Y cain Vmax F AC Vmin 十 Yun) — 21] mar 
aa Za F Yass F Yain 


(4. 21) 


: 
(2) 当 tae > M VexYnin > Cap) € UL I ;其 中 
i=l 


Il, = (Ca, | a< p} = (Ca. | Bla) SPa Vi Ca.) <0}, 
II; = {(a.B) | B< min(a.fCG) V Ca. B) <0}, 
函数 P (a,B) Plab) 和 Bo) 定义 如 式 (4. 19) 3X (4. 21). 
则 我 们 有 下 面 的 不 等 式 成 立 
p(M(a.B)) ola.) — 1. 
HU GPNSS 迭代 方法 收敛 到 线性 系统 (4.1) 的 精确 解 x* . 
证 明 : 利用 相似 矩阵 的 性 质 以 及 谱 范 数 的 性 质 , 我 们 有 
(CP +S)M(BP + $)?) 
=p((BP—N)(aP+N) (aP — S) (aP - S)) 
p((BI — P> N) (aI +P>N)™ (al — P? SY(BlI + P?^S)7) 
x || (BI — P No (al + P^ N)? || * || CI — P^ S) (BI + P^ 55? |l; 


B—A . a +e 

ath| «ee |B e 

= (8— A! + 3" a^ +e 
nine (aty) +o | eo [A B. +e 


因为 a— B+ 2y>0, E REV ((8—30* + 9?) / CC yy? +9?) 关于 变量 7 递增 ,因此 ， 


— max 
à,€6, 


有 
max Ponti 一 max (B— Ye)" ses 
nee A Ca HV HNE numm rs N Ca F Yn)? F Nax 
— max E (B= Yu one (4, 22) 
(a+ Yuin)? + pine N (a+ Vmax)? TF 
类 似 的 ,我 们 有 下 列 不 等 式 
1 
z z ln r3 AL ax 
a Te _ Ë Cmax 
max |z = (4. 23) 
wea VB te le +e, 
zo a 
M Emin 
下 面 ,我 们 分 两 种 情况 讨论 . 


(mY 
利用 函数 BCa) 的 定义 , 式 (4. 22) 的 右 端 可 化 为 
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B= Yous)? + bar 
(B= Tous) 有 < fa) 
QT Ymax N max (4. 24) 
(B= V Ea > 

la F or PPPO 


下 面 我 们 把 区 域 
D = (GB |a>0.8> 0}. D-lJ D, 


分 成 四 个 部 分 (如 图 4. 1 所 示 sao = I Yenax Vmin — 2x) ,分 别 记 作 
D, = ((a, |a<B< Bla)}, D; = (Co, D | B — min{a,Pla)}} 
D; = (C2, | Bla) x; B a), D, = {(a,8) | B> maxia,Pa)}} 


pi 


图 4.1 DD 的 四 个 子 区 域 


(a) 当 (a,B) € D, 时 ,由 式 (4.22)、 式 (4.23) 和 式 (4. 24) ,我 们 有 


(B= Yeux)? F ous la” + ea 
, . 1 
pMa < JETE EX EN ETE 


(b) 当 (e,B) € D; 时 , 谱 半 径 满足 


pM p) < [E Tu 7- . ERA 
上 述 不 等 式 严 格 小 于 1 当 且 仅 当 
V, (a8) = (Pag — ela (a — B) — 2087 maz — ets Youn (a — B) Par 0 
CO 当 (a,B) € Da: 时 , 谱 半径 满足 


(8— Y.) DE. JE Cmax 
Map < fE NET Vi tee <! 


当 且 仅 当 


V, (a) = OL. — ela) (a — B) — 2afY mas — Debs Ys + Ca — [ha 0 
(d) 当 (a,B) ED, 时 ,我 们 有 


mom 


mon 。 增 量 未 知 元 及 其 预 处 理 迭 代 算法 


(8— Dat Br C Ru. 
Map < fE 1 


自然 成 立 . 
因此 , 当 mus V Ymax Yin 时 ,有 


4 
o(M(a-B))<1, Wap EU a, 
i=l 


(2). taux A nux uin 
利用 函数 BCa) 的 定义 , 式 (4. 22) 的 右 端 可 化 为 


B= Xa qus. 二 


(at You)? ah. Nea “SF 
e Em s um «2 p> pla) (4. 25) 
P= ys) E ga a en ] 
Co ga Js rue. B< minfa GO) 
fi ELIO RIEKIE D= (Caf) 12770,2770) HUS HE SHIN D= UD? i1. 


2,3. 其 中 
Di = (Glas. Di = {(a,p) | a >B BG). 
Di = (CG | B= min{a,Bla)}} 
(a) Ha, p € Dy 时 ,由 式 (4.22)、 式 (4.23) 和 式 (4.25) ,我 们 有 
(B= Yuin)” F Yow, Js + ehar 
PM (asf) < Ee Fie NP tea = 
(b) 当 (a,B) € Ds 时 , 谱 半 径 满足 
(B= Ya)? E Tinax M a? + emin 
RUMINUS. ex Yun)! + qux NE + euin 
不 等 式 严格 小 于 1 当 且 仅 当 VG.) <0. 
CO Hap) € Dj 时 , 谱 半 径 满足 


o(M(a, m < [Bata Í à 和 一 1 
当 且 仅 当 V, (.)—0. 
因此 , 当 Nmax > V Ymax Vmin MF e RATA 


3 
Map —1. VG. €) m 
i 一 1 


4.3.3 最 优化 策略 


下 面 我 们 引入 记号 
(a* .B' ) = arg mince; i 
定理 4.4 给 出 使 得 GPNSS 和 迭代 矩阵 谱 半 径 的 上 界 取得 最 小 时 参数 (at ,8" ) 的 计算 
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方法 . 
定理 4.4 FF three nx Vain * Vs (a) 70. V, Ca 770 3E AA f f GPNSS 迭代 矩阵 
谱 半 径 oCMCG B). 的 上 界 取得 最 小 时 参数 (a,B) 值 为 
(ay BC 2) + Tuis Tmax — oak ein 
(a* ,B* ) = 4 (a0 »Blao))s ein < Tuin Taux = Tas e Cree 
(az ,6(az)) ， Tu Tmas — Niar > Cex 
其 中 ， 
V, (as) =— emin Tmin — BY sax t B Tuin — BY int [V nin 
V, (a, B) == ek Tuis — BY Snax + D. Tuis — BY st [Ein 


a. em disc Christ = Tui oua Di 
i Tmin 十 Tmax 

p aT Cmax — Tuin Tmax + Ao 
: Tus Tus. 


Ai = M Cen + T uis) Chin FT max) F [HN HRS Zekin — 2 Tuis Tmax) 


Ae letar + T in) Cemax + Thaw) 十 T] max C] + Zehar — 2 Twin Tmax) 
除 此 之 外 ,利用 不 等 式 (4. 23)、(4. 220 ,我 们 有 
ola1), — Tuis Tmax — Wax qnin 
ala” ,B' ) = jc(ao)， elis Tuis Vmax — Nax ekax 
o(a2)5 — Tuin Tmax — inex Cmax 
其 中 ,cla) MRA. 27) 定 义 . 
WEBB. 利用 不 等 式 (4. 23) 和 不 等 式 (4. 24) ,有 


(B= Ymax)? F 7max ? rela 
JE ee AN Em. ， (asp) € D, 


(B= Ya Hha ,— Ja F emin » (asp) € Dr 
o(a,p) = NGF Ye e NE 十 em (4, 26) 
(B= Yuin)? + pma , fa” + enin 
(estes Fre.’ (a.B) € D; 
(B= Vota)” gue, [°F Cae 
Se Pren?’ (a,B) € D, 
通过 计算 ,我们 发 现 当 (a,8) ED UD, 时 ,函数 ol(a,B) 关 于 B 的 偏 导数 og(a,B) 二 0 利用 定 
H 4.3 的 条 件 V Ca 770 F, Ca 770 TA 4 (a) € D, UD, 时 ,函数 ola,B) 关于 有 
FY fi FBC og Ca 0 70 AE « ERE Ca p) AY Se 7] ET TK HE DC dD, UD, 和 DS U D, 的 交 线 
B— BG) EU. 为 了 找到 函数 函数 c(e,B) 的 最 小 值 ,只 须 找到 函数 


(BC@) = Yuin? F pm a + el, 
{a= Bue dote 


(Bla) 一 Yaad F Tex Pa eh 
. ; 8X m. 
SERE, Ba? +e, to" 


def 
ala) =alasPla)) = 


Sih 
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的 最 小 值 . 利用 求 导数 公式 ,可 以 得 到 
Ca(o)m(la), a> ao 
o (a) = 
teas a Xa, 
HB ci 020.0 0270 分 别 为 两 个 正 函数 ,并 且 
大 
b = (Torin Tm) + 2( Yoox Tuis — Chex — aer E — nas (Tuin 
注意 到 ,六 (ec) 和 疡 (Co) 的 均 为 一 正 一 负 . 记 这 两 个 函数 的 正 根 分 别 为 wm ,az. 由 
人 
ign = (Tuis + Tax + 2 Tras Tuin ) (Lus Tmin — Nhar — Chex) 
我 们 有 下 面 的 结果 成 立 ， 
D 如 果 Tus Tuus — haces ,那么 ola,B) 的 最 小 值 是 oa), 即 ala BC) «a >a. 
@ 如 果 ehin <T min Tmar ~ Nax Seha s ABA ola,B) 的 最 小 值 是 alao), BI o Cas ,BCao)). 
@ MWA Tus Tuus Drax ora ,那么 ola,B) 的 最 小 值 是 oCas), 即 alaz ,BCas)) saz <a. 
定理 4.5 若 mwx 宇 VYmxYan ,使 得 GPNSS 和 迭代 矩阵 谱 半径 oMa, p) 的 上 界 取得 


最 小 时 参数 (ep) 值 为 


Ton 
T; 


m) 


azp), Tint ha Chin 
(a* p" = [97 77») , ein < Tat Tix S hax 


(ai Bi ) ， Tiat Ninax > ebes 


其 中 ， 
= (Dect gla + eben + V a F em) Thin) Cian — Cran)? F Chin) 
"M 2T, 
= (Viet rin) + erin + M CC] + erm) + Thain) Cena — min? + Thin) 
a2 2T. 
2 (Tint Nax) — ehar H A Omas F Cmax)? F Thin) (Cmax — Cmax)” F Thin) 
; 2T, 
B Pint ua deb Ew COnus Iun. E ee ee Lais) 
i 2T 
相应 地 ,函数 ol(a,B) 的 最 小 值 为 


ola: fs Ys That Pua rnin 
oa” ,B' ) — 4olpospo)s elis << Taint Prax Crax 


olai RO. Th pax > ela 


其 中 po = VT int pax - 
WEBB. HX (4.250 4 
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(一 0 并 十 这， a^ des. . 
ies Ya + im NF Fos’ COP € Di 


=) [B Yun + hms 。 fo 十 em ; 
ap . e (ap) € D 
dad Glial ae Wd me Pe 


一 PR } 2 十 
(B Yous)” Jw , [8 85 | (a,8) € D; 


显然 , 当 (a,B)E D3 WE PAB o Ca DKF B BU ii FB op Ca p) <0, We ok E SCR eE SE 
价 于 求 下 列 方程 的 最 小 值 : 


(8 一 yam) F Nmax Zur" . 
CIE wu Gu € Di 


B= Ton)? og fol dens " 
(a 十 Ymn)? 十 ass EEA dio d 
为 了 求 出 最 小 值 ,我 们 首先 计算 函数 o Ca D. KIAR o Ca D «og Ca [D ,事实 上 ， 
cı Ca BD mi (a) (asf) Dé 
atap = l D pe 
€; Ca Dj; (a) (a) € D? 
其 中 ,a Cap) ,cso(a,B) 分 别 为 两 个 正 函 数 ,并 且 gi ,7: 为 关于 a 的 两 个 函数 ， 


H (a) ym 


+ CL iin j — Cmax 20 — Torin Cmax 
, 
| 


Ne (a) 人 ( Tis Tinax Cain da Tmin m 


类 似 地 ,有 
C Ca Dg. (D Ca) € Dt 


og Ca. = L a 
&G«3.(QD (D € D; 
FE C3 (a) ,cla,B) 分 别 为 两 个 正 函 数 , 并 且 了 s qp WRF 8 的 两 个 函数 ， 
TD = Torin? + Gl, — (TEint 05:008 — Trin Cree 
ja = Tuff + Gli — (Lint 21,))8— Train Cin 
注意 到 71 (a) +792 GO g (QD. 都 各 有 一 个 正 根 和 负 根 . 我 们 把 上 面 4 个 函数 的 
正 根 分 别 记 作 xi ,zs Ba ,5，, 我 们 有 下 面 的 结果 ， 
Ca) 如 果 Phin + quels RAT a0 Bi a2 Ba ,因此 ,函数 cola,P) 在 区 域 Di 内 取 
不 到 最 小 值 ,而 在 区 域 {(c,8) la >p). RZ ola PTEN Ca: s Bs) 处 可 以 取 到 最 小 值 . 事实 


Esos Bit [C Tee cde [E 在 区 域 {e,B)1o>0,p> 0 内 的 最 小 
值 ,因此 


(B= Yu 2? + qax a^ + etin 
(GB) < : . 
«af i Fyn) qu. NF da 


(B— Yan). + Nax e+e, 
EA ES DES: [EP ax 


最 后 ,我 们 要 说 明 的 是 ,由 于 Bs 之 BC(az) 点 (aa Bi) EX BR DI VI DAC H Tiia tya 
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elis IST s PR Co B) YE DX BR Cap) 1770 8770) PEEL MEE olab). 

(b) A elin — T2 + ax Keun IBA as Bi a2 Bs ,容易 看 出 点 (zi,8,) 不 在 D? . 
i=1,2 内 ,因此 ,最 小 值 只 能 在 直线 =8 上 取得 最 小 值 , 即 最 小 值 是 使 得 cCa,a) 导 数 为 零 
的 点 , 记 为 po, 最 后 ,可 以 计算 出 po = VT hin F Mear ,函数 在 {(a,B) |a>0,8>0} 内 的 最 小 
值 为 ouo ,po). 


Co) 如 果 T2, queas TEE @ 1 <B 1 ass ,因此 ,函数 cla,B) 在 区 域 D; 内 取 
不 到 最 小 值 , 而 在 区 域 Di ,函数 ola,B) 在 点 (zi ,Bi) 处 可 以 取 到 最 小 值 . 事实 上 ,elai， 


Gaia tg. Te ín 
TT. ER. NL. 在 区 域 {(a,B) [a0 8770) PS AR MET EL 


- 万 (B= Yuin)? + nux a^ + en 
Gi B) < : . 
MULA E ES Jas N F Ee 


(8— Y yi 1 oe? +e, 
x LE. E 
= ETRE Wea’ «ce 


IE, H Tein HN has > esa FY, PR oC, 8) 在 在 区 域 {(a,B) 1a770, 8770). 内 的 最 小 值 是 
o(a Bi). 


4.3.4 不 精确 两 参数 预 处 理 NSS 和 迭代 策略 


为 了 提高 两 参数 预 处 理 NSS 迭代 方法 的 效率 ,我 们 可 以 利用 迭代 方法 对 两 个 子 线性 
代数 方程 组 近似 求解 ,于 是 得 到 下 面 的 不 精确 两 参数 预 处 理 NSS 迭代 方法 . 

不 精确 的 两 参数 预 处 理 NSSCGPNSS) 和 迭代 方法 : 给 定 一 个 初 值 xY CC" ,对 于 人 一 0， 
1.2. HA Ce rS, A x ? 为 初 值 , 利 用 某 种 迭代 方法 (例如 CG 方法 ) 近 似 计算 


GaP--NDx U 3) x aP SEO Ho MA ER r CHH RADA x CH X LEERLO 
种 迭代 方法 (比如 Krylov FA li] Jr JE BUT PHS x ~ (BP — ND x (3) +b 的 
fit 4B) xt 

为 了 简化 数值 计算 和 收敛 性 分 析 , 将 GPNSS 迭代 格式 写成 如 下 等 价 形式 ; 

给 定 一 个 初 值 EO EC', 对 于 k= 二 0,1,2,…, 直 到 {5 中 收敛 . 

QD EWR TR ZR TE); BUB aP 2- N)z zer (F^ —b — Ax ^ ) isst p? =F 一 
(aP+N)z ? ,满足 

Ip? I «e llz |, 

然后 计算 ECE) 一 元 LE®, 

© 近似 求解 线性 方程 组 (BP 十 9)z 2) p UD (UD Lp Ag CD ) 直 到 残 
at 2) =F CO?) — (ppt sy?) ee gD pen | pO Lo 
ee =y (+) + ECT) RM, || . || 为 向 量 范 数 . 

FEMA 


A =M; — N; = (aP +N) — (aP — S) 
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M; — N: = (BP +S) — (BP — N). 
引入 下 面 的 范 数 ， 
ll x Illa = I| Mx ll. 对 于 所 有 的 xecC* 
其 诱导 范 数 为 
ll X lle = || MXM™ ll. FRAK X E CP” 

这 时 , 若 令 M==BP 十 $, 我 们 有 下 面 的 收敛 性 定理 : 

定理 4.6 令 AEC"* 是 一 个 正定 矩阵 ,NEC”" 是 一 个 正规 矩阵 ,SE Cr BFR 
对 称 和 矩阵 ,并 且 满 足 4 王 N 十 S. PEC 是 一 个 埃 尔 米 特 正定 矩阵 x” EC" 为 线性 方程 
28 (4. D AYA. (x 9) 为 GPNSS 方法 得 到 的 迭代 序列 , 则 有 下 式 成 立 ， 

||| x**? — x* [l| Colas D + dog 1 + &,) ||| x? —x" Hl» 有 一 0,1,2，， 
其 中 [ll x lll ARE I xe llo I| (BPS) aP +HN)™ || ,,0= || AP +HS)™ || 2. 

特别 地 ,如 果 (oCa,B) 十 Gomwax) (1 十 emox) 二 1, 那么 迭代 序列 {站} 收敛 到 线性 方程 组 
DME x EC" ,这 里 Emax = max (e, } Has = max{ Me E 


4.3.5 数值 算 例 


例 4.2 求解 下 面 的 反应 扩散 方程 
— Qu, us) bu, = fs (rey € 
luco, (x,y) € AN 
JEP 020 ESAERA (0.107. 使 用 网 格 步 长 为 为 h 的 五 点 离散 格式 


b : 
(Aug — ui — Wy — Uia — uia) d gg Ctm uia? fü: 


h* 
这 里 wj A fi; 分 别 表示 w 和 了 的 近似 值 ,得 到 下 面 的 线性 方程 组 
AU =f 
和 和 矩阵 的 HS 分 裂 不 同 ,对 于 给 定 的 矩阵 A EE NS 分 裂 不 是 唯一 的 ,这 里 选取 
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Kp i= yoI, ESAR. 

在 计算 中 ,选用 IU 预 处 理 和 ILU 分 解 预 处理 , 用 4 种 迭代 方法 计算 ,这 4 种 迭代 方 
法 分 别 是 NSS 迭代 、GNSS 和 迭代.GPNSSGIU) 和 GPNSSCOLUO.. 比较 迭代 和 矩阵 的 谱 半径 
和 收敛 速度 ,结果 如 下 . 


1. 谱 半 径 和 收敛 域 

我 们 首先 考察 GPNSS 方法 的 谱 半径 和 收敛 域 . 为 了 有 效 地 证 明理 论 的 正确 性 ,数值 
实验 分 两 步 进 行 : 

@ 取 c=1,6 二 1, 可 以 验证 ,在 NSS 迭代 、GNSS 迭代 和 GPNSS(IU) 三 种 方法 中 


Tea < V Ymar Vain * Toin Tax — Nax > Coax ;而 在 GPNSS(ILU) 方 法 中 ,por > V Yoan Vain Coin Coax 一 
Th 之 ex ;利用 定理 4.3, 我 们 可 以 算出 GPNSS(IU) 的 收敛 域 ( 见 图 4. 2) 
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4 
C= {(a, | a> 0.827 Bi GO ) =U 0. 


Bh 
pra 


B=B(@) 


B-B.(o) 


图 4.2 GPNSS(IU) 方 法 的 收敛 域 
Job m DE 33. 根据 定理 4. 5 可 以 计算 出 参数 (e,B) 的 值 为 (ws ,B(wz 0. iif 
半径 的 比较 结果 见 表 4.5,mo 代 表 IU 离散 的 粗 网 格 数 . 数据 表明 ,GPNSS(CIU) 和 迭代 和 矩阵 
的 谱 半径 要 远 小 于 其 他 三 种 迭代 方法 . 
表 4.5 0 一 1 时 ,NSS 方法 和 GNSS 方法 收敛 半径 的 比较 


No 3 4 5 6 

a 1. 7321 1. 1589 0. 7265 0. 2680 
NSS 

e(M(a)) 0. 8045 0. 8620 0.9077 0. 9596 

a 0. 6180 0. 8172 0.9421 1. 0213 
GNSS B 1.1849 0. 9574 0. 8692 0. 8278 

EOM Ga) 0. 7687 0. 8520 0. 9058 0. 9336 

a 0. 9042 1.0168 0.9306 1. 2488 
GPNSSCILU) B 1. 2272 1.0963 1.1969 0. 8905 

pOM Ca « 0) 0. 7609 0. 8283 0. 8727 0.9077 

a 0.5911 0. 3595 0. 4799 0. 6063 
GPNSS(CIU) B 2.0320 1.5571 1. 2626 1.0700 

pOlM Ca.) 0.5678 0. 6562 0. 7275 0. 7855 


Q 取 c=3,6= 二 10, 可 以 验证 ,在 四 种 方法 中 均 有 加 ss 之 V Ymax Vmin * Ehin < T min + T] max S 
eh » All FH aE 4. 3 ,我 们 可 以 算出 GPNSS(IU) 的 收敛 域 ( 见 图 4. 3) 


3 
C — (G4 | a> 0,877 GO) = UJ mr 
i=1 
8. 13a— 1. 10 


B: — 
Fb G)— 999078. 13 


注意 ,集合 DL BER. 根据 定理 4.5, 可 以 计算 出 参数 (a,pB) 
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图 4.3 GPNSS(IU) 方 法 的 收敛 域 


的 值 为 (yo ,wm ) , 谱 半 径 的 比较 结果 见 表 4. 6. 注意 到 ,此 种 情况 下 ,GNSS 方法 即 为 NSS 
方法 . 


表 4.6 b—10 时 ,NSS 方法 和 GNSS 方法 谱 半 径 的 比较 


No 3 4 5 6 

a 3.0473 3.0154 3. 0064 3. 0031 
NSSCGNSS) 

p(M(a)) 0. 810 0. 8972 0.9292 0. 9498 

a 3. 0697 3.0352 3.0178 3. 0094 
GPNSSCILU) B 3. 0697 3.0352 3.0178 3. 0094 

pOMCG «0 ) 0.8877 0. 9101 0.9271 0. 9401 

a 3. 2280 3. 1370 3.0740 3. 0403 
GPNSS(CIU) B 3. 2280 3. 1370 3.0740 3. 0403 

pOM Ca) 0. 7934 0. 7927 0.9141 0. 8726 


2. 收敛 速度 的 比较 


本 节 主 要 比较 4 种 迭代 方法 的 收敛 速度 , 表 4.7 和 表 4. 8 分 别 给 出 在 O=1 和 0 一 10 
时 的 迭代 步 数 (IT) 和 CPU 时 间 的 比较 ,在 精度 为 | ~ || :<10 一 时 ,无 论 IT 上 还 是 CPU 
上 ,GPNSS(IU) 都 具有 一 定 的 优越 性 . 
表 4.7 b=1 时 ,CPU 时 间 和 和 迭代 步 数 (IT) 的 比较 


No 3 4 5 6 


IT 41 62 98 264 
NSS 


CPU 0. 2968 1. 0938 3. 4063 20. 9060 
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续 表 
No 3 4 5 6 
Er 37 62 94 134 
GNSS 
CPU 0. 2968 1. 0313 2. 9688 7. 4375 
n 35 50 69 95 
GPNSS(ILU) 
CPU 0. 0937 0. 2343 1.0625 2.5469 
IT 17 23 30 40 
GPNSS(CIU) 
CPU 0. 0937 0. 3281 0. 9843 2. 2344 
表 4.8 b—10 时 ,CPU 时 间 和 和 迭代 步 数 (IT) 的 比较 
no 3 4 5 6 
IT 61 94 139 195 
NSS 
CPU 0. 3130 1.1875 3. 7344 8. 4063 
IT TI 99 126 160 
GPNSS(CILU) 
CPU 0.1875 0.5312 1. 5000 4. 2188 
IT 25 40 60 71 
GPNSS(IU) 
CPU 0. 2185 0. 4843 1.5156 2.7500 
例 4.3 考虑 线性 代数 方程 组 Ax —b. BOB IE 
[| 
A= 
一 BT 0.51 
其 中 
pa 9781 0 le — 
0 IGT-TGI 
一 E. E Ri 
FGI 
T—tridiag( —1.2, —1) € R"*", F=ôh + tridiag( —1,1,0) € Re" n= 为 网 格 步 长 . 
由 于 HSS 迭代 为 NSS 迭代 的 一 种 特殊 情况 ,所 以 ,在 数值 实验 中 ,我 们 选择 
N= laca j s= a-a’) 


通过 计算 ,我 们 发 现 
Wo < a/ i s ihm < Ta Ti < max) 
由 定理 4. 3,GPNSS 迭代 方法 无 条 件 收敛 ,上 且 两 个 参数 a p 取 值 相同 . 事实 上 ,a=B= 


Tuis Cmax , 表 4.9 给 出 了 和 迭代 矩阵 谱 半径 的 比较 , 表 4. 10 给 出 了 和 迭代 步 数 和 CPU 时 间 
的 比较 . 
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表 4.9 8 —10 时. 谱 半径 的 比较 


m 5 10 15 20 
a 1.9319 1.1268 0. 7806 0.5963 
NSSCGNSS) 
pM GO) 0. 7099 0.6148 0. 8335 0. 9086 
a 1.8157 1.0646 0. 8519 0. 7763 
GPNSSCILU) B 1.8157 1.0646 0. 8519 0.7763 
EOM Ca B)» 0. 3646 0. 5287 0. 6034 0. 6318 
表 4.10 6 一 10 时 ,和 迭代 步 数 和 CPU 时 间 的 比较 
m 5 10 15 20 
IT 34 60 93 123 
NSS(GNSS) 
CPU 0.3405 7.1903 78. 6130 428. 5062 
IT 0.6180 0. 8172 0. 9421 1.0213 
GPNSSCILU) 
CPU 0. 1102 2.9242 44. 0634 196.5727 


注 4.5 对 于 给 定 的 矩阵 4, 它 的 正规 /反对 称 和 矩阵 分 裂 不 是 唯一 的 ,为 了 方便 ,在 数 
值 试 验 中 ,我 们 选择 了 对 称 / 反 对 称 分 裂 . 为 了 获得 更 快 的 收敛 速率 ,当然 可 以 选择 其 他 
类 型 的 正规 /反对 称 和 矩阵 分 裂 , 这 也 是 我 们 未 来 要 研究 的 问题 之 一 . 
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基于 SOR 和 迭代 的 复 对 称 线性 系统 的 
MHSS 加 速 方法 


5.1 简介 


复线 性 代数 系统 广泛 出 现在 科学 与 工程 计算 的 许多 领域 ,例如 光学 层 析 成 像 , 量 子 力 
学 ,分 子 散射 和 结构 力学 等 .其 中 复 对 称 线性 系统 ( 即 系数 矩阵 的 实 部 和 虚 部 均 为 对 称 矩 
阵 ) 由 于 其 特殊 的 结构 和 广泛 的 应 用 背景 ,吸引 着 越 来 越 多 的 学 者 的 关注 . 针对 复线 性 方 
程 组 的 直接 解法 是 最 早出 现 的 求解 方法 ,但 是 当 矩 阵 规模 较 大 时 ,由 于 直接 解法 的 存储 量 
和 计算 量 一 般 很 大 ,所 以 迭代 法 在 近年 来 越 来 越 受 到 重视 . 目前 ,迭代 法 已 成 为 求解 大 型 
稀疏 复 对 称 线 性 方程 组 的 主要 方法 . 然而 ,由 于 复 系数 和 矩阵 在 方程 组 的 规模 较 大 或 系数 
矩阵 的 条 件数 较 大 时 ,系数 矩阵 容易 呈现 病态 特征 ,所 以 许多 适用 于 实 线性 代数 方程 组 的 
TEAR ATH «üt A Ee BE 7r 1 CCG), HHE d be BE Jr i (PCG) ,对 于 复线 性 代数 方程 组 的 求解 
效果 并 不 理想 ,存在 着 不 收敛 或 收敛 速度 慢 的 问题 . 于 是 ,很 多 学 者 对 原 线性 代数 方程 采 
用 适当 的 预 处 理 技术 ,例如 SSOR 预 处 理 、 多 项 式 预 处 理 \ 不 完全 分 解 预 处 理 .C-to-R 等 
预 处 理 方法 ,用 来 降低 系数 和 矩阵 的 条 件数 ,改善 矩阵 病态 特性 ,以 加 快 迭代 方法 的 收敛 
速度 . 

对 于 复线 性 代数 方程 组 的 求解 , 另 一 种 方法 就 是 将 其 转化 为 等 价 的 实 系数 线性 代数 
方程 组 ,这 样 就 可 以 将 实 系数 线性 代数 方程 组 的 解法 用 在 复 系数 线性 代数 方程 组 的 求 
解 上 . 

例如 ,我 们 考虑 下 面 的 复 系数 线性 代数 方程 组 的 求解 问题 ， 

Ax =b (5.1) 
Hp ASWHiT.x=ytiz.b=dt+ig.iX H EE W.T HS OBI. yz. deg 都 是 实 


向 量 ,i 一 V 一 1 为 虚数 单位 . 方程 组 (5.1) 可 以 写成 下 列 形式 : 


[e wle Ls] B 


这 里 需要 强调 的 是 ,方程 组 (5. 1) 的 实 系数 形式 并 不 唯一 ,例如 ,也 可 以 将 其 写 为 


MEM BALIH - 


然而 普通 迭代 法 求解 式 (5. 2) 或 者 式 (5. 3) ,其 收敛 速度 往往 要 慢 于 最 初始 的 复 系 数 线性 
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代数 方程 组 (5. 1). 尽管 如 此 ,我 们 依然 可 以 借助 复 系数 线性 方程 组 的 实 系数 形式 ,构造 
合适 的 预 处 理子 ,来 加 快 迭代 方法 收敛 的 速度 . 

近来 ,基于 矩阵 的 对 称 和 反对 称 分 裂 (HSS) 技 术 , 白 中 治 ,Golub,Ng 等 提出 了 修正 
的 对 称 和 反对 称 分 裂 (MHSS) 方 法 ,MHSS 迭代 方法 保持 了 HSS 迭代 方法 的 优点 , 即 对 
于 格式 中 Ya>0,MHSS 迭代 均 无 条 件 收敛 到 复线 性 代数 方程 组 (5. 1) 的 解 . 同时 又 克 
服 了 HSS 方 法 中 的 不 足 , 即 避免 了 在 迭代 中 计算 系数 和 矩阵 为 复数 矩阵 的 情形 . MHSS 方 
TAFE £O B E W «T 为 对 称 和 矩阵 情形 下 对 HSS 方法 的 推广 。 很 多 学 者 对 MHSS 方法 进 
行 了 推广 . 王 坤 针对 复 系数 线性 方程 组 的 特点 ,对 MHSS 算法 进行 了 修正 ,提出 了 
MLHSS 和 迭代 方法 . 2013 年 , 郭 晓 霞 、 王 坤 讨论 了 矩阵 W,T 为 非 对 称 情形 的 情况 ,并 给 出 
了 收敛 定理 ,进行 了 数值 实验 ,证 明了 方法 的 有 效 性 . 

本 章 首先 回顾 修正 的 HSS 迭代 方法 及 其 收敛 定 理 , 为 了 提高 复线 性 代数 方程 组 的 求 
解 效率 ,基于 NSS 超 松 弛 加 速 迭代 格式 ,提出 一 类 求解 复线 性 代数 方程 组 的 超 松 弛 迭代 
MHSS 加 速 迭代 格式 . 然后 推导 该 加 速 迭代 方法 的 收敛 条 件 以 及 加 速 迭代 格式 中 参数 
co 的 选取 办 法 ,最 后 用 数值 实验 证 明 方法 的 正确 性 和 有 效 性 . 


5.2 MHSS 迭代 方法 


A5 MGR f AH RE A TE; RAL 
Ax —b x.be€c' (5.4) 
其 中 4EC” "为 大 型 稀疏 复 对 称 和 矩阵 , 即 4 具有 下 述 形式 : 
A—W-riT (5.5) 
AE W T 4029 Sz YE E AE [ERR Sc IE 2E zE AB Pe Fi WT —W T7 — T. i= /—1 为 虚数 单 
位 . 为 了 保证 A 为 非 Hermit 矩阵 ,我 们 假定 TAO. 为 了 避免 在 迭代 中 计算 系数 和 矩阵 为 
复数 矩阵 , 白 中 治 ,Golub,Ng 提出 了 修正 的 对 称 和 反对 称 分 裂 (MHSS) 方 法 ,迭代 格式 
定义 如 下 : 
MHSS 迭代 方法 : 给 定 一 个 初 值 x € C" MET 80.1.22 ERI x | Se A 


Girma O =a in +, 
(af + Tx = (af +iW)x C3) — ib, 
其 中 a 二 0,T 为 与 矩阵 A 同 阶 的 单位 矩阵 . 
对 于 MHSS 迭代 方法 的 收敛 性 ,我 们 有 下 面 两 个 收敛 定理 . 
定理 5.1 假设 A 二 W 二 iTEC” ,WER””, TER” ,分 别 为 实 正定 矩阵 和 实 正 半 定 
AVE ELL WT WT" — T.i— / — 18 E URL. EP Vac 0. MHSS 迭代 矩阵 的 谱 半 
径 pMa) KRI E 


(5.6) 


p(M(a)) < ola) 
其 中 


Ma) = (al 4- T) ? (aI -- iW) (aI +W) (al — iD. 
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EM Jæ +A} 
ola) = max 一 一 一 一 一 


yew a +A; 
FEA ACW) de as FE WASTE AL. 
pW(a)) Sola) <1, Va>0, 

即 对 a7-0. MHSS 和 迭代 方法 收敛 到 线性 方程 组 (5. 4) 的 唯一 解 x* € c7. 

定理 5.1 说 明了 MHSS 方 法 的 收敛 速率 只 与 对 称 正定 矩阵 W 的 特征 值 有 关 , 而 与 
和 矩阵 TT 的 特征 值 以 及 矩阵 W,T 的 特征 向 量 无 关 . 若 已 知 矩 阵风 的 最 大 和 最 小 特征 值 , 便 
可 以 得 到 使 得 谱 半径 oMa) ER ol(a) 取 得 极 小 时 参数 的 计算 方法 . 

定理 5.2 (Ri A—W--iT€ C"".W€ R"".T€ R”" 分 别 为 实 正定 矩阵 和 实 正 半 定 
和 矩阵 且 满足 WT —W TT —T.i— V 一 1 为 虚数 单位 . 令 Tw 和 了 ,代表 对 称 正定 矩阵 W 的 
特征 值 的 最 大 值 和 最 小 值 ,那么 


24 32 
a” = arg min max a |= AJ Vestn T cane 


Tuis SAST max ata 
相应 地 
a Lat NE KW) +1 
ola") 
Tmin 十 V Tmax K(W) +1 


Kp aeo = Tonos / Trin BEAD RR PE W 的 谱 条 件数 . 


5.3 基于 SOR 迭代 的 MHSS 加 速 方法 


5.3.1 MHSS 加 速 方法 的 建立 


下 面 我 们 基于 NSS 超 松 弛 加 速 迭 代 格 式 ,提出 一 类 求解 复线 性 代数 方程 组 的 超 松弛 
MHSS 加 速 迭代 格式 . 为 了 给 出 MHSS 加 速 迭代 方法 ,我 们 首先 考虑 下 面 两 个 不 动 点 
方程 

(al +W)x = (al — iT) y +b, 
(al +T)y = (al + iW)x—ib, (5.7) 

容易 得 到 ,这 两 个 不 动 点 方程 的 解 与 方程 组 (5. 4) 的 解 有 如 下 关系 . 

定理 5.3 WIR x” EC" 是 (5.4) 的 解 ,那么 同时 也 是 式 (5.7) 中 两 个 方程 的 不 动 点 
反 过 来 ,如 果 x" EC" 是 (5.7) 中 任 一 个 方程 的 不 动 点 ,那么 也 是 式 (5.4) 的 解 . 

WEBB; 设 x* € C" 是 (5.4) 的 解 ,由 方程 Ax 二 b 和 一 iAx 二 一 ip 同 解 ,以 及 

Ax' =(W+iT)x’ 

GI 十 WDOr — (al —iT)x" = b 
iAx* =(—iW+T)x* 

(al +T)x" — (al +iW)x* ib 


可 知 ,定理 结论 成 立 . 
下 面 , 我 们 将 两 个 不 动 点 方程 等 价 转换 写成 具有 2X2 分 块 结构 的 线性 方程 组 ， 
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x af 十 —(al —iD x b 
4u(a)| |= : =|. (5.8) 
y —(al--iW) a+T fly — ib 


定理 5.4 证 明了 A GO KIERRE, EM 5.5 证 明了 求解 式 (5. 1) 和 求解 式 (5. 8) 是 等 
价 的 .从 而 我 们 可 以 从 式 (5. 8) 出 发 建立 MHSS 迭代 加 速 方法 . 

定理 5.4 (RI A—W--iT€ C"".WC€R'"".T€ R""4H9lJ Sc iE xE £g [4e Rl Sz 1E ^E E 
AE HW ig W —W.T! —T.i— / —128 He GR. XT VY a>, ABI 

did | al +W B "— 
— (al + iW) of +T 

为 非 奇异 矩阵 . 

证 明 由 于 WER”".TER "分 别 为 对 称 正 定 和 对 称 半 正 定 和 矩阵, 因此 ,对 于 Yao>0， 
al +T fll al HW 为 正定 的 非 奇异 矩阵 , 除 此 之 外 ,成 立 


al +W 一 (af — iT) 
Ao (a) = . 
— Cal + iW) al +T 
I 0 al --W | — (al —iT) 
一 (5.9) 
—(al+iW) (al --W)'!I 0 Ia) 


HPP (a) = (GI HT) — GI +iW) (al HW) ^! Xa — iT) HEER Schur 补 且 
Ia) = (al -- TU — Mta) 
Ma) H MHSS a (Jr (f 35 FOR P « rh CR 113 JP EM 2.1 RE Va 0. 90M G0 — 
1. 因此 ,4。(a) 为 非 奇 异 和 矩阵 ,定理 证 毕 . 
利用 定理 5. 4 ,两 个 不 动 点 方程 (5.7) 构 成 了 一 个 2X2 块 状 线性 代数 方程 组 ,这 个 方 
程 组 与 (5.4) 等 价 ,于 是 得 到 下 面 定理 : 
定理 5.5 假设 A 满足 定理 5.4 的 条 件 . 如 果 x” EC" 是 (5.4) 的 解 ,那么 z* = 
上 上 ce 加 性 广 本 组 (5. o tuia y -l. |]ec* 为 线性 方程 组 (5. 8) 的 解 ， 
那么 x =y 必然 成 立 , 且 x" 为 式 (5.4) 的 解 . 
WEBB. 利用 定理 5. 4 和 简单 计算 可 以 得 到 定理 的 结论 ,证 明 过 程 略 . 
下 面 ,我 们 由 式 (5.8) 出 发 ,给 出 分 块 Jacobi 迭代 和 超 松弛 (SOR) 和 迭代 算法 : 
KR Jacobi 迭代 方法 : 给 定 一 个 初 值 x” ,yo EC" 对 于 k= 二 0,1,2…， 
xP = (dW [Cal — iT y'? 4b]. 
{ae = (al + T) [Cal + iW) x — ib], 
其 中 o> 0.1 558A 同 阶 的 单位 矩阵 . 
Jacobi 迭代 方法 也 可 以 等 价 地 写 为 
ZAP = J(a)z® +C(a) 


其 中 
0 (al 23-W) ^! (al — iT) 
JG) = » . (5. 10) 
— (al + T)? Cal + iW) 0 
(al +W)'b 
Cla) = . 
— Cal + T)! ib 


63. 


增 量 未 知 元 及 其 预 处 理 选 代 算法 


2X2 分 块 结构 的 线性 方程 组 (5. 8) 的 块 SOR 迭代 格式 定义 如 下 : 
块 SOR 和 迭代 方法 ”给 定 一 个 初 值 xo ,yw EC" ,对 于 一 0,1,2…， 


Qon 一 (1 一 o)x + o(aI +W) [CI —iT) y'? +b], (6.11) 
yr (—o)y® + (eI 4- T) ! [ CaI +iW)x — ib], 
其 中 a0. o 为 松弛 参数 . 工 为 与 矩阵 4 同 阶 的 单位 矩阵 . 
式 (5.11) 也 可 以 等 价 地 写 为 
z= = Tyla F (a) 

其 中 

Ta) = | Amo w(al +W) gag (5.12 

wl — o) Cal + T)! CaI + iW) (1— 9I +w MCa) 


这 里 ,M(a) 为 MHSS 和 迭代 方法 的 迭代 和 矩阵. 
注 5.1 当 w=1 时, 块 SOR 退化 为 块 2X2 线性 方程 组 (5. 8) 的 块 G-S 迭代 ,相应 地 
T, GO T 3B [Eg G-S EREE Ti Ca). 


5.3.2 MHSS 加 速 方法 收敛 性 的 证 明 


由 块 Jacobi 矩阵 J Ca) , 块 G-S 迭代 矩阵 N (a) Al MHSS 迭代 和 矩 , 容 易 得 到 下 面 的 定理 ， 

定理 5.6 O 如 果 y 是 J(a) 的 特征 值 ,那么 jy? JE T. CO WEA. WR ON Ti CO 
的 一 个 非 零 特征 值 且 pp? = 二 0, 那 么 y 为 J(a) 的 特征 值 ， 因此 块 Jacobi 迭代 方法 收敛 的 充 
分 必要 条 件 是 块 G-S 迭代 方法 收敛 , 且 如 果 二 者 均 收 和 敛 , 那 么 pT Ca) =p a). 

© Jt G-S 和 迭代 矩阵 Ty (a) 的 非 零 特征 值 集合 等 于 MHSS 3E (OB IE M GO 的 特征 值 集 
fr. 因此 ,G-S 迭代 收敛 当 且 仅 当 MHSS 迭代 收敛 ;如 果 二 者 同时 收敛 ,我 们 有 p(T, G0 = 
pu MG) <1. 

利用 定理 5.6 及 相关 文献 ,我 们 马上 得 出 下 面 的 推论 : 

推论 5.1 假设 A 满足 定理 5.4 的 条 件 . 块 G-S 迭代 方法 和 MHSS 迭代 方法 均 无 条 
件 收敛, 且 p(T C0) — (MG) =p? C0) —1. 

注 5.2 推论 5.1 说 明了 G-S 3E [CRI MHSS 和 迭代 的 渐 近 收敛 速率 是 一 样 的 , 且 二 者 
的 收敛 速率 是 块 Jacobi 和 迭代 收敛 速率 的 2 fir. 

下 面 我 们 进一步 讨论 块 SOR 方法 的 收敛 条 件 . 由 于 系数 矩阵 A (c) 具 有 块 状 性 质 A. 
我 们 可 以 得 出 SOR 迭代 和 矩阵 T, Ca) Al Jacobi 矩阵 J(a) 的 特征 值 之 间 有 具有 下 面 的 关系 : 

定理 5.7 假设 A 满足 定理 5. 4 的 条 件 . J CO ,7T,(a) 分 别 由 式 (5.10) 和 式 (5.12) 所 


定义 , 则 : 
@ 设 4 为 T。(a) 的 非 零 特征 值 , 则 由 关系 式 
A+o- 1)? = wp? (5.13) 
所 确定 的 4 为 J(a) 的 特征 值 . 


© Hv 为 J(a) 的 特征 值 , 则 由 关系 式 (5. 13) 所 确定 的 1 为 Y。(a) 的 特征 值 . 

WEBB mk. 

最 后 ,由 推论 5.1 和 定理 5. 7. (8 BR SOR 迭代 方法 的 收敛 定理 : 

定理 5.8 假设 A 满足 定理 5.4 的 条 件 . J(a) ,7,(a) 分 别 由 式 (5.10) 和 式 (5.12) 所 
定义 , 则 ， 
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(D 24 J Co) CS BRIE (E. 79 SEB oC T, (a) <1 的 充分 必要 条 件 是 0 —2. 
© 当 jJ(a) 的 某 些 特征 值 为 虚数 时 ,如 果 对 于 某 个 整数 rE (0,1) ,7,(a) 的 每 个 特征 
fi w= 王 9 二 ip, 点 (0,8) 都 落 在 椭圆 


«p (6.2): e+e = 1} 


KARE Oo 3 pT.) 1. 
WEBB OM. 
5.3.3 数值 实验 


下 面 ,我 们 用 数值 实验 来 检验 SOR 方法 对 MHSS 的 加 速效 果 . 
考虑 求解 下 列 复线 性 代数 方程 组 


IC | 38, (x | =n =b 
HP r RREH K KER” = I@V,, +V,, O1.V,, =h tridiag(— 1,2, —1) € R”*”" X 
在 [0,1]X[o,1] 上 满足 齐 次 边 值 条 件 的 拉 普 拉 斯 算 子 经 过 五 点 差分 离散 所 得 到 的 矩阵， 
因此 ,K H X ERMEE nm? ,空间 步 长 4 一 二 二, 右 端 向 量 5 的 元 素 


a—ij 
roD? 


XE i= VY 一 1 为 虚数 单位 . 经 过 计算 ,得 到 
W —K4 SEA, T=K4 $8, 


b = = 


为 了 方便 ,在 SOR 加 速算 法 3 中 ,参数 a 与 MHSS 方法 中 参数 的 选取 一 致 ,另外 ,在 
所 有 情形 下 , 取 o—1.2. 从 表 5. 1 中 可 以 看 出 ,在 各 种 测试 下 ,基于 SOR 迭代 的 MHSS 
加 速 方法 在 迭代 步 数 (IT) 和 时 间 (CPU) 上 都 要 优 于 MHSS 迭代 方法 . 


R51 和 迭代 步 数 和 迭代 时 间 的 比较 


— MHSS 方法 MHSS 加 速 方法 
m m 
IT CPU IT CPU 

20? x 20 iii 24. 273 71 9. 921 
25? X25? 134 82. 540 87 37. 752 
30? x 30? 159 176.51 103 114. 86 
40? X40? 207 966. 72 135 645. 12 

5.3.4 总 结 


本 章 提 出 了 一 类 求解 大 型 稀 复 玻 线 性 代数 方程 组 加 速 MHSS 方法 ,并 分 析 了 方法 的 
收敛 性 质 ,数值 实验 证 明了 加 速 MHSS 方法 在 迭代 步 数 和 和 迭代 时 间 上 都 优 于 MHSS 方 
法 . 最 后 ,我 们 要 强调 的 是 ,为 了 方便 ,本 文中 加 速 方法 中 松弛 参数 o 的 值 选 为 1. 2 ,事实 
上 ,w 最 优 值 的 讨论 是 一 个 很 有 意义 和 难度 的 问题 ,将 在 以 后 继续 研究 . 


第 6 
非 线性 系统 的 迭代; 


6.1 简介 


本 章 我 们 主要 研究 形 如 

F(x) =0 (6.1) 

fS CHER ii AE Be VE AR BOT FE A RUR RIA. XE, F: DCC >C ERAK D 上 的 一 个 非 
线性 连续 可 微 函数 . F(x) 的 Jacobi 矩阵 F GO FE fiib AE PIE ERI. 对 于 方程 (6.1) 的 
数值 求解 ,已 经 有 很 多 学 者 做 了 大 量 广泛 深入 而 富有 成 效 的 研究 工作 ,并 且 已 经 得 到 了 一 
系列 实用 高 效 的 计算 方法 . 在 迄今 发 展 起 来 的 各 种 求解 非 线性 方程 组 (6. 1) 的 数值 方法 
中 ,Newton 法 是 最 基本 和 最 重要 的 方法 之 一 .Newton 法 的 优点 是 收敛 速度 快 ,其 不 足 
则 是 好 的 初 值 难以 选 到 ,并 且 在 迭代 的 每 步 均 需 计算 和 形成 Jacobi 矩阵 ,以 及 需要 精确 
求解 Newton FFE. 当 Jacobi 矩阵 难以 形成 ,或 者 当 问 题 的 规模 较 大 时 ,Newton 法 的 计 
算 代 价 就 很 昂贵 . 此 外 ,不 具有 全 局 收敛 性 质 则 是 Newton 法 的 另 一 个 缺点 . 为 降低 
Newton 法 的 使 用 成 本 ,在 1982 年 ,Dembo 和 Eisenstat 提出 了 求解 非 线 性 方程 组 的 不 精 
确 Newton 法 . 不 精确 Newton 法 就 是 在 Newton 法 的 每 步 迭 代 只 对 Newton 方程 进行 近 
似 求解 . 不 精确 Newton 法 的 本 质 就 是 一 类 内 外 迭代 算法 . 外 迭代 为 经 典 Newton 法 ,内 
迭代 则 可 采用 任何 能 够 准确 而 有 效 地 计算 Newton 方程 解 的 线性 迭代 法 . 由 此 ,许多 求 
解 线性 方程 组 的 迭代 方法 ,例如 GMRES 方法 .CG 方法 ,SOR 方法 等 ,都 可 以 用 来 求解 
Newton 方程 ,从 而 得 到 相应 的 Newton-GMRES 77 ik , Newton-CG 方法 、Newton-SOR 
方法 等 理论 分 析 表 明 ,不 精确 Newton IX FE 2K ft A JUBERR LAE ZR TEZ; Fe £H BY BOTT. 
2010 年 , 白 中 治 、 郭 学 萍 等 将 求解 大 型 系数 非 埃 尔 米 特 正定 线性 方程 组 的 对 称 / 反 对 
称 分 裂 HSS 和 迭代 方法 与 不 精确 Newton 法 相 结合 ,得 到 了 一 类 Newton-HSS(NHSS) 方 
法 ,数值 实验 证 明了 NHSS 的 稳健 性 和 有 效 性 . 然而 NHSS 方法 只 具有 局 部 收敛 性 质 . 
2011 年 , 郭 学 萍 等 对 NHSS 方法 进行 了 改进 , 即 在 不 精确 Newton 方向 上 配 以 线性 搜索 
后 退 策略 ,提出 了 一 种 Newton-HSS 后 退 方 法 ,证 明了 它 的 全 局 收敛 性 . 杨 爱 利 、 伍 渝 江 
等 将 Newton-HSS 方法 进一步 推广 ,提出 了 一 类 求解 非 线性 方程 组 的 Newton-PSS 和 迭代 
法 . 除 此 之 外 , 当 弱 非 线性 代数 方程 组 F(x) =Ax—@ (x) =0 时 , 白 中 治 等 利用 Picard 3€ 
代 作 为 外 迭代 ,HSS 方法 作为 内 迭代 ,提出 了 求解 弱 非 线性 代数 方程 组 的 Picard-HSS 和 
非 线性 HSS-like 迭代 方法 ,这 两 类 方法 在 计算 的 过 程 中 均 不 需要 在 当前 迭代 步 x 处 计 
算 和 形成 Jacobi 4E E , AEA ETE HSS 方法 中 出 现 的 两 个 子 线性 代数 方程 组 的 系数 矩阵 
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不 随 迭 代步 数 的 变化 而 变化 . 和 Newton-HSS 方法 相 比 ,这 两 类 方法 大 大 减少 了 算法 的 
工作 量 , 提 高 了 求解 速度 . 


然而 当 弱 非 线性 代数 方程 组 
F(x) = Ax —6 (0 — 0 (6. 2) 
的 系数 矩阵 具有 下 述 形式 


A-—-W-riT 

矩阵 妈 T 分 别 为 实 正定 矩阵 和 正 半 定 和 矩阵 且 满 足 W —W. T" —T.6 : DOC> 为 定 
义 在 复线 性 空间 C" 的 一 个 开 凸 子 集 D 上 的 连续 可 微 函 数 , 且 非 线性 函数 更 (x) AY Jacobi 
EED’ GOTERE x € D 处 是 非 Hermite 和 负 半 定 的 . 此 种 情况 下 ,Picard-HSS 和 非 线性 
HSS-like 方 法 收敛 速度 会 大 大 降低 ,这 是 因为 在 Picard 迭代 的 每 一 步 利 用 HSS 方法 求 
解 时 ,会 出 现 两 个 系数 矩阵 分 别 为 aI HW R aI HiT 的 线性 代数 方程 组 ,对 于 后 一 个 方程 
组 ,往往 需要 利用 重复 法 求解 . 大 大 影响 了 求解 非 线 性 方程 组 的 求解 速度 . 

对 于 非 线 性 代数 方程 组 (6. 10 ,我 们 将 LHSS 方法 与 不 精确 Newton 法 相 结 合 ,并 在 
不 精确 Newton 方向 上 配 以 线性 搜索 策略 ,提出 了 具有 全 局 收敛 性 质 的 不 精确 Newton- 
LHSS 后 退 方法 ,对 于 具有 特殊 结构 的 弱 线 性 代数 方程 组 (6. 2) ,我们 提出 了 Picard- 
AHSS 方法 和 AHSS-like 方 法 以 及 Picard-MHSS 方法 ,对 于 每 种 方法 都 仔细 分 析 了 其 收 
敛 性 ,给 出 相应 的 收敛 定理 . 大量 的 数值 实验 证 明了 理论 的 正确 性 及 方法 的 有 效 性 . 


6.2 关于 Newton-LHSS 后 退 方 法 及 其 全 局 收敛 性 的 研究 


本 章 采用 的 范 数 均 为 欧式 范 数 . 记 BO (yl || y—x || <r} 是 一 个 以 x 为 中 心 ， 
以 r>0 为 半径 的 开 球 . B(x,r) 表示 B.r) WAE. AT 表示 A 的 转 置 . 


6.2.1 Newton-LHSS 后 退 方法 的 提出 


不 精确 Newton 方法 仅 具 有 局 部 收敛 性 质 , 如 果 使 用 LHSS 方法 求解 Newton 方程 ， 
并 在 其 不 精确 Newton 方向 s 上 配 以 线性 搜索 策略 , 则 可 以 得 到 具有 全 局 收敛 性 质 的 不 
精确 Newton-LHSS 后 退 方法 (NLHSSB 方法 ). 该 方法 可 以 简 述 如 下 . 
给 定 Xo s max € L041) tE (0,1) ERR a, EE BM {Ly} 0 0 Duis < Omar <1 tol770. 
如 果 || Fox) || >tol * mint || FG) || n) BU x, RB xa AP PETI: 
第 一 步 : 选取 € [0 ms] di. —0. 
对 于 10.1.2. 4 —1. H LHSS Wik 
H(x, dui —— SC day — FG), 
(af + S(x.)) din = (al — Hx) dint = F(x) 


R de ,使 得 | FG) HE Ga dia, lp || EC I 
BL H Gu) — TOF Gi) - FG" )SG T TU" FG" ) 


第 二 步 : 沿 着 不 精确 Newton 方向 di 配 以 线性 搜索 后 退 策略 ,具体 方法 如 下 : 


m68- Hu X hgnAGUBOLE GN UE 


(D $t di =de 57, =m. 
Q MF || FG, td) ll 1:0 3,51 || FG) I 
选取 O= [Omin + Omax] 
© d,—0d,.5,—1—00 —3)0. 
WEP: Xem =d. 
值得 指出 的 是 ,在 不 精确 Newton P , ARA x, 远离 非 线性 方程 组 (6. 1) 的 真 解 
时 ,不 宜 选用 过 小 的 强制 项 . 否则 ,很 可 能 会 导致 oversolving 现象 , 即 对 于 Newton 方程 
的 求解 达到 适当 的 精度 时 , | PCr.) || 就 会 得 到 最 好 的 下 降 . 若 这 时 再 对 Newton 方程 一 
味 地 去 精确 求解 , 则 可 能 会 造成 大 量 的 计算 开销 ,并 且 对 FG ‖ 的 下 降 也 起 不 到 好 


的 作用 ,更 有 可 能 会 导致 中 断 现 象 . 
本 书 选取 下 面 两 种 强制 序列 ,这 两 种 方法 现 已 得 到 广泛 而 成 功 的 应 用 . 
-1 二 V5 — |l FG) ll — ll FG +F Gad Ill 
选择 1. 给 定 5 二 一 2 一，'Y mh EL0,1), 选 取 pA Tro.) 4 
k—0,1,2,-- 


为 了 防止 me 过 小 ,进一步 计算 p , 即 采 取 “ 保 护 措施 ”: 
保护 措施 1. 如 果 ^ 0. 1, JU TE p= maxi hg, 
选择 2. 给 定 A€ [0,1] Hl p€ (1,2], EF. Vg € [0,10 ER p= 


| FG) I y 
À ,b—0,1,2,7 
( | FG) I 


保护 措施 2. 如 果 Aot 770. 1, I EE — max p Anos}. 

除 此 之 外 ,为 了 防止 出 现 oversolving 现象 ,文献 [143] 提 出 了 两 种 “保护 措施 ”, 即 
保护 措施 3: m 7 minUp equ) 

保护 措施 4: 如 果 pe <2e/ || FG) || , 则 令 peo. 8e/ || Fœ |l . 


6.2.2 Newton-LHSS 后 退 方法 的 全 局 收敛 性 定理 


为 了 得 到 NLHSSB 方法 的 全 局 收敛 性 定理 ,首先 做 出 如 下 假设 ,并 给 出 引 理 6.1. 
假设 6.1 AEB x €C'.F:DCC'—C"'fE— 4 JF SBR N CD 内 是 G- 可 微 的 . 


FOEN 中 是 连续 .正定 的 . 假设 F(x) 二 HH(x) 十 S(x), 其 中 HG) 一 二 (Cr 十 


F(x)" SQ =F UF) —F (0 ) 分 别 为 Jacobi 矩阵 F(x) 的 对 称 部 分 和 反对 称 部 分 ， 除 


此 之 外 ,假设 下 面条 件 成 立 : 
(A1) LHSS 方法 收敛 . 
(Az) (Lipschitz 条 件 ) 存 在 非 负 常 数 La 和 工 , 使 得 对 所 有 的 x,yE B(xo,r)CNo, 有 
|| Hi) —HCOQpl <Lallx—-yll. lSco—Sool < L, x y Il 
在 假设 6.1 下 ,应 用 扰动 引 理 ,我 们 得 到 引 理 6. 1. 
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引 理 6.1 在 假设 6.1 下 ,我 们 有 下 面 三 个 结论 成 立 : 
(1) |F œF o) || <L || x—y |] 3x B.L—L;-L.. 
(2) || Fœ || <L ll x—xp || +28 


d P 4 3g a if 
(3) 如 果 r< AF (x) 是 非 奇 异 的 ,并 且 满 足 E GO 7 |l SIL l| xx: T 


假设 6. 1 保证 了 引 理 6. 1 成立. 这 就 意味 着 F 是 Lipschitz 连续 的 ,Lipschitz 常数 


为 上 ,存在 一 个 正常 数 m HA HERA Adn) D= U ixl xx <r) 上 满足 
LF’ Ce)? || Sang. 利用 相关 文献 中 的 定理 ,我 们 有 下 面 的 两 个 全 局 收敛 定理 . 

定理 6.1 给 定 初 值 xxEc",tE(0,1). 假设 (x,) 是 由 NLHSSB 算法 得 到 的 序列 ， 
p 是 由 选择 1 和 所 有 的 保护 措施 给 定 . 假设 序列 {xe} 有 界 , 假 设 6. 1 成 立 . WBA (xs BK 
到 F(x) —0 的 一 个 根 x* , 且 有 下 面 两 个 结论 成 立 : 

(D 如 果 了 办 ~0, 则 NLHSSB 方 法 具有 强 - 超 线性 收敛 速度 . 

© 如 果 对 于 某 个 K,>0, WER pK, || Fx) | ^.9 NLHSSB 方法 具有 强 1 十 p 阶 
收敛 速度 . 

定理 6.2. 给 定 初 值 xxEC",E(0,1). 假设 (xi) 是 由 NHSSB 算法 得 到 的 序列 ,办 
是 由 选择 2 和 所 有 的 保护 措施 给 定 . 假设 序列 {xx} 有 界 ,假设 6.1 成 立 . 那么 (x) 收敛 到 
F(x) —0 的 一 个 根 x* , 且 有 下 面 两 个 结论 成 立 : 

(D 如 果 4 二 1, 则 NLHSSB 算法 是 强 -o 阶 收敛 的 . 

© 如 果 )=1, 则 NLHSSB 算法 是 弱 -p 阶 收敛 的 ;对 于 任意 pE[1,p) ,NLHSSB 算法 
是 强 -p 阶 收敛 的 . 


6.2.3 数值 实验 
考虑 下 面 的 二 维 非 线 性 反应 扩散 方程 


Cue +t) Fu, T us) e, (zy) EN 
Poe =0, (x,y) € OD 


Jt O= (041) x 0.0.20 为 9 的 边界 ,gq 为 正常 数 ， 应 用 五 点 差分 格式 离散 扩散 项 ,中 
心 差分 格式 离散 对 流 项 . h= ,Re 一 公分 别 表示 网 格 步 长 和 网 格雷 诺 数 ， 我 们 得 到 


N+1 
下 面 的 非 线性 方程 组 : 
» Tae 一 


u = nsus sun)" su; = (üa stzs* sun)" si = 1,2, 7, N 
这 里 ASTLOQIHOT, ,四 表示 克 罗 内 克 内 积 ,T- AT, 表示 三 角 和 矩阵 . 
T. = tridiag(tz ,ti ,ts) 
T, = tridiag(t, 0.45). 
t —4. t l1—Re,. ú 1 十 Re 
选择 下 面 4 种 强制 序列 : 
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O 选择 1, 保 护 措施 1, 保 护 措施 3. 保 护 措施 4. 
@ 选择 2, 保 护 措施 2, 保 护 措施 3 ,保护 措施 4. 
© 对 于 所 有 的 &, 选 取 o, —0. 01. 


FG) Il — IF OF Gade , 
ve [Fœ Tl ,R—0,1,2,-, 
保护 措施 1 ,保护 措施 3 ,保护 措施 4. 


注意 ,第 一 种 强制 序列 和 第 四 种 强制 序列 计算 方法 在 分 母 上 稍 有 不 同 , 因 为 F(x) 
比 || Fox) | 更 接近 于 当前 迭代 的 值 ,所 以 在 强制 序列 的 计算 中 采用 此 值 是 合理 的 . 对 
于 算法 中 参数 的 选择 RIED OD .田中 ,选取 p=0. 5, hm =0. 9. 在 线性 搜索 中 ,选取 
1—107*,0,,—0. 1,0 一 0. 5. 在 算法 NLHSSB 中 ,选择 0€ [Onin Om] f 49 g (0) = 


| Fon todo ll RH. Am Lp 15. Witz He EPI qm Los m See C1951, 


N=20. K 6.1—6. 4 给 出 了 在 4 种 强制 序列 的 选择 下 , NLHSSB 方法 和 NHSSB 方法 的 
比较 结果 . 从 表 6. 1 一 表 6. 4 中 可 以 看 出 ,NLHSSB 方法 在 CPU 时 间 和 迭代 步 数 上 都 优 
于 NHSSB 方法 . 


表 6.1 强制 序列 1 情形 ,两 种 方法 的 比较 


方 法 外 迭代 数 内 迭代 数 | FG) | CPU 时 间 
2 15. 541 0. 1402 
3 6. 385 0. 2804 
7 2. 2402 0. 5608 
9 0. 7663 0. 9213 
Newton-HSS 方法 8 
13 0. 1257 1. 422 
20 0. 005 84 2.1932 
32 3.5753X le 一 5 3.3448 
50 1.0917 X1e—8 4. 7669 
2 9. 3892 0. 2303 
2 1. 7029 0. 4607 
Newton-LHSS ; 2 0. 3509 0. 7110 
方法 3 0. 006 83 1. 3319 
4 2.0841X le—5 2.073 
6 1.1206 X le 一 8 2.6738 
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表 6.2 强制 序列 2 情形 .两 种 方法 的 比较 
A & Shik AR RL Vi (CR | FG |l CPU 时 间 
2 15. 541 0. 1803 
3 6.385 0. 3205 
7 2.399 0. 5808 
9 0. 8315 0.9914 
Newton-HSS 方法 8 
13 0.1792 1.4721 
20 0.001 58 2.1331 
32 0.000 31 3.0744 
50 4. 871 X 1e—7 4. 2561 
2 9. 3892 0. 2303 
2 1. 7029 0. 4607 
Neros taSs ; 2 0. 3509 0. 7110 
方法 3 0. 006 83 1.382 
4 2. 0839 X 1e—5 2.1331 
5 1.8612X1e—8 2.6538 
56.3 强制 序列 3 情形 ,两 种 方法 的 比较 
方 法 外 和 迭代 数 内 迭代 数 | F(x.) |l CPU 时 间 
47 9. 3732 1.5422 
65 1. 6263 3. 8856 
Newton-HSS 方法 5 90 0. 2305 7. 2004 
89 0. 002 10. 485 
61 2. 0245 X 1e—7 12. 238 
8 9. 3732 0. 8702 
9 1. 6263 2. 2733 
fpem 5 11 0. 2305 4. 6267 
11 0. 002 7.0301 
8 1.2112X1e—7 8. 2819 
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表 6.4 强制 序列 4 情形 ,两 种 方法 的 比较 


» * AIR Vj3k EBC | FG |l CPU 时 间 
2 15.541 0. 1603 
3 6. 385 0. 3105 
6 2. 2402 0. 5808 
5 1. 2216 0.8111 
Newton-HSS 方法 9 8 0. 4194 1.1016 
12 0.072 01 1.5823 
17 0. 0052 2. 2532 
23 0. 000 14 3.0544 
33 7.0045 X1e—7 4. 0458 
2 9. 3892 0. 2704 
2 1. 7029 0.5107 
2 0, 3509 0. 7510 
2 0. 0555 1. 1016 
a 9 2 0, 0081 1.392 
2 0. 0013 1.6924 
2 0. 0002 1. 9428 
2 3. 0765 X le—5 2.1431 
3 7. 3496 X 1e—7 2.4836 


6.3 Picard-AHSS 方法 及 其 局 部 收敛 定理 


假设 矩阵 A E CAECIS LIE REB PE. PR : DCC" 一 C" 为 定义 在 复线 性 空间 
C" 的 一 个 开 凸 子 集 D 上 的 连续 可 微 函数 ,上 且 ®@ (x) 的 Jacobi EED’ GO ZEA x* ED 处 是 
dE Hermite 和 负 半 定 的 . 求解 弱 非 线性 代数 方程 组 (6. 2) 

我 们 这 里 说 的 弱 非 线性 方程 组 指 的 是 在 某 种 范 数 意义 下 ,线性 项 Ax 强占 优 于 非 线 
性 项 ® (x). 一 个 矩阵 GEC”" 为 负 半 定 的 、 负 定 、 正 半 定 和 正定 的 当 且 仅 当 它 的 Hermite 


部 分 H(G) 空 证 (G 二 G7) 分 别 是 负 半 定 的 、 负 定 、 正 半 定 和 正定 的 ,一 个 矩阵 G 的 


Hermite/ 反 Hermite 分 裂 指 的 是 对 于 任何 一 个 矩阵 GE CH AT LA BA A 
Hermite/ 反 Hermite 矩阵 之 和 , 即 
G = H(G) + S(G) 


def 
其 中 H(G) G 的 Hermite BF. s (7 6-6 A G 的 反 Hermite 部 分 . 
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6.3.1 AHSS 方法 


考虑 下 列 大 型 稀 玻 线性 方程 组 
Ax =b 
其 中 AEC"" 为 大 型 稀疏 矩阵 ,x BEC". 
非 对 称 的 HSS 分裂 迭 代 格 式 的 定义 如 下 : 
AHSS RAE: 给 定 一 个 初 值 xO EC”, 对 于 k 二 0,1,2,… ,直到 {x } 收 敛 ,计算 
GI + HDx (€?) = (aI — Sx? +b, 


(6. 3) 
(BI 4- S)x'*? = (BI — Hx C) +b, 
其 中 a,8 为 给 定 正 整数 ,LI 为 与 矩阵 4 同 阶 的 单位 和 矩阵. 
上 述 迭 代 过 程 可 以 等 价 于 
x^? —M(G.Dx^? + GG. b 
=M(a, fp x 十 Y MG PG. Db (6.4) 


l= 0,1,2, 
其 中 


Ma.) = (BI 4- S) ! (BI — H) (al +H)" (al —S) 
G(a.B) = (a+) (BI+S) (aI +H)" (6.5) 
事实 上 ,AHSS 和 迭代 过 程 也 可 以 看 成 是 由 系数 和 矩阵 A 按 下 述 分 裂 
A = Bla, P) — Cla.) 


得 到 ,其 中 
Bs) = 3 (al +H) (BI +8) (6.6) 
Ca. Del H) (BI — H) («I +H)” (aI —S) (6.7) 
Ma,p) = Ba.) Cf) GG f) = BGB? (6. 8) 


MAE ACC" AAR Hermite 正定 矩阵 时 , 若 a,B 满足 一 定 的 条 件 ,AHSS 方法 的 迭代 矩 
阵 的 谱 半径 小 于 1. 杨 爱 利 、 伍 渝 江 等 对 该 方法 的 收敛 性 进行 了 更 加 细致 的 讨论 ,并 详细 
给 出 了 两 个 参数 的 计算 方法 . 数值 实验 证 明了 AHSS 方法 是 较 HSS 方法 更 加 有 效 稳健 
的 迭代 方法 . 


6.3.2 Picard-AHSS 方法 


下 面 我 们 给 出 求解 非 线性 代数 方程 组 的 Picard-AHSS 迭代 方法 . 
对 于 非 线性 方程 组 (6. 2). 当 线性 项 Ax 强占 优 于 非 线 性 项 ® (x) 时 ,我 们 可 以 利用 
Picard 迭代 
Ax" = G(x"), e= ay (6.9) 
作为 外 迭代 计算 (6. 2) 的 解 . 当 和 矩阵 4 DAKUN TE XE AB IAE. TES. k Ib 3E IC x PURIS 
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我 们 利用 AHSS 迭代 格式 
Bla, Dx? = Cla, p) xt? +O (x?), 1—0,1,--.4 —1,k = 0,1,2, 
(6. 10) 
作为 Picard 迭代 的 内 迭代 算法 ,这 里 Bap) 和 CCo. ils C6. 6) SX C6. 05E Mea. 为 
给 定 的 正常 数 ,{2e} 电 为 一 指定 的 正 整数 数列 ,在 第 & 步 Picard 迭代 时 ,选取 内 迭代 初 值 
x^? 二 x 中 ,由 此 ,我 们 可 以 构造 出 一 类 新 的 方法 , 称 为 Picard-AHSS 方法 . 
Picard-AHSS 迭代 方法 : HH: DCCC’ 为 一 连续 可 微 函 数 ,和 矩阵 A 为 一 正定 矩 


阵 ,H i (A+A"* ),S i (A—A" ) 分 别 为 4 的 Hermite 部 分 和 反 Hermite 部 分 . 给 定 


初始 估计 x" € D 以 及 正 整 数 序列 { 之。 对 于 & 一 0,1,2,…', 利 用 下 面 的 步骤 计算 
x**» ,直到 近似 解 x 中 达到 给 定 精度 . 
DOI= 
© 对 于 /二 0,1,2,… ,一 1, 求 解 下 述 线性 方程 组 得 到 x7?， 
(al + HDx (^3) = (af — S)x*? +® (x®), 
(BI + S)x ^? = (BI — HxH) + 6 (x), 
其 中 a. 为 给 定 的 正常 数 . 
Q 4 oye, 
由 上 述 和 迭代 方法 可 知 , 当 参数 a — i. Picard-AHSS 退化 为 Picard-HSS 方法 ,也 就 
是 说 ,Picard-AHSS 方法 实际 上 是 Picard-HSS 方法 的 一 个 推广 , 它 既 保持 了 Picard-HSS 
方法 的 优点 ,不 需要 在 当前 迭代 步 x 处 计算 和 储存 Jacobi 4E PE, HL Tk TAD p 3x (se 
中 ,出 现 的 两 个 子 线性 方程 组 的 系数 矩阵 均 为 常数 和 矩阵 ,又 将 Picard-HSS 方法 进一步 优 
化 ,在 内 和 迭代 过 程 中 ,对 于 不 同 的 子 线性 代数 方程 组 的 系数 矩阵 选取 不 同 的 参数 a 和 B. 
提高 了 非 线 性 代数 方程 组 的 求解 效率 . 


6.3.3 Picard-AHSS 方法 收敛 定理 


接 下 来 我 们 给 出 Picard-AHSS 方法 的 局 部 收敛 定理 . 
首先 引入 下 面 的 记号 . 对 于 任意 的 向 量 x CC" ARE XO, FR gg CIS EE A 
及 其 诱导 范 数 
lel S eoru sxll,, Exil S ear sxcar s? Il. 
我 们 有 下 面 的 局 部 收敛 定理 : 
定理 6.3 Wx ED 是 非 线性 方程 组 (1) 的 解 . DC C" C" fk x" 的 开 邻 域 


NoCD EG TEAD’ (x) 连续 . 假设 矩阵 A Wy TRL LE HE PE = (ALA » 


(6.11) 


S— 1 (A—A' ) 分 别 为 A 的 Hermite 部 分 和 反 Hermite 部 分 . 记 
Kap) = || Mel. a= || ATI], v= | A? 9,' (x") Il 


EAEN AHSS 迁 代 满足 收 全 条 件 且 0<1, td (12) [maa | met x 
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的 一 个 开 邻 域 NCNoe ,使 得 对 于 任何 初 值 x"E N RIETI IE EZ 3 1, 6—0.1.2, 7. 
由 Picard-AHSS 迭代 方法 生成 的 序列 (x } 亿 "总 是 适 定 的 且 收 和 敛 于 非 线 性 方程 组 (6. 2) 
的 解 x* . 除 此 之 外 ,有 下 式 成 立 : 
limsup || x® —x* || 二 之 v 十 (1 十 vj0(asB%*， 其 中 1 = liminfl, 
特别 地 ,如果 lim, = 0° Jf A XE fU JUR eR A R 线性 收敛 的 , 且 R 收敛 因子 至 多 为 
v, 即 
limsup Il x? —x" Il * <v 


WEBB 由 Picard-AHSS 和 迭代 格式 (6. 11) ,经 计算 可 以 得 到 


4-1 
x? = Mla, Phx +Y Ma DGR B(x), k=0,1,2,.%, (6.12) 
j=0 
其 中 Ma, DM Gap) 由 (6.5) 定 义 . hF Fa ) 二 0, 故 有 
471 
x' =M Gofyix* + XMa. pG, B(x"), k—0,1,2,. — (6.13) 
j=0 


式 (6. 12) 减 去 式 (6. 13) ,我 们 得 到 

ġa 
x) — x" = M(B (x'? —x* ) + >) Ma pG aM (x) — a )], 
j=0 
k=0,1,2,. (6. 14) 
id 
def , 
E(x,x" ) =@ (x) —0(x'0—0 (x' )(x— x") 
则 (6. 14) 式 可 以 等 价 地 写成 
x"? —xy* 
4) 4-1 


= [M Gap) + >) M GGG D (x" )] G9 — x") + >) M aG DE (x? a ) 
T md j=0 


= [M(a. f^ + UI Ma At ®'(x* JG? — x") + Q — MG ATE ox" ) 
= [A7 ®'(x* ) -- Mia (1 — A? ®(x" )) Ja? —x* ) + G— MG 4A? EG? x"), 
(6.15) 


—x' || «EI A? 'G^»5l + IMG» |] 4+ IA? e'o»1»]Ix*? —x 1 + 
A++ || MG: I^»1 A? II EGIT? 5x") I (6. 16) 

—[v4- 1+ 90 ^ ] || x? —x* || +a 43-6Ca 0^) || Ex ix? ) II 
HAAS, DC CC" 在 rz ED 的 某 个 开 邻 域 NoCD AG ARMO’ E x” € D fib 
连续 可 以 得 出 : 对 于 任意 给 定 的 s 盖 0 和 所 有 的 x € N, 都 存在 x 的 一 个 开 邻 域 NCNo， 


1 
使 得 上 E(x,x* ) || ell x—x* ||. 由 已 知 条 件 v<1， >| (Ics) incon.» 可 知 


vt (14 v)8G.B) <1 
进一步 假设 e 满足 
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v+ A+, PD +a(l+O(a,B)e — 1 
则 当 a € N 并 且 1, — limZ, 时 ,根据 式 (6. 16), 利 用 数学 归纳 法 ,我 们 立刻 得 到 
| ECX9,x* y el xt —x*: | 
l| x? —x* || <[v+ C1 3- 00a .)^ +a + 0Ca 0^)2e] || x” —x* || (6.17) 
并 且 对 于 所 有 的 b—0.1.2.-- MBA x? € N mr. 由 式 (6.16) 和 式 (6.17) 可 以 得 到 


limx® = x", 


lim sup || x^ — x" | * <v (4- v6.) Hall + Ca fe 

由 6270 的 任意 性 和 R 收敛 因子 范 数 独立 可 知 ,由 Picard-AHSS 35 (6r p^ ^E AY Fe 9] 
{x }AY R 收敛 因子 至 多 为 v 十 (1 十 v)0(a,B)%. Hab lim4 二 ,由 式 (6. 17) 可 知 , 序 
列 x? FER REKA FLUC DE AE ev. 

定理 6. 3 表明 ,Picard-AHSS 和 迭代 的 收敛 速率 实质 上 由 wv 和 0(a,B) 决 定 . 通常 ,v 和 
O(a.) 越 小 ,Picard-AHSS 迭代 的 收敛 速率 越 快 . 当 弱 非 线性 方程 组 中 的 系数 矩阵 A 为 
病态 时 ,0(a,B) 接 近 于 1,Picard-AHSS 迭代 的 收敛 速率 会 变 得 非常 慢 以 至 于 内 迭代 步 数 
很 多 ,在 这 种 情况 下 ,可 以 考虑 使 用 预 处 理 的 AHSSCPAHSS) 和 迭代 代替 AHSS 迭代 方法 ， 
这 样 会 使 得 0(a,B) 小 一 些 , 不 精确 的 Picard 和 迭代 的 收敛 速率 更 快 一 些 . 

本 节 最 后 ,为 了 方便 数值 计算 ,我 们 将 Picard-AHSS 和 迭代 写成 等 价 的 残 量 校正 

Picard-AHSS 的 残 量 校正 形式 : 设 ® DC CC 为 一 连续 可 微 函 数 ,矩阵 4 为 大 型 
稀疏 正定 矩阵 ,H La FA*),S : (A—A" ) 分 别 为 4 的 Hermite 部 分 和 反 Hermite 
部 分 . 给 定 初始 估计 x € D 以 及 正 整数 序列 {4} 守 ,对 于 A=0,1,2,…, 利 用 下 面 的 步 
又 计算 x*+? ,直到 近似 解 x? 达到 给 定 精度 . 

[6] so :; 00,5? :=@ (x )—Ax('? 

@ 对 于 /==0,1,2,…,l 一 1, 利 用 下 面 的 迭代 格式 求解 线性 方程 组 得 到 snm : 


1 
(al -- IDs = (gp — ss +b”, 


T 


(6. 18) 


(BI -- S)s**? = (BI— ms C» +b”; 


@ ^ xD : = se. 
6.3.4 非 线 性 AHSS-like 迭代 方法 及 其 收敛 性 定理 


本 节 我 们 提出 另 一 种 求解 (6. 2) 的 和 迭代 方法 ,叫做 非 线性 AHSS-like 迭代 方法 . 这 种 
方法 是 基于 非 线 性 不 动 点 方程 
(al Hox = (GE — S)x-- (x). 和 C8EL -- S) x = (BI— Hx - 6 (x) 
得 到 的 . 它 既 可 以 保持 Picard- AHSS 方法 的 优点 ,又 可 以 避免 使 用 显 式 的 内 迭代 方程 . 
非 线性 AHSS-like 迭代 方法 : 给 定 初始 估计 x^? € D. 利用 下 面 的 步骤 计算 x“ > , 直 
到 近似 解 zx 达到 给 定 精度 
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(al +H)x (=) = (al —S)x® + (2) 
(BI -- S)? = (BI —H)x (#2) +ø (x €) 
这 里 ,a,B 为 给 定 的 正常 数 , 工 表示 单位 矩阵 . 
下 面 给 出 非 线 性 AHSS-like 迭代 方法 的 局 部 收敛 定理 . 
定理 6.4 Kx HARARE BAG. DWM. D: DC C" 在 x* 处 下 可 导 . 矩阵 


4 为 一 正定 矩阵 ,五 L (A*A' )，S la A’) 分 别 为 A 的 Hermite 部 分 和 反 


(6.19) 


Hermite 部 分 . id 
M(a.Bix' ) = (BI+S) (BI —H+@'(x" )) (al +H) (al —S+@'(x" )). 
如 果 pCM(a,B;x* <1, Il x* ED 为 非 线性 AHSS-like 迭代 的 一 个 吸引 子 . 
证 明 首先 ,定义 函数 
U(x) = (GI +H) (Cal — x 6)». 
Vix) = (I +S) (al — Hox - (x)) 
AW (x) —V*U GO : —VQU GOD , 则 非 线性 AHSS-like 迭代 格式 (6. 19) 等 价 于 
y*"—4.*(Q?), Ek—0;152,- 
由 Ostrowski 定理 ,如 果 oC (xt 00 —1.. IBA x" 为 非 线 性 AHSS-like 和 迭代 序列 的 吸引 
T. 下 面 我 们 证 明 PCR ct 0) —1. 
因为 x"* 为 非 线 性 方程 组 的 解 ,所 以 有 
x' =U(x"), x' =V(x") 


(6. 20) 


和 


U' (x' ) = (aI +H)" (aI —S+@’ (x' )) 
V’ (x' ) = (BI - S)! (BI—H--9' (x' )) 
显然 成 立 . 利用 链 式 求 导 法 则 ,我 们 可 以 得 到 
Y(x)-—V (x' JU (x' ) = M (e fix! ) 
由 此 ,定理 结果 显然 成 立 . 证 毕 . 
下 面 的 定理 给 出 了 pCM(a,B;x* <1 的 充分 条 件 . 
定理 6.5 (Bit. DCC" 一 C" 满足 定理 (6.4) 的 条 件 . 引入 下 列 记号 
ô = max (ll 6' (x' ) (BI --S)? l3. | Ø’ (x* ) (aI +H) ls) 
n= |l (ef 一 S) (BI+S) ll; 0a.) = || Map) I 
其 中 矩阵 M(a,B) 如 (6.5) 所 定义 . 若 内 迭代 AHSS 迭代 满足 收敛 条 件 且 
2(1 — a.) 


« : > 
E | ») | J(53) 4 (a. f) —1) 
TW o(MCa.Bsx" 2) 1. Bl x* ED 为 非 线 性 AHSS-like 和 迭代 序列 的 吸引 子 . 
证 明 经 过 计算 ,可 以 得 到 下 面 的 等 式 
(BI - S)M (a«fx" ) (BI - S) 
= (BI +S)M(a.8) (BI --S)? + (BI — H) (aI +H)? 9' (x ) (BI -S)? 十 


(6.21) 


(6. 22) 


mm 
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©’ (x' ) (al -- H)! (aI — S) (BI - S)? +0’ (x' ) (a1 +H)" 6' (x' ) (BI +S)” 
因此 


M (a.83x" ) || = || (BE+S)M (asB;x* ) (GE - S) |; 
< || (BI 2- S)M(.) (BI 4-S)? |; + ll (8I — H) (a1 +H)? 9' (x* ) (8E +S)7 ll; + 
©’ (x' ) (aI +H)" («aI — S) (BI - S)? |;d (6. 23) 


O’ (x' ) (GI +H)" 6' (x' ) (BI - S)" ll; 

Mta.) || + || (8E — H) (aI -- H)? ll; || Ø’ (x^ ) (CT 十 S) ll; + 
$' (x' ) (al +H)" ll; ll (e1 —S) (81 - S)! li 

$' (x' ) (aI +H)” |l; || Ø’ (x' ) (gp - S)? ll; 


IN 


< 0G. Hit pcs 
利用 定理 中 条 件 (6. 210 ,我 们 有 
(aD TP 一 1 


因此 ,我们 有 
p (M (asBsx* )) < ll M (a. Bx" ) ll <1 
x' € D 为 非 线性 AHSS-like 迭代 序列 的 吸引 子 . 
注 6.1 在 证 明定 理 6. 5 中 公式 (6. 23) 的 最 后 一 个 不 等 式 中 ,我 们 用 到 了 
| (BI+-H) (aI +H) || 2, || (al —S)(BI+S) || ;—1. 这 是 因为 AHSS 迭代 收敛 ,a,B 需 要 
在 一 定 范围 内 取 值 ,而 这 个 范围 正好 保证 了 此 式 成 立 . 


6.3.5 数值 结果 


为 了 验证 Picard-AHSS 方法 和 非 线性 AHSS-like 迭代 方法 的 可 行 性 与 有 效 性 ,我 们 
考虑 下 面 的 二 维 非 线性 偏 微分 方程 
— Cus us) Faltis Fu) — us e* - sin /1o-u Fui. lay) EO 
u(x,y) — 0. (x.y) € IN, 
其 中 Q==(0,1)X(0,1),90 HKE O 的 边界 . 将 求解 区 域 进行 等 距 分 割 ,z,y 方向 步 长 


HH h= 利用 迎风 差分 格式 ,可 以 得 到 下 面 的 非 线 性 方程 组 


F(x) —Ax—$(G0—0, 其 中 n=NXN 
在 实际 计算 中 ,选取 x" =0. 外 Picard 迭代 停止 标准 取 为 


外 dx? lls 
I Fx?) 


内 HSS fll AHSS 和 迭代 的 停止 标准 取 为 


F' (xt) s@@ + F(x) |; 
| F(x”) l = 


Hep n 表示 每 次 Picard EAC JH EI AHSS 的 和 迭代 步 数 . p 为 强制 项 ,用 来 控制 内 选 
代 求解 的 精确 程度 . 在 本 试验 中 ,选取 因为 常数 , 记 为 思 对 于 内 和 迭代 AHSS 和 迭代 中 出 现 


= 
N+1’ 


<= lo* 
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的 两 个 子 线性 代数 方程 组 ,我 们 利用 LU 分 解 进行 求解 . 

首先 ,对 于 不 同 的 N ,不同 的 g, 我 们 给 出 HSS 方 法 和 AHSS 方法 中 参数 的 值 和 对 应 
谱 半径 的 值 的 大 小 的 比较 . 在 这 两 类 方法 中 参数 的 选取 均 是 使 得 相应 谱 半径 上 界 取得 最 
小 时 参数 的 值 . MÆ 6. 5 中 可 以 看 出 ,利用 AHSS 和 迭代 方法 求解 线性 代数 方程 组 (6. 9) 是 
收敛 的 , 且 和 迭代 矩阵 的 谱 半径 明显 小 于 HSS 方法 中 系数 矩阵 的 谱 半 径 . 


X 6.5  Picard-AHSS 和 Picard-HSS 方法 中 的 最 优 参数 的 比较 


N 8 10 12 15 20 
a 1.4525 1.1735 0. 9856 0. 7956 0. 6029 
Picard-HSS 
p(M(a)) 0.7002 0.7486 0.7835 0.8207 0. 8606 
q=10 a 0. 0930 0. 0646 0. 0472 0.0317 0.0122 
Picard-AHSS B 0. 5800 0. 3892 0. 2798 0. 1854 0.0706 
e(M(a,f)) 0. 3424 0. 3652 0. 3792 0. 3718 0. 4078 
a 1.5370 1. 2201 1. 0139 0. 8108 0. 6097 
Picard- HSS 
e(M(a)) 0. 7002 0. 7486 0. 7835 0. 8207 0. 8606 
q=20 a 0. 2830 0. 1951 0. 1419 0. 0941 0. 0556 
Picard-AHSS B 0. 7608 0. 5156 0. 3721 0. 2475 0. 1446 
(OM Ga)» 0. 4907 0. 5271 0. 5504 0. 5722 0. 5916 
a 1.6214 1. 2666 1. 0422 0. 8261 0. 6164 
Picard-HSS 
e(M(a)) 0. 7003 0. 7486 0. 7835 0. 8207 0. 8606 
q=30 a 0. 5057 0. 3466 0. 2511 0. 1673 0.0977 
Picard-AHSS B 0.9581 0. 6531 0. 4732 0. 3159 0.1852 
e(M(a.B)) 0. 5569 0.5997 0. 6281 0.6557 0.6812 


下 面 比 较 内 迭代 步 数 和 CPU 时 间 . 此 处 只 给 出 Picard-AHSS, Picard- HSS, 4E ££ PE 
HSS-like 和 非 线性 AHSS-like 四 种 方法 的 效率 对 比 结果 . K 6.6 一 表 6. 8 分 别 给 出 了 四 
种 方法 内 迭代 数 IT im 、 外 和 迭 代数 IT 以 及 CPU 时 间 的 比较 . 从 表 6. 6 一 表 6. 8 中 可 以 看 
出 ,与 Picard-HSS 和 非 线性 HSS-like 方法 相 比 ,Picard-AHSS 和 非 线 性 AHSS-like 迭代 
在 内 、 外 迭代 数目 以 及 和 迭代 时 间 上 ,都 具有 一 定 的 优越 性 . 


表 6.6 当 〗 —0. 1.9—10 时 4 种 方法 迭代 步 数 和 和 迭代 时 间 的 比较 


N 8 10 12 15 20 
I 33 41 48 58 74 
非 线性 HSS-like 
CPU 0.9672 | 3.1980 | 8.3461 | 29.437 82. 883 
"m 9 9 9 9 9 


非 线性 AHSS-like 
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805 。 增 量 未 知 元 及 其 预 处 理 和 迭代 算法 
续 表 
N 8 10 12 15 20 
IT. 31 38 45 54 70 
Picard-HSS I Tu 4 4 4 4 4 
CPU 3.9312 | 13.026 | 32.745 | 111.17 236.28 
Jure 9 9 9 9 9 
Picard-AHSS [i S 4 4 4 4 3 
CPU 1.17 2.8392 | 6.5052 | 18.159 30. 716 
36.7. ig —0.1.q—20 R} 4 BRAT IAD AVIA ACER I8] B Ee t 
N 8 10 12 15 20 
IT 31 39 46 56 73 
非 线 性 HSS-like 
CPU 0.9360 | 3.0264 | 8.0497 | 27.987 82. 181 
IT 13 14 14 14 14 
非 线性 AHSS-like 
CPU 0.4680 | 1.0920 | 2.3712 | 7.1136 15.725 
IH 30 37 43 53 69 
Picard-HSS IT 4 4 4 4 3 
CPU 3.9156 14. 04 31.871 | 109.61 232, 660 
Ia 13 14 14 14 15 
Picard-AHSS IT; 4 4 4 4 3 
CPU 1.7004 | 4.5240 | 10.343 | 22.854 50. 295 
36.8. Yq —0.1.q0—30 时 4 种 方法 和 迭代 步 数 和 迭代 时 间 的 比较 
N 8 10 12 15 20 
IT 29 36 44 54 71 
JEPE HSS-like 
CPU 0.9048 | 7.7924 | 7.7376 | 27.097 89. 903 
IT 16 18 19 20 20 
非 线 性 AHSS-like 
CPU 0.5304 | 1.5288 | 3.2916 | 10.171 25. 241 
Ih 28 34 4l 51 67 
Picard-HSS fT. 3 3 3 3 3 
CPU 2.8263 | 8.8765 | 23.634 | 82.930 227.418 
IT, 3 3 3 3 20 
Picard-AHSS Hus 4 4 4 4 3 
CPU 1.4508 | 4.1184 | 9.7189 | 34.726 66.971 
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6.4 一 类 弱 非 线性 方程 组 的 Picard-MHSS 迭代 方法 


6.4.1 Picard-MHSS 方法 及 其 局 部 收敛 定理 


考虑 弱 非 线性 代数 方程 组 (6.2). 4@: DCC" 一 C" 为 一 个 常数 向 量 , 即 ®@ (x) = 
b 时 , 非 线 性 代数 方程 组 (6.2) 退 化 为 
Ax 一 5，x EC" (6. 24) 
HEP A-—W-HiT ABE W.T 4r 9I SE IE oe PE AS E OF s AB Ie Hii i£ W* =W, T — T. 
修正 的 HSS 迭代 C(MHSS) 方 法 定义 如 下 . 
MHSS 和 迭代 方法 : 给 定 初始 估计 x" € DXF 10.1.2. 利用 下 面 的 格式 计算 
x^*? ,直到 满足 停止 条 件 : 


GI +W)x(*2) = (af —iT) x +b 
(6. 25) 


(al + Tx ^? = (al +iW)x (^2) — ib 
其 中 a iE CE SE A 同 阶 的 单位 矩阵 . 
由 于 MHSS 迭代 可 以 改写 为 
L 
x) = M(@)x® 4- GG)b = MG)! x? + $3 M (GYGG)b,. 1 = 0,1,2, 


j-o 


(6. 26) 
其 中 
M(a) = (al 4- T)? (aI + iW) (al +W) (aI — iT) 
Gla) = (1— Da (al -- T) ! (aI --W)! (6. 27) 
若 记 
B(a) = Lor +W)(al +T) 
Cla) = E ‘(al + iW) (al — iT) (6.28) 
显然 成 立 
A = B(a) — C(a) 
经 过 计算 ,可 以 得 到 
M(a) = Bla) iC(a), Gla) = Bla) (6. 29) 


因此 ,MHSS 迭代 格式 也 可 以 看 作 是 由 和 矩阵 AY PA = BC) —CCa) 13 8] RR PE 
B(a) 也 可 以 看 作为 矩阵 4 B8 TUAE FEAR PE UY flt MHSS 预 处 理子 . HARE A=W + iT EE 
W.T r3 Jy Sc iE XE E P RISE IE 2E SE E PE EL ii WT —W.T* — T BE. H P8 S AGIEWI T : 
对 于 任意 正 数 a MHSS 方法 的 迭代 和 矩阵 M(a) 的 谱 半径 oCM(Ca)) 小 于 1, 即 MHSS 方法 
对 于 任意 正 数 a 和 初始 估计 x” 均 收敛 到 线性 方程 组 (6. 24) 的 精确 解 . 大 量 数值 实验 证 
明了 MHSS 方法 对 于 求解 复 对 称 线性 代数 方程 组 是 有 效 稳 健 的 . 对 于 非 线性 方程 组 
(6. 2) , 当 线 性 项 Ax 强占 优 于 非 线 性 项 ® (x) 时 ,我 们 可 以 利用 Picard 迭代 


81 
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Ax" d 二 0 

HAG. DWI. ERE A=W HiT, EE W, T 分 别 为 实 正定 和 实 正 半 定 和 矩阵 且 满 足 
W! —W.T'—T Ht. f] MHSS 迭代 对 上 述 方程 组 进行 近似 求解 , 即 利用 MHSS 迭代 
格式 

B(a)x*"? = C(a)x*? + @(x), k—0.1,2,-- (6. 30) 
作为 Picard 迭代 的 内 迭代 算法 ,这 里 BCa) FII Cla) 如 式 (6. 28) 定 义 ,a 为 给 定 的 正常 数 ， 
{Le} eo A FB EY TE SB HES k 步 Picard 迭代 时 ,选取 内 迭代 初 值 x ^? 一 x ,由 
此 ,我 们 提出 一 类 新 的 方法 , 称 为 Picard-MHSS 方法 . 

Picard-MHSS 迭代 方法 iø: DC C"— C" 为 一 连续 可 微 函 数 ,和 矩阵 A=WH+IT, HE 
HIERE WT 分别 为 实 正定 矩阵 和 正 半 定 矩阵 且 满 足 W* —W.T" — T. 给 定 初始 估计 
x € D 以 及 正 整数 序列 (ho ,对 于 k= 二 0,1,2,…, 利 用 下 面 的 步骤 计算 OP ,直到 近 
似 解 x* 达 到 给 定 精度 . 

D Rr =”, 

© 对 于 10.1.2. 4, —1 RA FRR Ey EAA 8B x77 ， 

(al +W)x (^t) = (al —iT) x? +O (x) 
(al + T)x ^? = (al +iW)x (#2) —id (x?) 
其 中 a 为 一 给 定 的 正常 数 . 

@ F xFtP =x, 

注 6.2 由 上 述 算 法 可 知 ,Picard-MHSS 方法 一 方面 保持 了 Picard-HSS 方法 中 的 优 
点 ,不 需要 在 当前 迭代 步 x 处 计算 和 存储 Jacobi 矩阵 , 且 在 所 有 内 迭代 过 程 中 ,出 现 的 
两 个 子 线性 方程 组 的 系数 矩阵 均 为 常数 矩阵 , 另 一 方面 ,又 避免 了 在 内 迭代 过 程 中 出 现 系 
数 和 矩阵 为 复 矩 阵 aI HiT 的 线性 代数 方程 组 . 

下 面 ,我 们 分 析 Picard-MHSS 方法 的 局 部 收敛 性 质 . 首先 引入 下 面 的 定义 . 对 于 任 
意 的 向 量 x © CREE X € C" ,我 们 定义 向 量 范 数 及 其 诱导 范 数 


1z1 苦 1 cer Della. UX S d rmx car m2 I 
可 以 得 到 下 述 局 部 收敛 定理 : 
定理 6.6 Bx ED 是 非 线性 方程 组 (6.2) 的 解 . D: DOCCE x^ RF SBI 
NED EG 可 导 且 ® ‘(x* EB. 假设 A—W--iT. WT —W.T? — T. 矩阵 W,T 分 别 
为 实 正定 和 正 半 定 矩阵 . 记 
Ha) = || M@ Il, w= |] ATI. B= | Aag cat) ll 


(6, 31) 


#5 <1 sto>| (155) [into |, 则 存在 x 的 一 个 开 令 域 NCN,, 使 得 对 于 任何 
初 值 x € N 和 任何 正 整数 序列 1 .k= 二 0,1,2,… ,由 Picard-MHSS 迭代 方法 生成 的 序列 
{Xe bee BA AE Te FE AY AO PAE ZR EENT FEH A x” 除 此 之 外 ,有 下 式 成 立 : 

limsup || x® —x* || * < 8+ Q 0 G^ 
其 中 4 =liminfl. 特别 地 ,如 果 limi, = 09 ,那么 欠 代 序列 的 收敛 率 为 尺 ERPE SA, E. R 
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收敛 因子 至 多 为 8, 即 limsup || x® —x* || *<p, 
证 明 : 由 Picard-MHSS 迭代 格式 (6. 31) ,经 计算 可 以 得 到 


ka 
xP — M ("xP + 9) M GGG) (9), k= 0.12. (6.32) 
其 中 M(a) Fil Ga) H C6. 29) EX. 由 于 x'€DX FGo-0 的 解 , 故 有 
x” =M "x" + 5 M (oiG(a) B(x"), bk — 0.1.2, (6. 33) 
式 (6. 32) 减 去 式 (6. 33) ,我 们 得 到 ” 


xD — x" = M a)" (x — x* 2 + b» M aG CDO (x?) — B(x") J, k—0.1.2.- 


j=0 
(6. 34) 
id 
EGx' 03,0) 6G 0-9'GoG- x») 
则 式 (6. 34) 可 以 等 价 地 写成 
xew y’ 
4a 4-1 
= [M co + >) M (GG Go] G9 一 xz) 十 > M GYGGOEG? x") 
j=0 j=0 


= [M (a)* + (I — MG) At $'G'*0] GÓ — x") + I Mla) ATEC” yx") 
= [AT Ø x") + Mta) I— A7 $'(xc )) (x — x") + OG — MG) DA? E” ,x* ) 


(6. 35) 
因此 
| x? —x* | «EI A? ,'G51 + IMG | Ft AT D ao’ DOJ —x + 
+ Il Mæ l^» | A? ll | EGO sx» ll 
= [BHA +AA] || x? — x* || 2 pO 2-0€0^). || EG? sx ) II (6. 36) 


由 函数 ® : DC CC" 在 x ED WET BIN, © DAG n] SRI "Ex" € D 处 连续 
可 以 得 出 对 于 任意 给 定 的 s 盖 0, 对 于 所 有 的 x € N, 都 存在 * 的 一 个 开 邻 域 NCN, ,使 得 
| EGx,x*) |] <e || x—x' || 


ie Ae ge] In{ =) tmr frm 
B-- OQ -- B0 Ca)^ <1 
进一步 假设 e 满足 
B+ A+8)0 (a) +p (14-0 (a) )e 1 

则 当 x, € N JFB lo — liminfZ, 时 ,根据 式 (6. 36) ,利用 数学 归纳 法 ,我 们 立刻 得 到 

Bw} <e | x*?-—x* || 

|| x? —x* | <[B+04+ 0 C^ +4146 Ca) de] || x'? —x* || (6.37) 
并 且 对 于 所 有 的 &=0,1,2,… 都 有 x' EN 成 立 . 由 式 (6. 36) 和 式 (6. 37) 可 以 得 到 


limx® = x 


keo 
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和 
limsup | x® —x* | * s g-- Q1 2- B6 Ca^» +p 2-0 (a) de 

由 se>0 的 任意 性 和 RR 收敛 因子 范 数 独 立 可 知 ,由 Picard-MHSS 迭代 所 产生 的 序列 x ag 
R 收敛 因子 至 多 为 (1 十 6)0 (a)^. 另外 ,车 liml 二 ,由 式 (6. 37) 可 知 ,序列 x 是 尺 线 
性 收敛 的 且 R 收敛 因子 至 多 为 B. 

定理 6.6 表明 Picard-MHSS 迭代 的 收敛 速率 实质 上 由 B 和 0(a) 决 定 . 通常 ,8 和 
0(a) 越 小 ,Picard-MHSS 和 迭代 的 收敛 速率 越 快 ， 当 弱 非 线性 方程 组 中 的 系数 矩阵 4 为 病 
态 时 ,0(a) 接 近 于 1,Picard-MHSS 迭代 的 收敛 速率 会 变 得 非常 慢 以 至 于 内 迭代 步 数 很 
多 . 在 这 种 情况 下 ,可 以 考虑 使 用 预 处 理 的 MHSS(CPMHSS) 和 迭代 代替 ,这 样 会 使 得 0(a) 
小 一 些 , 不 精确 的 Picard 迭代 的 收敛 速率 更 快 一 些 . 

本 节 最 后 ,为 了 方便 数值 计算 ,将 Picard-MHSS 和 迭代 写成 下 面 等 价 的 残 量 校正 形式 . 

Picard-MHSS 迭代 残 量 校正 形式 : HO. DCC" 一 C" 为 一 连续 可 微 函 数 ,矩阵 A= 
WHIT HEPER WT 41523 SziE E E [RUE 2E oe AB I FLU WT —W. T! — T. 给 定 初 始 
估计 x© € D WARE RUF Uo XIF 8—0.1.2.-- A F EMRA xd, 
到 近似 解 x** 达到 给 定 精度 . 

(D di 5*9 ,—9, p : =0 (x'? ) — Ax? 

© 对 于 (一 0,1,2,… 必 一 1, 求 解 下 述 线性 方程 组 得 到 sm s 


(al --W)s ^) = (al — iT)s ^? +b” 
(al + T)s^*? = (al +iW)s (t) — yo 


Q) 4 xc, —xQ sone, 
6.4.2 非 线 性 MHSS-like 和 迭代 方法 及 其 收敛 定理 


在 Picard-MHSS 算法 中 ,内 和 迭代 数 0.1.2. 通常 依赖 于 所 考虑 的 线性 方程 
组 ,在 数值 计算 中 很 难 决定 其 大 小 ,为 了 避免 使 用 显 式 的 内 迭代 过 程 ,但 是 仍然 保持 
Picard-MHSS 方法 的 优点 ,基于 非 线 性 固定 点 方程 
(GI +W)x = (aI —iDx-- d) 和 (al +T)x = (al +iW)x—i@ (x) 
我 们 提出 下 面 的 非 线性 MHSS-like 迭代 方法 . 
非 线 性 MHSS-like 迭代 方法 : 给 定 初始 估计 x" ED, 利 用 下 面 的 步骤 计算 x7 , 直 
到 近似 解 x? 达到 给 定 精 度 . 
(aI --W)x (2) = (af —iT)x® +6 (x?) 
(al +T)x* = (al --i:W)x 2) — i (x C22) 
这 里 ,a 为 给 定 的 常数 , 工 表示 单位 矩阵 . 下 面 给 出 非 线性 MHSS-like 迭代 方法 的 局 部 收 
SER 
定理 6.7. x HIEREN EHG. DARED: DOCCE x^ lib F nT 5t ABI A= 
W 十 并, 其 中 和 矩 阵 W,T 分 别 为 实 正定 矩阵 和 正 半 定 矩 阵 且 满足 W —W.T" — T. 
id M(a3x*) = (al+D GI iW — iB'(x* )) (GI HW) (al —iT+® (x’)) 


(6. 38) 


(6. 39) 
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如 果 pCMCGasx " ))<1, M] x* € D HARA HE MHSS-like 迭代 的 一 个 吸引 子 . 
证 明 : 首先 ,定义 函数 
U(x) = (al +W) (Cal — iT) x 4-6 (x)) 
V(x) = (GI +T) 1 CI + iW)x — i 6 (x)) 


(6.40) 


和 
Ya) =V °- U(x): = V(U(x)) 
则 非 线性 MHSS-like 迭代 格式 (6. 39) 等 价 于 
0 
由 Ostrowski 定理 可 知 , 如 果 pC’ Ce 00 —1.J A x^ 为 非 线性 MHSS-like 迭代 序列 的 
吸引 子 . 因此 ,下 面 我 们 证 明 CW Cr" )) 志 1. 因为 x "为 非 线性 方程 组 的 解 ,所 以 有 
x” =U(x"), x' =V(x") 


和 


U(x") = (GI +W) (Cal — iDx^ +® G5) 
Vix") = (I +T)? CaI +iW)x* — id (c 0) 
显然 成 立 . 利用 链 式 求 导 法 则 ,我 们 可 以 得 到 
W(x") =V'(x" U(x") = MGix^) 
由 此 ,定理 结果 显然 成 立 . 证 毕 . 
下 面 的 定理 给 出 了 pCM(a;x" <1 的 充分 条 件 . 
定理 6.8 假设 ® : DCC'—C"'il EEH 6.7 的 条 件 . 引入 下 列 记号 
à = max{ || Bx (al HW |, ll B(x" Ca +7)" l Oa) = | MCa) Ils 
其 中 矩阵 M(o) 如 (6.27) 所 定义 . 车 
0 一 一 1 十 V2 一 6Ca) (6.41) 
TW oCMCa ix * <1. Hl x* ED 为 非 线性 MHSS-like 和 迭代 序列 的 吸引 子 . 
证 明 : 经 过 计算 可 以 得 到 下 面 的 等 式 
(af + T)M(Ca ix^ (al + T)? (6. 42) 
= (al + T)M(a) laI +T) + (al + iW) (al 4-W)? $'(x' (al + T?! 一 
i ' (x' )(al +W)? (al — iT) (al 4- T)? —i®’(x* )(aI --W)? ®'(x" )(aI 4- T) 
两 边 同时 取 范 数 , 得 到 
| MCasx*) || = || (aI +T)M(a;x* (aI + T)? ||, 
x | (aI +T)M (a) (aI +T)7 |lz || (aI +iW) (aI +W) $' (x^ ) (aI - T)? llo + 
| iØ’ (x* ) (eI À-W)? (aI —iT) (aI + T)? |; 
| i^ (x^ ) (a1 --W)? ®’ (x" ) (aI -- T)? ll; 
« Il M (a) ll + | (aI +iW) (eI --W)? lla || ©’ (x^ ) (aI - T)? l2 十 
| ^ (x* ) (aI -W)? ls || («E — iT) (aI - T) la 
| ^ (x^ ) (ap -W)? ls || D’ (x^ ) (er T)? Ils 
x) + 26+ 8 (6. 43) 
事实 上 ,在 上 面 的 证 明 过 程 中 我 们 利用 了 
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l| CI +W) (aI --W)! |; — 1. |l (eE — iT) (aI - T)? ll; — 1. 
利用 定理 中 条 件 (6. 4D ,我 们 有 
Wa) +26+8 <1 
因此 ,我们 有 
e(M(a;x")) x || M (aix^ ) || <1 
x” ED 为 非 线 性 MHSS-like 迭代 序列 的 吸引 子 . 


6.4.3 数值 结果 
为 了 验证 Picard-MHSS 方法 的 可 行 性 与 有 效 性 ,我 们 考虑 下 面 的 二 维 非 线性 方程 
— (a o ii) (Qus + us, ) + pu = (a2 Fig )ue" 4 sin V1 +u +u, (Gr. EN 
u(xsy) = 0, (x.y) € 20 
其 中 A=(0,1)X 0.1), 20 为 区 域 2 的 边界 ,po 为 正常 数 . MAM a =p — 1.0; — [5 — 


0.5. 在 等 分 布 网 格 上 应 用 中 心 差分 格式 h— 二, 可 以 得 到 下 面 的 非 线性 方程 组 


F(x) =Ax—@(x)=0, 其 中 n=NXN 
在 我 们 的 计算 中 ,选取 x — 0. 9b Picard 迭代 停止 标准 取 为 


| FG?» |l; 
| Fo) |, 


内 HSS 和 MHSS 和 迭代 的 停止 标准 都 取 为 

l| F'(x'?5s*"? + Fc 

| FG) |, 
由 于 Picard-HSS 与 Newton-HSS, Newton-GMRES 等 方法 的 效率 已 经 在 相关 文献 

中 进行 了 详细 的 对 比 , 此 处 只 给 出 Picard-MHSS Picard-HSS, dEZE E HSS-like 和 非 线性 
MHSS-like 四 种 方法 的 效率 对 比 结果 . 在 内 迭代 MHSS 和 HSS 中 ,我 们 利用 LU 分 解 来 
求解 方法 中 的 子 线性 系统 . 除 此 之 外 ,选取 六 为 常数 , 记 为 x. 选择 不 同 的 nn 值 ,不 同 的 p 
值 以 及 不 同 的 内 迭代 停止 标准 ,利用 4 种 方法 分 别 求解 非 线性 方程 组 . 表 6. 9 和 表 6. 10 
分 别 给 出 了 Picard-MHSS 和 Picard-HSS 两 种 方法 中 最 优 参数 的 值 ,从 表格 中 可 以 看 出 ， 
两 种 方法 中 最 优 参数 的 值 完 全 一 致 . 

表 6.9 Picard-MHSS 方法 的 最 优 参数 c 


= o 


a 
) 
l c 


N p=1 e=10 0-50 0-100 
10 1.1559 1. 3921 2.1789 2.9265 
15 0. 8004 0. 9607 1. 4865 1.9714 
20 0.6112 0. 7330 1.1293 1. 4903 
25 0. 4924 0. 5924 0. 9108 1.1991 
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表 6.10  Picard-HSS 方法 的 最 优 参数 c 
N pal p=10 p=50 p= 100 
10 1.1559 1.3921 2.1789 2.9265 
15 0. 8004 0. 9607 1. 4865 1.9714 
20 0.6112 0. 7330 1.1293 1. 4903 
25 0. 4924 0. 5924 0. 9108 1.1991 


K 6. 11 和 表 6.13 分 别 给 出 了 4 BEATE PEAR AG PC Tina Sb EAR AE PC IT 的 比较 ， 
表 6.12 和 表 6. 14 给 出 了 4 种 方法 CPU 时 间 的 比较 . 从 表格 中 可 以 看 出 ,与 Picard-HSS 
和 非 线 性 HSS-like 方 法 相 比 ,Picard-MHSS 和 非 线 性 MHSS-like 迭代 在 内 、 外 迭代 数目 
上 相差 不 多 ,但 是 在 迭代 时 间 上 随 着 求解 规模 n 的 不 断 增加 ,Picard-MHSS 和 非 线 性 


MHSS-like 的 求解 时 间 有 了 很 大 程度 的 减少 . 


ROU 当 %==0.1 时 4 种 方法 迭代 步 数 的 比较 
N 5 10 12 15 
非 线 性 HSS-like IT 25 46 55 70 
非 线性 MHSS-like IT 30 49 56 68 
IT. 24 41 45 57 
e=1 | Picard-HSS 
IT. 6 7 7 8 
ITine 31 48 55 64 
Picard-MHSS 
ITA, 6 f 7 8 
非 线性 HSS-like iT 25 48 43 54 
非 线 性 MHSS-like IT 30 36 55 65 
ITa 20 34 38 56 
p=10 | Picard-HSS 
Wa 6 j 7 7 
It. 31 47 39 53 
Picard-MHSS 
ea 6 7 7 7 
R612 Hq =0.1 时 4 种 方法 迭代 时 间 的 比较 
N 5 10 12 15 
非 线 性 HSS-like 0. 1092 4.7112 13. 400 49. 874 
非 线 性 MHSS-like 0. 1404 3. 8680 10. 639 35. 116 
p=1 
Picard-HSS 1. 2324 28. 267 74. 865 318. 68 
Picard-MHSS 0. 7020 17.597 57.05 258. 32 
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BR 
N 5 10 12 15 
非 线性 HSS-like 0. 1092 3. 8220 10. 546 39. 094 
非 线性 MHSS-like 0. 1248 2. 5100 8. 1560 33. 431 
Picard-HSS 1.0140 23. 931 62. 619 250. 68 
Picard-MHSS 0.7020 13.712 46. 301 197. 69 
R613 Aq =0.01 时 4 种 方法 迭代 步 数 的 比较 
N 5 10 12 15 
非 线性 HSS-like IT 25 46 55 70 
非 线性 MHSS-like IT 30 49 56 68 
Mis 28 47 56 67 
p=1 |Picard-HSS 
if. 6 7 8 8 
Ii 36 55 63 75 
Picard-MHSS 
Wi 6 7 7 7 
非 线性 HSS-like IT 25 48 43 54 
非 线性 MHSS-like IT 30 36 55 65 
ITiw 24 40 46 58 
0710 | Picard-HSS 
Ta 6 7 7 7 
ITiw 36 55 63 66 
Picard-MHSS 
Ir 6 y 7 7 
56.14. 349 —0.01 时 4 种 方法 迭代 时 间 的 比较 
N 5 10 12 15 
非 线性 HSS-like 0. 1092 4. 7112 13. 400 49. 874 
非 线性 MHSS-like 0. 1404 3. 8680 10. 639 35.116 
da Picard-HSS 1. 2948 32. 339 102. 02 381. 98 
Picard-MHSS 0. 8424 28.077 76.968 300. 145 
非 线性 HSS-like 0. 1092 3. 8220 10.546 39. 094 
非 线性 MHSS-like 0. 1248 2.5100 8. 1560 33. 431 
ate Picard-HSS 1. 0608 26.599 89. 270 279. 76 
Picard-MHSS 0. 8241 22.718 78. 440 214, 09 
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